arxiv:math/886v [math.ca] 8 Aug CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION V. S. ADAMCHIK Abstract. This paper presents a family of new integral representations and asymptotic series of the multiple gamma function. The numerical schemes for high-precision computation of the Barnes gamma function and Glaisher s constant are also discussed.. Introduction In a sequence of papers published between 899-94, Barnes introduced and studied see [9,,, ] a generalization of the classical Euler gamma function, called the multiple gamma function Γ n z. The function Γ n z satisfies the following recurrence-functional equation [, ]: Γ n+ z + = Γ n+z, z C, n N, Γ n z Γ z = Γz, Γ n =, where Γz is the Euler gamma function. The system has a unique solution if we add the condition of convexity [] n+ dn+ d z log Γ nz, z >. n+ 99 Mathematics Subject Classification. Primary E, M5; Secondary Y. Key words and phrases. Barnes function, gamma function, Riemann zeta function, Hurwitz zeta function, Stirling numbers, harmonic numbers, Glaisher s constant. This work was supported by grant CCR-4 from the National Science Foundation.
V. S. ADAMCHIK In this paper we aim at the most interesting case Gz = /Γ z, being the so-called double gamma function or the Barnes G-function: Gz + = Γz Gz, z C, G = G = G =, d log Gz, z >. d z The G-function has several equivalent forms including the Weierstrass canonical product Gz + = π z exp z + z + γ + z k z exp k k z and the functional relationship with the Hurwitz zeta function ζt, z, due to Gosper [] and Vardi [] log Gz + z log Γz = ζ ζ, z, Rz >, where γ denotes the Euler-Mascheroni constant and the derivative ζ t, z = d ζt, z is defined as an analytic continuation provided by dt the Hermite integral 8: ζ, z = z z log z 4 z log z + arctan x z z e πx dx + x logx + z e πx dx. Rz >, In particular, 4 ζ = ζ x log x, = e πx dx. For z, the G-function has the Stirling asymptotic expansion log Gz + = z z z logπ log A 4 + log z + O, z where A is the Glaisher-Kinkelin constant given by 5 log A = ζ.
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION Combining this definition with the representation 4, we obtain a new integral for log A: 6 log A = + logπ x log x π e x dx. The constant A originally appeared in papers by Kinkelin [] and Glaisher [7, 8, 9] on the asymptotic expansion when n of the following product p p...n np, p N. Similar to the Euler gamma function, the G-function satisfies the multiplication formula []: 7 Gn z = e ζ n n n z / nz+5/ n π n nz i= n j= G z + i + j n, n N 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 [], and in entries 6.44, 8. by Gradshteyn and Ryzhik []. Recently, Vardi [], Quine and Choi [7], and Kumagai [4] computed the functional determinants of Laplacians of the n-sphere in terms of the Γ n -function. Choi, Srivastava [, 4] and Adamchik [, 6] presented relationships including integrals and series between other mathematical functions and the multiple gamma function. The Γ n function has several interesting applications in pure and applied mathematics and theoretical physics. The most intriguing one is the application of G-function to the Riemann Hypothesis. Montgomery [6] and Sarnak [9] see [] for additional references have conjectured that the non-trivial zeros of the Riemann zeta function are pairwise distributed like eigenvalues of matrices in the Gaussian Unitary 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 matrices, 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 M n = B B... B n B B... B n+... B n B n+... B n where B n are the Bell numbers, defined by B n = n k= { } n k and { n k} are the Stirling subset numbers [8] 8 { } n = k k { } n + k It was proved in [5], that { } n, k det M n = Gn +, n N. { } { n, n =, =, n. 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]: 9 log G = 8 + log 4 log π log A, 4 where A is the Glaisher-Kinkelin constant 5. The next two special cases are due to Srivastava and Choi [4]: log G = 4 G 4π 4 log Γ 9 log A, 4 8 log G = log G + 4 4 G π log log π 8 4 where G denotes Catalan s constant: + log Γ, 4 G = k= k k +. Identities and are proven by using the duplication formula 7 setting n = together with 9. More special cases for z = 6,,
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 5 and 5 can be derived as well. Here is one of such formulas see [] for 6 the others: Proposition.. log G = 9 + log 7 + π log Γ 8 4 log A ψ π, where ψ z = z log Γz is the polygamma function. Proof. We recall the Lerch functional equation for the Hurwitz zeta function see [, 8] ζ n, z + n ζ n, z = e πin/ n! π n Li n+ πi n + B n+z+ e πiz, < z <, n N, where ζ t, z = d dt ζt, z, B nz are the Bernoulli polynomials and Li n z is the polylogarithm. Setting n = and using yields the reflexion formula for the G-function see also [4, 6]: log G + z = z log G z sinπz where Cl z is the Clausen function defined by π π Cl πz, < z <, 4 Cl z = ILi e iz. Upon setting z = we obtain to, and making use of Cl π = ψ π 9, 5 G G = π 6 Γ exp π ψ 8π.
6 V. S. ADAMCHIK On the other hand, using the multiplication formula 7 with z = and n =, we find 6 G G = 7/ e / πa 8 Γ. Now, combining 5 and 6, 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 functions and constants for other rational arguments. The most general result available for the multiple gamma function Γ n z is a closed form at z = : Proposition.. For n N 7 n n! log Γ n!! log π = n + n n n log P k,n B k+ k + k + k k P k,n ζ k, where P k,n are coefficients by xk, k n in expansion of n x + j j= Proof. The proof follows directly from [4], formula 7, by setting z =. Here are a few particular cases: the gamma function n = : the Barnes function n = : Γ = π G = /4 e /8 A / π /4 the triple gamma function n = : Γ = A/ π /6 /4 exp 7 ζ π 8
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 7. Hermite integrals In this section we derive a few integral representations for the Barnes function. We start with recalling the Hermite integral see [7, 5] for the Hurwitz zeta function: 8 ζs, z = z s + z s s + sins arctanx/z x + z s/ e πx dx, which provides analytic continuation of ζs, z to the domain s C {}. Differentiating both sides of 8 with respect to s, letting s = and s =, and using the second Binet formula for log Γz see [7] 9 log Γz = z logπ arctanx/z log z z + + dx e πx together with, we readily obtain: log Gz + = z log z H z ζ ζ x logx + z e πx dx, Rz >, where H p = p /k are the harmonic numbers and ζ z is a derivative of the Riemann zeta function. In the same manner, repeating the above steps n times, we derive the general integral representations. Proposition.. For Rz > n! log Γ n+ z + = z n+ log z H n+ n + n k k! n k= k n k { } n log Γ k+ z + + k n x n arctan x/z dx e πx ζ k z n k
8 V. S. ADAMCHIK and n! log Γ n z + = zn log z H n + n n n k ζ kz n k + k k= n { n k k! k } log Γ k+ z + n x n logx + z e πx where { n k} are the Stirling subset numbers, defined in 8. Upon substituting z = into formulas and, we obtain new integrals n 4 n x n arctanx e πx x n log + x e πx dx dx = H n n+ n n + + k k dx = H n n n + k= k= ζ k n k ζ k k Performing simple transformations of the integrands in and 4 will apparently lead to a large number of special integrals: x log x e πx dx = ζ x log x e πx + dx = ζ + log x log + x e πx + dx = 4 4 log + ζ arctan x z logπ dx = z log z z + log Γ z +, Rz > e πx +
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 9 arctanx 4 dx = + log e πx + xdx e πx + x + z = ψ z + π log z, argz xdx e πx x + z = log z ψz z, argz π 4. Binet-like representation In this section we derive the Binet integral representation for the Barnes function. Proposition 4.. The Barnes G-function admits the Binet integral representation: 5 log Gz + = z log Γz + z 4 log z B z log A e zx x e x x x dx, Rz >. Proof. Recall the well-known integral for the Hurwitz function [8] 6 ζs, z = Γs x s e zx dx, Rs >, Rz > e x The integral in the right hand-side of 6 can be analytically continued to the larger domain Rs > by subtracting the truncated Taylor series of e x at x =. Since 7 and e x = x + + x + Ox x s e zx Γs x + + x dx = s z s + z s z s s, Rs >, Rz > we therefore can rewrite 6 as
V. S. ADAMCHIK ζs, z = z s s + z s + s z s + x s e zx Γs e x x x dx, Rs >. Differentiating the above formula with respect to s and computing the limit at s =, we obtain ζ, z = z 4 + log z B z+ e zx x e x x x dx, Rz >, where B z is the second Bernoulli polynomial. The convergence of the above integral is ensured by 7. Finally, using, we conclude the proof. 5. Asymptotic expansion of Γ n z The asymptotic expansion for Γ n z when z follows straightforwardly from formulas and by expanding log + x and arctanx into the Taylor series, and then performing formal term-byterm integration. Lemma 5.. For Rz > 8 log Gz + = z log z log z z ζ + ζ n B k+ 4 k k + z + O k zn+, z. Proof. We start with the integral representation. By recalling the Taylor series for log + x we find that log + x z = n k= k k + xk+ n x y n+ z z n+ y + z dy,
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION n k= x log + x /z e πx k k + z k+ dx = x k+ n dx e πx z n+ t dt t e πt y n+ y + z dy which, by appealing to formula.6.4 in [7] 9 t k e π t dt = k+ B k 4 k, k N, reduces to n+ x log + x /z e πx dx = B k+ n 4 k k + zk t dt t/z e π t y n+ + y dy Substituting this into, we complete the proof. Lemma 5.. For Rz > log Γ z + = z log z 6 log z z ζ + ζ n + z log z + ζ + ζ + ζ B k+ n 4 k k + z B k+ k k + k z + k O zn+, z.
V. S. ADAMCHIK Proof. Employing the same technique as in getting 8, we begin with the integral representation of the triple gamma function: log Γ z + = z log z H + log Gz + + z ζ zζ + ζ + x arctanx/z e π x that follows directly from Proposition. with n =. Expanding arctan into the Taylor series arctan x z = n and making use of 9, we readily obtain x arctanx/z e π x k k xk n x y n + z z n y + z dy dx = n B k+ 4 k + k z k + n t dt y n z n e π t y + z dy. The conclusion of the lemma follows immediately by substituting this into formula and using Lemma 5.. t dx 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 log Γxdx = z z + z logπ + z log Γz log Gz +, where Rz >. Employing integration by parts, we derive z z log Gz+ = + z z logπ+ xψx dx, Rz > The integral representation provides an efficient numeric procedure for arbitrary precision evaluation of the G-function. This representation demonstrates that the complexity of computing Gz depends at
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION most on the computational complexity of the polygamma function. The later can be numerically computed by the Spouge approximation [] a modification of the Lanczos approximation. For the numerical integration in we shall use the Gaussian quadrature scheme. The restriction Rz > in can be easily removed by analyticity of the polygamma. For example, by resolving the singularity of the integrand at the pole x =, we continue log Gz+ to the wider area Rz > : log Gz+ = z z + z z logπ+logz++ xψx dx x + Another method of the analytic continuation of the G-function is to use the following well-known identity ψx = ψ x π cotπx x which, upon substituting it into, yields log G z = log G + z z logπ + z πx cotπx dx. The identity holds everywhere in a complex plane of z, except the real axes, where the integrand has simple poles. Therefore, combining with yields the following definition of the Barnes function 4 G n =, n N, Gz = π z / exp z z + xψx dx, γ argz π, where the contour of integration γ is a line between and z that does not cross the negative real axis; for example, γ could be the following path {, i, i + z, z }.
4 V. S. ADAMCHIK It remains to define Gz for z R. We do this by using the reflexion formula, which upon periodicity of the Clausen function 4 can be rewritten as follows: 5 G z = z/ Gz + sinπz z+ π exp π Cl πz z, z R. 7. Glaisher s constant In this section we discuss numeric computational schemes for Glaisher s constant: 6 A = exp ζ In light of the functional equation for the zeta function πs ζs = s π s Γ sζ s sin we can rewrite 6 as γ 7 A = exp ζ π /. π Exactly this formula is used in Mathematica V4. for computing Glaisher s constant. Another computational scheme for Glaisher s constant is based on the Barnes G-function identity: 8 log A = + log 6 log π 6 log G. Here G can be computed via integral. A more efficient algorithm is provided by the following identity 9 log A = log + 6 ζk + 8 + + k 6, + k which allows to approximate Glaisher s constant as
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 5 log A = log + 6 N ζk + 8+ + k 6 +O + k 4 N This approximation is suitable for arbitrary precision computation - it requires p log terms to achieve p decimal-digit accuracy. We close this section with a proof of formula 9. Recall that [8] 4 ζs = x s dx, Rs >. Γs e x e x Upon replacing ζk + in 9 by its integral representation 4 and inverting the order of summation and integration, we obtain 4 log A = log + dx 6 e x e x Next, noting that and t k k! k + s = t s we compute the inner sum t k k! = cosht t t k 8+ k! + k 6. + k x s coshx dx, t 4 t k 8 + k! + k 6 = 6 + t 4 t 4 + + k 6 + 5 t + 4 t 4 cosht 6 t + t sinht In the next step, we substitute this back to 4 and, employed by 6, evaluate the integral. Unfortunately, integration cannot be done term by term. The way around is to multiply the integrand by t λ, λ > and consider the limit as λ. This yields
6 V. S. ADAMCHIK 6 t λ 4 e t + 6 tλ 4 e t e t 6 t λ 5 tλ + e t e t e t + t λ 4 tλ + e t e t e t + 7 tλ 4 6 tλ e t e t e t + 5 tλ e t e t + 4 t λ e t e t 6 tλ e t e t tλ e t + 8 tλ e t e t dt = 6 Γλ 9 5 λ Γλ 6 Γλ 9 4 λ Γλ Γλ 5 λ Γλ Γλ Γλ λ + 4 Γλ + 7 Γλ + λ + 6 Γλ ζλ. Computing the limit as λ, we find that the right hand-side of the above expression evaluates to This completes the proof of 9. log 6 ζ.
CONTRIBUTIONS TO THE THEORY OF THE BARNES FUNCTION 7 References [] V. S. Adamchik, H. M. Srivastava, Some series of the zeta and related functions, Analysis 998 44. [] V. S. Adamchik, Polygamma functions of negative order, J. Comp. and Appl. Math. 998 9 99. [] V. S. Adamchik, On the Barnes function, in Proceedings of the International Symposium on Symbolic and Algebraic Computation, London, Ontario; July -5,, ACM Press, New York,, 5. [4] V. S. Adamchik, Multiple Gamma Function and Its Application to Computation of Series, preprint,, http://lanl.arxiv.org/abs/math.ca/874. [5] M. Aigner, A characterization of the Bell numbers, Discrete Math. 5999 7. [6] W. P. Alexejewsky, Ueber eine Classe von Funktionen, die der Gammafunktion analog sind, Leipzig Weidmanncshe Buchhandluns 46894 68 75. [7] G. Andrews, R. Askey, R. Roy, Special Functions, Cambridge University Press, Cambridge, 999. [8] H. Bateman, A. Erdelyi, Higher Transcendental Functions, Vol., McGraw- Hill, New York, Toronto, London, 95. [9] E. W. Barnes, The theory of the G-finction, Quart. J. Math. 899 64 4. [] E. W. Barnes, Genesis of the double gamma function, Proc. London Math. Soc. 9 58 8. [] E. W. Barnes, The theory of the double gamma function, Philos. Trans. Roy. Soc. London, Ser. A 969 65 88. [] E. W. Barnes, On the theory of the multiple gamma function, Trans. Cambridge Philos. Soc. 9 94 74 45. [] J. Choi, H. M. Srivastava, J. R. Quine, Some series involving the zeta function, Bull. Austral. Math. Soc. 5 995 8 9. [4] J. Choi, H. M. Srivastava, Certain classes of series involving the zeta function, J. Math. Anal. and Appl. 999 9 7. [5] J. Choi, H. M. Srivastava, Certain classes of series associated with the Zeta function and the multiple Gamma function, J. Comput. Appl. Math. 8 87 9. [6] J. Choi, H. M. Srivastava, V. S. Adamchik, Multiple Gamma and related functions, Appl. Math. and Comp. 4 55 5. [7] J. W. L. Glaisher, On a numerical continued product, Messenger of Math. 6 877 7 76. [8] J. W. L. Glaisher, On the product,,,...,n n, Messenger of Math. 7 887 4 47. [9] J. W. L. Glaisher, On certain numerical products, Messenger of Math. 89 45 75. [] R. W. Gosper, m 6 n log Γzdz, in: M. Ismail, D. Masson, M. Rahman editors, 4 Special functions, q-series and related topics. The Fields Institute Communications, Amer. Math. Soc., 4, 997 7 76. [] I. S. Gradshteyn, I. M. Ryzhik, Table of Integrals, Series, and Products, Sixth Ed. Academic Press, San Diego,.
8 V. S. ADAMCHIK [] J. P. Keating, N. C. Snaith, Random matrix theory and ζ, + it Commun. Math. Phys. 4 57 89. [] Kinkelin, Ueber eine mit der Gammafunction verwandte Transcendente und deren Anwendung auf die Integralrechnung, J. Reine Angew. Math. 57 86 58. [4] H. Kumagai, The determinant of the Laplacian on the n-sphere, Acta Arithmetica XCI. 999 99 8. [5] W. Magnus, F. Oberhettinger, R. P. Soni, Formulas and theorems for the special functions of mathematical physics Springer, Berlin, 966. [6] H. L. Montgomery, The pair correlation of the zeta function, Proc. Symp. Pure Math. 497 8 9. [7] J. R. Quine, J. Choi, Zeta regularized products and functional determinants on spheres, Rocky Mountain Journal 6 996 79 79. [8] K. N. Rosen and et.., Handbook of Discrete and Combinatorial Mathematics CRC Press, New York,. [9] P. Sarnak, Quantum chaos, symmetry and zeta functions, Curr. Dev. Math. 997 84 5. [] J. L. Spouge, Computation of the gamma, digamma, and trigamma functions, SIAM J. Numer. Anal. 9949 944. [] I. Vardi, Determinants of Laplacians and multipe gamma functions, SIAM J. Math. Anal. 9 988 49 57. [] M. F. Vignéras, L équation fonctionelle de la fonction zêta de Selberg du groupe mudulaire SL, Z, Astérisque 6 9795 49. [] 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, 96. Carnegie Mellon University, Pittsburgh, USA E-mail address: adamchik@cs.cmu.edu