arXiv:math/0308086v1 [math.CA] 8 Aug 2003 in uwt eafnto,Siln ubr,hroi nu harmonic numbers, Stirling function, zeta constant. Hurwitz tion, Foundation. hr Γ( where ouini eadtecniino ovxt [32] convexity of condition the add we if solution 11Y. ue n tde se[,1,1,1] eeaiaino th of functio generalization gamma a multiple the 12]) called 11, function, 10, gamma [9, Euler (see cal studied and duced (1) 32]: [31, tion h ucinΓ function The hswr a upre ygatCR0003fo h Nation the from CCR-0204003 grant by supported was work This 1991 e od n phrases. and words Key nasqec fppr ulse ewe 8910,Barnes 1899-1904, between published papers of sequence a In OTIUIN OTETER FTHE OF THEORY THE TO CONTRIBUTIONS am ucinadGase’ osataeas discussed. also Ba are the constant Glaisher’s of and computation function high-precision gamma functio for gamma schemes multiple the numerical of series asymptotic and tations Abstract. ahmtc ujc Classification. Subject Mathematics z steElrgmafnto.Tesse 1 a unique a has (1) system The function. gamma Euler the is ) Γ Γ Γ n 1 n n ( ( ( +1 1 1 = (1) − z z hspprpeet aiyo e nerlrepresen- new of family a presents paper This aifistefloigrcrec-ucinlequa- recurrence-functional following the satisfies ) Γ( = ) ( 1) z n )= 1) + ANSFUNCTION BARNES +1 ansfnto,gmafnto,Reanzt func- zeta Riemann function, , Barnes , z d z d 1. ) n .S ADAMCHIK S. V. , n +1 +1 Introduction Γ Γ n o Γ log +1 n ( z ( 1 z ) ) n z , ( z rmr 32,1M5 Secondary 11M35; 33E20, Primary ) ≥ ∈ 0 C > z , n , ∈ 0 . N br,Glaisher’s mbers, , .The n. rnes lScience al Γ n classi- e intro- n ( z ). 2 V.S.ADAMCHIK

In this paper we aim at the most interesting case G(z)=1/Γ2(z), being the so-called double gamma function or the Barnes G-function: G(z +1)= Γ(z) G(z), z C, ∈ G(1) = G(2) = G(3) = 1,

d3 log G(z) 0, z> 0. d z3 ≥ The G-function has several equivalent forms including the Weier- strass canonical product

2 ∞ 2 z z + z (1 + γ) z k z G(z +1) =(2π) 2 exp 1+ exp z − 2 k 2k −   kY=1     and the functional relationship with the ζ(t, z), due to Gosper [20] and Vardi [31]

(2) log G(z + 1) z log Γ(z)= ζ′( 1) ζ′( 1, z), (z) > 0, − − − − ℜ where γ denotes the Euler-Mascheroni constant and the derivative ′ d ζ (t, z) = dt ζ(t, z) is defined as an provided by the Hermite integral (18):

z2 z2 z ζ′( 1, z)= log z log z + − 2 − 4 − 2 (3) ∞ arctan( x ) ∞ x log(x2 + z2) 2 z z dx + dx. (z) > 0, e2πx 1 e2πx 1 ℜ Z0 − Z0 − In particular,

∞ x log x (4) ζ′( 1) = ζ′( 1, 0)=2 dx. − − e2πx 1 Z0 − For z , the G-function has the Stirling asymptotic expansion →∞ z 3z2 z2 1 1 log G(z +1)= log(2π) log A + log z + O , 2 − − 4 2 − 12 z     where A is the Glaisher-Kinkelin constant given by 1 (5) log A = ζ′( 1). 12 − − CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 3

Combining this definition with the representation (4), we obtain a new integral for log A:

1 + log(2π) 1 ∞ x log x (6) log A = dx. 12 − 2π2 ex 1 Z0 − The constant A originally appeared in papers by Kinkelin [23] and Glaisher [17, 18, 19] on the asymptotic expansion (when n ) of the following product → ∞

p p p 11 22 ...nn , p N. ∈ Similar to the Euler gamma function, the G-function satisfies the mul- tiplication formula [31]:

′ 2 2 2 G(n z)= eζ (−1)(1−n )nn z /2−nz+5/12 n−1 n−1 (7) n−1 i + j (2π) 2 (1−nz) G z + , n N n ∈ i=0 j=0 Y Y   The Barnes function is not listed in the tables of the most well-known special functions. However, it is cited in the exercises by Whittaker and Watson [33], and in entries 6.441, 8.333 by Gradshteyn and Ryzhik [21]. Recently, Vardi [31], Quine and Choi [27], and Kumagai [24] computed the functional determinants of Laplacians of the n-sphere in terms of the Γn-function. Choi, Srivastava [13, 14] and Adamchik [3, 16] presented relationships (including and series) between other mathematical functions and the multiple gamma function.

The Γn function has several interesting applications in pure and ap- plied mathematics and theoretical physics. The most intriguing one is the application of G-function to the Riemann Hypothesis. Mont- gomery [26] and Sarnak [29] (see [22] for additional references) have conjectured that the non-trivial zeros of the are pairwise distributed like eigenvalues of matrices in the Gaussian Uni- tary Ensemble (GUE) of random matrix theory. It has been shown in works by Mehta, Sarnak, Conrey, Keating, and Snaith that a closed representation for statistical averages over GUE of N N unitary ma- trices, when N can be expressed in terms of the Barnes× functions. →∞ Another interesting appearance of the G-function is in its relation to the determinants of the Hankel matrices. Consider a matrix M of the Bell numbers: 4 V.S.ADAMCHIK

B1 B2 ... Bn B B ... B M = 2 3 n+1 n  ...   Bn Bn+1 ... B2n−1    where Bn are the Bell numbers, defined by 

n n B = n k k=0   n X and k are the Stirling subset numbers [28]  n n 1 n 1 n 1, n =0, (8) = k − + − , = k k k 1 0 0, n =0      −     6 It was proved in [5], that

det M = G(n + 1), n N. n ∈

2. Special cases of the Barnes Function For some particular arguments, the Barnes function can be expressed in terms of the known special functions and constants. The simplest (but not-trivial) special case is due to Barnes [9]:

1 1 log2 log π 3 (9) logG = + log A, 2 8 24 − 4 − 2   where A is the Glaisher-Kinkelin constant (5). The next two special cases are due to Srivastava and Choi [14]:

1 3 G 3 1 9 (10) log G = logΓ log A, 4 32 − 4π − 4 4 − 8     3 1 G log2 log π 1 (11) log G = logG + +logΓ , 4 4 2π − 8 − 4 4 where G denotes  Catalan’s constant:  

∞ ( 1)k G = − . (2k + 1)2 Xk=0 Identities (10) and (11) are proven by using the duplication formula 1 1 2 (7) (setting n = 2) together with (9). More special cases for z = 6 , 3 , 3 CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 5

5 and 6 can be derived as well. Here is one of such formulas (see [3] for the others): Proposition 2.1. 1 1 log3 π log G = + + 3 9 72 18√3 − (12)   2 1 4 ψ(1) 1 logΓ log A 3 , 3 3 − 3 − 12√3π    2 (1) ∂ where ψ (z)= ∂z2 log Γ(z) is the .

Proof. We recall the Lerch functional equation for the Hurwitz zeta function (see [2, 8])

πi ζ′( n, z)+( 1)nζ′( n, 1 z)= B (z)+ − − − − n +1 n+1 n! e−πin/2 Li e2πiz , 0

G(1 + z) sin(πz) 1 (13) log = z log Cl (2 πz), 0

(14) Cl (z)= (Li (e−iz)). 2 −ℑ 2 1 Upon setting z = 3 to (13), and making use of

(1) 1 2 2π ψ 3 2π Cl2 = , 3 3√3 − 9√3    we obtain

G 1 √3 2π 2π2 3ψ(1) 1 (15) 3 = exp − 3 . G 2 √6 3Γ 1 18π√3 3  3     6 V.S.ADAMCHIK

1 On the other hand, using the multiplication formula (7) with z = 3 and n = 3, we find

7/12 2/3 1 2 3 3 e (16) G G = 1 . 3 3 s2πA8Γ     3 Now, combining (15) and (16), we conclude the  proof.  It remains to be seen if the Barnes function (or the multiple gamma function) can be expressed as a finite combination of elementary func- tions and constants for other rational arguments. The most general result available for the multiple gamma function Γn(z) is a closed form 1 at z = 2 : Proposition 2.2. For n N ∈ 1 (2n 3)!! log π ( 1)n(n 1)!logΓ ( )= − + − − n 2 − 2n (17) n 1 n k Pk,n( ) Bk+1 2 1 1 log2 2 + − P ( )ζ′( k), (k +1)2k 2k k,n 2 − Xk=1 Xk=1 where P ( 1 ) are coefficients by xk, k n in expansion of k,n 2 ≤ n−1 1 (x + j ) − 2 j=1 Y Proof. The proof follows directly from [4], formula (17), by setting 1  z = 2 . Here are a few particular cases: the gamma function (n = 1): 1 Γ( )= √π 2 the Barnes function (n = 2):

1 21/24 e1/8 G( )= 2 A3/2 π1/4 the triple gamma function (n = 3):

1 A3/2 π3/16 7 ζ(3) 1 Γ ( )= exp 3 2 21/24 32 π2 − 8   CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 7

3. Hermite integrals In this section we derive a few integral representations for the Barnes function. We start with recalling the Hermite integral (see [7, 25]) for the Hurwitz zeta function:

z−s z1−s ∞ sin(s arctan(x/z)) (18) ζ(s, z)= + +2 s/2 dx, 2 s 1 0 (x2 + z2) (e2πx 1) − Z − which provides analytic continuation of ζ(s, z) to the domain s C 1 . Differentiating both sides of (18) with respect to s, letting s∈= −1 {(and} s = 2), and using the second Binet formula for log Γ(z) (see [7])− − 1 log(2π) ∞ arctan(x/z) (19) logΓ(z)=(z )log z z + +2 dx − 2 − 2 e2πx 1 Z0 − together with (2), we readily obtain:

z2 log G(z +1)= log z H zζ′(0) ζ′( 1) 2 − 2 − − − − (20)     ∞ x log(x2 + z2) dx, (z) > 0, e2πx 1 ℜ Z0 − p ′ where Hp = k=1 1/k are the harmonic numbers and ζ (z) is a deriv- ative of the Riemann zeta function. P In the same manner, repeating the above steps n times, we derive the general integral representations.

Proposition 3.1. For (z) > 0 ℜ

(2n)! logΓ2n+1(z +1)=

z2n+1(log z H ) 2n 2n − 2n+1 ( 1)k ζ′( k) z2n−k 2n +1 − − k − − Xk=0   (21) 2n−1 2n ( 1)kk! logΓ (z +1)+ − k k+1 Xk=1   ∞ x2n arctan x/z 2( 1)n dx − e2πx 1 Z0 −  8 V.S.ADAMCHIK and

z2n(log z H ) (2n 1)!logΓ (z +1) = − 2n + − 2n − 2n

2n−1 2n 1 ( 1)k − ζ′( k)z2n−k−1 + − k − Xk=0   (22) 2n−2 2n 1 ( 1)kk! − logΓ (z + 1) − k k+1 − Xk=1   ∞ x2n−1 log(x2 + z2) ( 1)n dx − e2πx 1 Z0 − n where k are the Stirling subset numbers, defined in (8).  Upon substituting z = 1 into formulas (21) and (22), we obtain new integrals

∞ 2n 2n n x arctan(x) H2n+1 k 2n ′ (23) 2( 1) 2πx dx = + ( 1) ζ ( k) − 0 e 1 2n +1 − k − Z − Xk=0   (24) ∞ 2n−1 2 2n−1 n x log(1 + x ) H2n k 2n 1 ′ ( 1) 2πx dx = + ( 1) − ζ ( k) − 0 e 1 2n − k − Z − Xk=0   Performing simple transformations of the integrands in (23) and (24) will apparently lead to a large number of special integrals:

∞ x log x 2 dx = ζ′( 1) e2πx 1 − Z0 − ∞ x log x dx = 12ζ′( 1)+log2 e2πx +1 − Z0 ∞ x log(1 + x2) 3 23 1 dx = log2 + ζ′( 1) e2πx +1 4 − 24 2 − Z0 ∞ arctan x log(2π) 1 2 z dx = z log z z + logΓ z + , (z) > 0 e2πx +1 − 2 − 2 ℜ Z0   CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 9

∞ arctan(x) 4 dx = 2+3log2 e2πx +1 − Z0 ∞ x dx 1 π 2 = ψ z + log z, arg(z) = (e2πx + 1)(x2 + z2) 2 − 6 2 Z0  ∞ x dx 1 π 2 = log z ψ(z) , arg(z) = (e2πx 1)(x2 + z2) − − 2z 6 2 Z0 −

4. Binet-like representation In this section we derive the Binet integral representation for the Barnes function. Proposition 4.1. The Barnes G-function admits the Binet integral representation:

z2 log z log G(z +1)= z log Γ(z)+ B (z) log A 4 − 2 2 − − (25) ∞ e−zx 1 1 1 x 2 −x dx, (z) > 0. 0 x 1 e − x − 2 − 12 ℜ Z  −  Proof. Recall the well-known integral for the Hurwitz function [8]

1 ∞ xs−1e−zx (26) ζ(s, z)= dx, (s) > 1, (z) > 0 Γ(s) 1 e−x ℜ ℜ Z0 − The integral in the right hand-side of (26) can be analytically continued to the larger domain (s) > 2 by subtracting the truncated Taylor 1 ℜ − series of 1−e−x at x = 0. Since

1 1 1 x (27) = + + + O(x3) 1 e−x x 2 12 − and

1 ∞ 1 1 x s z−s z1−s xs−1e−zx + + dx = z−1−s + , Γ(s) 0 x 2 12 12 2 − 1 s Z   − (s) > 1, (z) > 0 ℜ ℜ we therefore can rewrite (26) as 10 V.S.ADAMCHIK

z1−s z−s s ζ(s, z)= + + z−s−1 + −1 s 2 12 − ∞ 1 s−1 −zx 1 1 1 x x e −x dx, (s) > 2. Γ(s) 0 1 e − x − 2 − 12 ℜ − Z  −  Differentiating the above formula with respect to s and computing the limit at s = 1, we obtain − 1 z2 log z ζ′( 1, z)= + B (z)+ − 12 − 4 2 2 ∞ e−zx 1 1 1 x 2 −x dx, (z) > 0, 0 x 1 e − x − 2 − 12 ℜ Z  −  where B2(z) is the second Bernoulli polynomial. The convergence of the above integral is ensured by (27). Finally, using (2), we conclude the proof. 

5. Asymptotic expansion of Γn(z)

The asymptotic expansion for Γn(z) when z follows straight- forwardly from formulas (21) and (22) by expanding→ ∞ log(1 + x2) and arctan(x) into the Taylor series, and then performing formal term-by- term integration. Lemma 5.1. For (z) > 0 ℜ z2 3 log z log G(z +1)= log z zζ′(0) + ζ′( 1) 2 − 2 − 12 − − −   (28) n B 1 2k+2 + O( ), z . 4 k (k + 1) z2k z2n+2 →∞ Xk=1 Proof. We start with the integral representation (20). By recalling the Taylor series for log(1 + x)

n x ( 1)k x k+1 ( 1)n x yn+1 log(1 + )= − ( ) −n+1 dy, z k +1 z − z 0 y + z Xk=0 Z we find that CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 11

∞ x log(1 + x2/z2) dx = e2 πx 1 Z0 − 2 n ( 1)k 1 ∞ x2k+3 ( 1)n ∞ t dt t yn+1 − 2k+2 2 πx dx −2n+2 2 πt 2 dy k +1 z 0 e 1 − z 0 e 1 0 y + z Xk=0 Z − Z − Z which, by appealing to (formula 1.6.4 in [7])

∞ t2k−1 B (29) dt =( 1)k+1 2k , k N, e2 π t 1 − 4 k ∈ Z0 − reduces to

∞ x log(1 + x2/z2) dx = e2 πx 1 Z0 − n+1 ∞ t/z 2n+3 B2k+2 n t dt y 2k 2( 1) 2 π t 2 dy − 4 k (k + 1)z − − 0 e 1 0 1+ y Xk=1 Z − Z

Substituting this into (20), we complete the proof. 

Lemma 5.2. For (z) > 0 ℜ

z3 11 z2 3 2logΓ (z +1) = log z + log z + ζ′(0) 3 − 3 − 6 2 − 2 −     log z z ζ′(0)+2ζ′( 1) + ζ′( 1) + ζ′( 2) 12 − − − − − (30)   n B n B 2k+2 2k+2 + 4 k (k + 1) z2k − 2(k + 1)(2k 1) z2k−1 Xk=1 Xk=1 − 1 O( ), z . z2n+1 →∞ 12 V.S.ADAMCHIK

Proof. Employing the same technique as in getting (28), we begin with the integral representation of the triple gamma function:

z3 2logΓ (z +1)= (log z H )+log G(z +1)+ 3 − 2 − 3 (31) ∞ 2 2 ′ ′ ′ x arctan(x/z) z ζ (0) 2zζ ( 1) + ζ ( 2) +2 2 πx dx − − − 0 e 1   Z − that follows directly from Proposition 3.1 with n = 1. Expanding arctan into the Taylor series

n x ( 1)k−1 x 2k−1 ( 1)n x y2n arctan( )= − ( ) + −2n−1 2 2 dy z 2k 1 z z 0 y + z Xk=1 − Z and making use of (29), we readily obtain

∞ 2 n x arctan(x/z) B2k+2 2 πx dx = 2k−1 + 0 e 1 − 4(k + 1)(2k 1) z Z − Xk=1 − ( 1)n ∞ t2 dt t y2n − dy. z2n−1 e2 π t 1 y2 + z2 Z0 − Z0 The conclusion of the lemma follows immediately by substituting this into formula (31) and using Lemma 5.1. 

6. Implementation In this section we discuss a numeric computational scheme for the Barnes function. Let us consider an integral representation dated back to Alexejewsky [6] and Barnes [9]:

z z(1 z) z log Γ(x)dx = − + log(2π)+ z log Γ(z) log G(z + 1), 2 2 − Z0 where (z) > 1. Employing integration by parts, we derive ℜ − z(1 z) z z (32) log G(z +1) = − + log(2π)+ x ψ(x) dx, (z) > 1 2 2 ℜ − Z0 The integral representation (32) provides an efficient numeric procedure for arbitrary precision evaluation of the G-function. This representa- tion demonstrates that the complexity of computing G(z) depends at CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 13 most on the computational complexity of the polygamma function. The later can be numerically computed by the Spouge approximation [30] (a modification of the Lanczos approximation). For the numerical integration in (32) we shall use the Gaussian quadrature scheme.

The restriction (z) > 1 in (32) can be easily removed by analyt- icity of the polygamma.ℜ For− example, by resolving the singularity of the integrand at the pole x = 1, we continue log G(z +1) to the wider area (z) > 2: − ℜ −

z(1 z) z z 1 log G(z+1) = − + log(2π)+log(z+1)+ xψ(x) dx 2 2 − x +1 Z0   Another method of the analytic continuation of the G-function is to use the following well-known identity

1 ψ(x)= ψ( x) π cot(πx) − − − x which, upon substituting it into (32), yields

z (33) log G(1 z) = log G(1 + z) z log(2π)+ πx cot(πx) dx. − − Z0 The identity (33) holds everywhere in a complex plane of z, except the real axes, where the integrand has simple poles. Therefore, combining (32) with (33) yields the following definition of the Barnes function

G( n)=0, n N, − ∈ (z 1)(z 2) (34) G(z)=(2π)(z−1)/2 exp − − + x ψ(x) dx , − 2  Z γ  arg(z) = π, 6 where the contour of integration γ is a line between 0 and z 1 that does not cross the negative real axis; for example, γ could be the− following path 0, i, i + z 1, z 1 . { − − } 14 V.S.ADAMCHIK

It remains to define G(z) for z R−. We do this by using the reflexion formula (13), which upon periodicity∈ of the Clausen function (14) can be rewritten as follows:

sin(πz) z+1 G( z)=( 1)⌊z/2⌋−1G(z + 2) − − π ∗ (35) 1 exp Cl 2π(z z ) , z R−. 2π 2 − ⌊ ⌋ ∈  

7. Glaisher’s constant In this section we discuss numeric computational schemes for Glaisher’s constant:

1 (36) A = exp ζ′( 1) 12 − − In light of the functional equation for the zeta function πs ζ(s)=2sπs−1Γ(1 s)ζ(1 s) sin − − 2 we can rewrite (36) as  

γ ζ′(2) (37) A = exp (2π)1/12. 12 − 2π2 Exactly this formula is used inMathematica V4.2 for computing Glaisher’s constant. Another computational scheme for Glaisher’s constant is based on the Barnes G-function identity:

1 log2 log π 2 1 (38) log A = + log G . 12 36 − 6 − 3 2 1   Here G( 2 ) can be computed via integral (32). A more efficient algorithm is provided by the following identity

log2 1 ∞ 3 6 (39) log A = + ζ(2k + 1) 1 28 + , 12 36 − 1+ k − 2+ k Xk=1    which allows to approximate Glaisher’s constant as CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 15

log2 1 N 3 6 1 log A = + ζ(2k + 1) 1 28+ + O 12 36 − 1+ k − 2+ k 4N Xk=1      This approximation is suitable for arbitrary precision computation - it requires p log 10 terms to achieve p decimal-digit accuracy. ⌈ 2 2 ⌉ We close this section with a proof of formula (39). Recall that [8]

1 ∞ xs−1 (40) ζ(s) 1= dx, (s) > 1. − Γ(s) ex(ex 1) ℜ Z0 − Upon replacing ζ(2k +1) 1 in (39) by its integral representation (40) and inverting the order of− summation and integration, we obtain

log2 1 ∞ dx ∞ t2k 3 6 (41) log A = + x x 28+ . 12 36 0 e (e 1) (2k)! 1+ k − 2+ k Z − Xk=1   Next, noting that ∞ t2k = cosh(t) 1 (2k)! − Xk=1 and

∞ 2k t t 2 2s−1 = 2s x (cosh(x) 1) dx, (2k)! (k + s) t 0 − Xk=1 Z we compute the inner sum

t4 ∞ t2k 3 6 28 + = 36+3 t2 14 t4 + 2 (2k)! 1+ k − 2+ k − − Xk=1   36+15 t2 + 14 t4 cosh(t) 36 t +3 t3 sinh(t) −   In the next step, we substitute this back to (41) and, employed by (26), evaluate the integral. Unfortunately, integration cannot be done term by term. The way around is to multiply the integrand by tλ,λ> 0 and consider the limit as λ 0. This yields → 16 V.S.ADAMCHIK

∞ 36 tλ−4 36 tλ−4 72 tλ−4 36 tλ−3 + + et 1 e2 t (et 1) − et (et 1) − et 1 Z0 − − − − 36 tλ−3 15 tλ−2 15 tλ−2 6 tλ−2 3 tλ−1 + + + + e2 t (et 1) et 1 e2 t (et 1) et (et 1) − et 1 − − − − − 3 tλ−1 14 tλ 14 tλ 28 tλ + + dt = e2 t (et 1) et 1 e2 t (et 1) − et (et 1) − − − − ! 36Γ(λ 3) 925−λ Γ(λ 3) 36Γ(λ 2) 924−λ Γ(λ 2) − − − − − − − − 3Γ(λ) 21Γ(λ 1) 1521−λ Γ(λ 1) 3Γ(λ) + 14Γ(λ + 1) − − − − − 2λ − 7Γ(λ + 1) + 36Γ(λ 1) ζ(λ 1). 2λ − − Computing the limit as λ 0, we find that the right hand-side of the above expression evaluates→ to

3 3log2 36 ζ′( 1). − − − This completes the proof of (39). CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 17

References [1] V. S. Adamchik, H. M. Srivastava, Some series of the zeta and related functions, Analysis (1998) 131–144. [2] V. S. Adamchik, Polygamma functions of negative order, J. Comp. and Appl. Math. 100 (1998) 191–199. [3] V. S. Adamchik, On the Barnes function, in Proceedings of the 2001 Interna- tional Symposium on Symbolic and Algebraic Computation, London, Ontario; July 22-25, 2001, ACM Press, New York, 2001, 15–20. [4] V. S. Adamchik, Multiple Gamma Function and Its Application to Computa- tion of Series, preprint, 2003, http://lanl.arxiv.org/abs/math.CA/0308074. [5] M. Aigner, A characterization of the Bell numbers, Discrete Math. 205 (1999) 207–210. [6] W. P. Alexejewsky, Ueber eine Classe von Funktionen, die der Gammafunktion analog sind, Leipzig Weidmanncshe Buchhandluns 46 (1894) 268–275. [7] G. Andrews, R. Askey, R. Roy, Special Functions, (Cambridge University Press, Cambridge, 1999). [8] H. Bateman, A. Erdelyi, Higher Transcendental Functions, Vol. 1, (McGraw- Hill, New York, Toronto, London, 1953). [9] E. W. Barnes, The theory of the G-finction, Quart. J. Math. 31 (1899) 264– 314. [10] E. W. Barnes, Genesis of the double gamma function, Proc. London Math. Soc. 31 (1900) 358–381. [11] E. W. Barnes, The theory of the double gamma function, Philos. Trans. Roy. Soc. London, Ser. A 196 (1901) 265–388. [12] E. W. Barnes, On the theory of the multiple gamma function, Trans. Cam- bridge Philos. Soc. 19 (1904) 374–425. [13] J. Choi, H. M. Srivastava, J. R. Quine, Some series involving the zeta function, Bull. Austral. Math. Soc. 51 (1995) 383–393. [14] J. Choi, H. M. Srivastava, Certain classes of series involving the zeta function, J. Math. Anal. and Appl. 231 (1999) 91–117. [15] J. Choi, H. M. Srivastava, Certain classes of series associated with the Zeta function and the multiple Gamma function, J. Comput. Appl. Math. 118 (2000) 87–109. [16] J. Choi, H. M. Srivastava, V. S. Adamchik, Multiple Gamma and related functions, Appl. Math. and Comp. 134 (2003) 515–533. [17] J. W. L. Glaisher, On a numerical continued product, Messenger of Math. 6 (1877) 71–76. [18] J. W. L. Glaisher, On the product 11, 22, 33,...,nn, Messenger of Math. 7 (1887) 43–47. [19] J. W. L. Glaisher, On certain numerical products, Messenger of Math. 23 (1893) 145–175. m 6 [20] R. W. Gosper, n log Γ(z) dz, in: M. Ismail, D. Masson, M. Rahman editors, 4 Special functions, q-series and related topics. The Fields Institute Communi- R cations, (Amer. Math. Soc., 14, 1997) 71–76. [21] I. S. Gradshteyn, I. M. Ryzhik, Table of Integrals, Series, and Products, Sixth Ed. (Academic Press, San Diego, 2000). 18 V.S.ADAMCHIK

1 [22] J. P. Keating, N. C. Snaith, Random matrix theory and ζ 2 + it , Commun. Math. Phys. 214 (2000) 57–89.   [23] Kinkelin, Ueber eine mit der Gammafunction verwandte Transcendente und deren Anwendung auf die Integralrechnung, J. Reine Angew. Math. 57 (1860) 122–158. [24] H. Kumagai, The determinant of the Laplacian on the n-sphere, Acta Arith- metica XCI.3 (1999) 199–208. [25] W. Magnus, F. Oberhettinger, R. P. Soni, Formulas and theorems for the special functions of mathematical physics (Springer, Berlin, 1966). [26] H. L. Montgomery, The pair correlation of the zeta function, Proc. Symp. Pure Math. 24 (1973) 181–193. [27] J. R. Quine, J. Choi, Zeta regularized products and functional determinants on spheres, Rocky Mountain Journal 26 (1996)719–729. [28] K. N. Rosen and et.., Handbook of Discrete and Combinatorial Mathematics (CRC Press, New York, 2000). [29] P. Sarnak, Quantum chaos, symmetry and zeta functions, Curr. Dev. Math. (1997) 84–115. [30] J. L. Spouge, Computation of the gamma, digamma, and trigamma functions, SIAM J. Numer. Anal. (1994)931–944. [31] I. Vardi, Determinants of Laplacians and multipe gamma functions, SIAM J. Math. Anal. 19 (1988)493–507. [32] M. F. Vign´eras, L‘´equation fonctionelle de la fonction zˆeta de Selberg du groupe mudulaire SL(2, Z), Ast´erisque 61 (1979)235–249. [33] E. T. Whittaker, G. N. Watson, A course of Modern Analysis: An Introduction to the General Theory of Infinite Processes and of Analytic Functions; With an Account of the Principal Transcendental Functions, Fourth Ed. (Cambridge University Press, Cambridge, 1963).

Carnegie Mellon University, Pittsburgh, USA E-mail address: [email protected]