On the coefficients of the Báez-Duarte criterion for the Riemann hypothesis and their extensions arxiv:math-ph/ v2 7 Sep 2006

Transcription

1 On the coefficients of the Báez-Duarte criterion for the Riemann hypothesis and their extensions arxiv:math-ph/ v2 7 Sep 2006 Mark W. Coffey Department of Physics Colorado School of Mines Golden, CO 8040 {mcoffey@mines.edu} (Received 2006 September 6, 2006 Abstract We present analytic properties and extensions of the constants c k appearing in the Báez-Duarte criterion for the Riemann hypothesis. These constants are the coefficients of Pochhammer polynomials in a series representation of the reciprocal of the Riemann zeta function. We present generalizations of this representation to the Hurwitz zeta and many other special functions. We relate the corresponding coefficients to other known constants including the Stieltes constants and present summatory relations. In addition, we generalize the Maślanka hypergeometric-like representation for the zeta function in several ways. Key words and phrases Báez-Duarte criterion, Riemann hypothesis, Stieltes constants, Pochhammer polynomial, Riemann and Hurwitz zeta functions, polygamma function, complete Bell polynomials, Maślanka representation, Dirichlet L function

2 Introduction As reformulated by Báez-Duarte [3], the Riemann hypothesis (RH is equivalent to a certain growth condition on constants c k that appear in a series representation of the reciprocal of the Riemann zeta function ζ(s. In this paper we analytically investigate these constants and introduce some parametrized extensions. We point out the importance of such extensions for future work. We relate parametrized coefficients c k (b,a to other important constants of analytic number theory including the Stieltes constants and present summatory relations for the former. With c k =0( ( k, k 0, ( ζ(2 +2 unconditionally c k = O(k /2 and on the RH we have c k = O(k 3/4+ǫ for any ǫ > 0 [3]. The constants c k are known to consist of a relatively rapidly decreasing term /k 2 and an oscillatory contribution, as can be shown by Rice s integrals or other asymptotic methods [4, 20]. Not surprisingly, it is the detailed behaviour of the oscillatory contribution upon which the validity of the RH depends. Numerical results for c k are presented in Refs. [5, 6, 4, 20] and the first billion values have been reported. These values are consistent with the Báez-Duarte criterion under the RH. If further c k = O(k 3/4, then the complex zeros of ζ(s are on the critical line Re s = /2 and are simple. Summatory relations for c k In this section we relate c k to the Stieltes constants γ k. The latter are the coeffi- 2

3 cients of the Laurent series of the Riemann zeta function about s =. In preparation we have Lemma concerning the derivatives of Pochhammer polynomials. Let P k (s ( s k / = ( k Γ(s/Γ(s k, where (a n = Γ(a + n/γ(a is the Pochhamer symbol. Let ψ = Γ /Γ be the digamma function, where Γ is the Gamma function, and ψ ( the polygamma function [2]. Lemma. Set g(s [ψ( s/2 ψ(k + s/2], 2 (2a with Then we have and g (l (s = ( l 2 l+ [ ψ (l ( s/2 ψ (l (k + s/2 ], (2b ( ( d s s ds P k = P k g(s, (3 2 2 ( ( ( d s s [ P k = P k Y g(s,g (s,...,g ( (s ], (4 ds 2 2 where Y are (exponential complete Bell polynomials []. Proof. Equation 3 follows from the definition of the Pochhammer symbol and Eq. (4 by Lemma of Ref. [7]. A very small subset of the relations between the constants c k and the Stieltes constants is contained in the following. Equation(5c presents how the Euler constant γ = ψ( may be written in terms of c. Proposition. We have Γ(k +/2 c k = 0, 3 (5a

4 and 2 π Γ(k +/2 ψ(k +/2c k =, (5b Γ(k +/2 [ψ 2 (k +/2+ψ (k +/2]c k = 4 π(γ 2ln2. (5c As prelude to the proof of Proposition, we know that the representation ζ(s = c k P k (s/2, (6 holds unconditionally in the half plane Re s >, converging uniformly on compact sets [3]. Proof of Proposition. We combine the definition of the Stieltes constants (e.g., [2] with Eq. (6 to write [ ] ζ(s = s + γ k (s k = c k P k (s/2 (7a = s γ(s 2 +(γ 2 +γ (s 3 +( γ 3 2γγ γ 2 /2(s 4 +O[(s 5 ]. (7b Here we used γ 0 = γ. We then Taylor expand the right side of Eq. (7a about s = using Lemma. In particular, we have ( ( d s Γ(k +/2 [ P k = Y g(s,g (s,...,g ( (s ] ds 2, (8 π s= s= and g( = 2 [ψ(/2 ψ(k +/2] = Γ(k +/2 2 π k l=0 l+/2, (9a g (l ( = ( l [ ψ (l (/2 ψ (l (k +/2 ] = l! k 2 l+ 2 l+ =0 ( +/2 l+. (9b 4

5 Equations (9a and (9b follow from the functional equations of the digamma and polygamma functions respectively. We next equate successive like powers of s on both sides of Eq. (7. Effectively, we evaluate the successive derivatives of the representation (6 as s +. We first obtain Eq. (5a. We then use it to obtain Eq. (5b. We then use both of these to obtain Eq. (5c and Proposition follows. Remarks. The continuation of the process ust described yields the explicit relation between γ k and sums over the constants c. The appearance of the low order Bell polynomials Y 0 =, Y (x = x, and Y 2 (x,x 2 = x 2 +x 2 is implicit in writing Eqs. (5a-(5c. They do not appear in the final results there since we performed successive manipulations. Equation (5a is a reflection of the simple pole of ζ(s at s = and of the relation n= µ(n/n = 0, where µ is the Möbius function. Once results such as Eqs. (5 have been derived, they may be directly verified using the alternative expression [3] For instance, we recover c k = n= Γ(k +/2 µ(n n 2 c k = π ( n 2 k. (0 n= µ(n n = 0. ( n2 Corollary. In Eq. (5a, the ratio Γ(k +/2/ k /2 for k while in Eq. (5b the factor ψ(k + /2 lnk as k. Therefore, the former equation shows that if the c k s did not change sign, they would have to decrease at least as fast as k /2 as k. 5

6 The following expression may be combined with the Stieltes constant expansion of /ζ(s given in Eq. (7 to write γ in terms of sums of the constants c k and Stirling numbers of the first kind s(k,l [, 6, 7]. Lemma 2. The coefficient of (s on the right side of Eq. (7a is given by ( kc k l= s(k, l 2 l ( l ( l. (2 Proof. We first re-express the Pochhammer polynomials using the Stirling numbers s(k, l, then binomially expand and reorder sums: c k ( s/2 k = = ( kc k l=0 = ( kc k =0 l= c k s(k, l 2 l ( k+l s(k,l( s/2 l l=0 l =0( l ( l (s ( l s(k, l 2 l ( l (s. (3 First extension of the constants c k The Hurwitz zeta function ζ(s,a = (k+a s, a / N 0 for Re s > extends to an analytic function with only a simple pole at s =. Here we demonstrate the representation Corollary 2. For Re s > we have ζ(s,a = c k (ap k (s/2, (4 with c k (a = =0( ( k, k 0. (5 ζ(2 +2,a 6

7 In particular, there results at a = /2 Corollary 3. ζ(s = (2s c k (/2P k (s/2. (6 Corollary 2 follows from Proposition 2. Define for b > +δ with δ > 0 and Re a > 0 the functions F(x,b,a ( k+ x k ζ(bk, a (k! (7 and ϕ(s,b,a 0 x (s/b+ F(x,b,adx, < Re s < b. (8 Put c k (b,a =0( ( k ζ(b +b,a. (9 Then we have the series and integral representations for Re s > ζ(s,a = In the proof we use the infinite series c k (b,ap k (s/b = ϕ(s,b,a Γ( s/b. (20 Lemma 3 and k= ( k x k = ( k + (k +! x k+ = x! ex, (2 Lemma 4 [ ] x F n=0 (n+a b,b,a = xe x. (22 The series (7 is uniformly convergent on compact sets of the complex x plane so that the interchange of sums used to show Eq. (22 is valid. 7

8 We have [ ] x (s/b+ x F 0 (n+a b,b,a dx = (n+a s ϕ(s,b,a. (23 Without assuming the RH, ϕ defined in Eq. (8 converges absolutely and uniformly in the strip specified. The behaviour F(x,b,a x/ζ(bk,a as x 0 dictates the condition Re s < b. The behaviour F(x,b,a = o(x /b as x without the RH gives the requirement Re s >. Then summing both sides of Eq. (23 on n from 0 to and using Lemma 4 gives ζ(s,aϕ(s,b,a = so that we have obtained the outer equality of Eq. (20. 0 x s/b e x dx = Γ( s/b, (24 We next re-express the function ϕ in terms of Pochhammer polynomials. We have from Eqs. (7 and (8 ϕ(s,b,a = = 0 =0 We next apply Lemma 3 so that ϕ(s,b,a = = 0 0 =0 0 =0 = ( ζ(b +b,a! x s/b dx (25a =0 ( e x ζ(b +b,a! e x x s/b dx. (25b ( ζ(b +b,a ( ζ(b +b,a k= ( k e x x k s/b dx ( k e x x k s/b dx (26a (26b c k(b,aγ(k + s/b. (26c 8

9 In obtaining Eq. (26b from (26a we reordered the double series and (26c used the definition (9. We have therefore found that ϕ(s,b,a = Γ( s/b c k (b,ap k (s/b, (27 and Proposition 2 is completed. We have shown a way to directly relate the RH criteria of Báez-Duarte [3] and of Riesz [5]. The summatory function appearing in the Riesz criterion corresponds to R(x F(x,2, in Eq. (7 and under the RH it is O(x /4+ǫ for ǫ > 0. A recent construction for the function /ζ similar to Proposition 2 has been given in Ref. [5]. The authors of that reference used the Möbius function in that development, whereas we have proceeded differently and obtained a result also applying to the reciprocal of the Hurwitz zeta function. Arguing as we have in Proposition 2 gives many extensions. An example is Corollary 4. Putting G(x ϕ G (x = ( k 0 x k ζ(2k+, (28 x (s+/2 G(xdx, (29 and we have for Re s > the representations ζ(s = c G k =( ( k ζ(2 +, (30 c G kp k [(s+/2] = ϕ G (s Γ[( s/2]. (3 9

10 Given Corollary 2 and Proposition 2 Eqs. (7 extend to [ ζ(s,a = s + ] γ k (a (s k = c k (b,ap k (s/b = s +ψ(a(s 2 +[ψ 2 (a+γ (a](s 3 +[ ψ 3 (a+2ψ(aγ (a γ 2 (a/2](s 4 +O[(s 5 ], wherein γ 0 (a = ψ(a has been used. The Stieltes constants γ k (a may be written in terms of sums containing the Bernoulli numbers B and elementary constants such as powers of ln2. Further properties of γ k (a are given in the very recent Refs. [8] and [9]. (32 Remark. Analogous to Eqs. (4 and (5 it is not difficult to show that the Maślanka representation for ζ(s [4] may be extended to where ζ(s,a = s A k (a ( s/2 k, (33 A k (a =0( ( k (2 +ζ(2 +2,a. (34 As a Corollary, we obtain as special cases a representation for Bernoulli polynomials B n (x = nζ( n,x and polygamma functions ψ (m (z = ( m+ m!ζ(m+,z. Additional cases include representations for the alternating Hurwitz zeta function, the digamma function, and the function β(x = (/2[ψ[(x + /2] ψ(x/2] = ( k /(x+k. Remark 2. We believe that it is very useful to have the constants c k extended to include one or more parameters. In that case manipulations on c k (a and /ζ(s,a 0

11 for instance may be performed with respect to a and/or s and then the Riemann zeta function case recovered by putting a = or /2. For instance, we have from Eq. (5 d da c k(a = =0( ( k (2 +2ζ(2 +3,a. (35 ζ 2 (2 +2,a Remark 3. The representation (33 may be further extended to the Hurwitz-Lerch zeta function Φ(z,s,a = n=0 z n /(n+a s, where s C for z <, Re s > when z =. In this case we have Proposition 3 where Φ(z,s,a = s A k (z,a ( s/2 k, (36 A k (z,a =0( ( k (2 +Φ(z,2 +2,a. (37 The proof of Proposition 3 and an extension of it is contained as a special case of Proposition 4 proved below. As a Corollary, we obtain Maślanka-type representations for polylogarithm functions (Jonquière s function Li s (z = zφ(z,s,, (38 where s C for z <, Re s > when z =. Remark 4. We have Conecture. For a class M of analytic functions expressible as a Dirichlet series and possessing at most polar singularities in the complex plane such that for f M and q > we have the representation f(s = s A k P k (s/q, (39

12 holding in a half plane of C with A k =0( ( k (q +q f(q +q. (40 We believe that a demonstration of some form of this Conecture is possible by an appropriate application of sampling theory (cf. the Appendix. However, a more expedient approach may be to apply known results for approximating analytic functions in terms of the zeta function. Then given the representation (33 Conecture would follow. In support of Conecture we have the following Proposition 4. Let f(s,a = n=0 f n (n+a s, (4 where it is assumed that {f n } 0 is such that the series converges in a half-plane Re s > σ > and a C/{0,,...}. Then we have for p > with f(s,a = s A k (ap k (s/p, (42 A k (a =0( ( k (p +p f(p +p,a. (43 Proof. We proceed as in Ref. [4] in the case f n =, a =, and p = 2, forming α s (s f(s,a = α α s n=0 f n (n+a s. (44 After term-by-term differentiation of the series we evaluate at α =. The interchange of various operations is ustified by the assumption on the sequence f n and the estimates of Ref. [4]. 2

13 As a Corollary, we obtain generalized Maślanka-type representations for other important special functions. These include the multiple zeta function ( ζ n (s,z = (k +k k n +z = k +n, (45 s (k +z s n k =0 k n=0 and the multiple Gamma function Γ n. The latter function has a product representation and may be defined by the recurrence-functional equation Γ n+ (z + = Γ n+ (z/γ n (z, Γ (z = Γ(z, Γ n ( =,forz C andn N +. ThemultipleGamma function may be expressed in terms of derivatives of the multiple zeta function [9]: ( ζ n (s,z n z lnγ n (z = lim + ( k R n+ k, (46a s 0 s k where n R n = lim s 0 Second extension of the constants c k ζ n (s,. (46b s Another possible extension of the constants c k would be to write c k (a,b =0( ( k, k 0, (47 ζ(a +b wherein we considered a = b in Eq. (9. This extension has also been observed in Ref. [5] and numerical experiments presented. However, we remain with the case a = b = 2 and instead note that ζ(2 +2 = (2π2+2 ( B 2+2 2(2 +2! Then we may consider = (2π2+2 ( 2(2 +2! B 2+2 (x x=0. (48 c k (x = =0 ( k 2(2 +2! (2π 2+2 B 2+2 (x, (49 3

14 such that B n(x = nb n (x. If f(x = /B 2+2 (x, then we have d dx f(x = f(x(2 +2B 2+(x B 2+2 (x, (50 and Lemma of Ref. [7] applies for the higher order derivatives of f(x in terms of (exponential complete Bell polynomials. When evaluated at x = 0, B 2+ (0 = 0 unless = 0 when B = /2, and Ref. [8] describes the Bell polynomials when the odd-indexed variables are set to zero. Connection with the η constants We first mention in passing the following that recovers a result of [3], but in a different way. We have Lemma 3 lim P k(s(k + s = k Γ( s. (5 Proof. We write P k (s = ( s k / = Γ(k + s/γ( s, apply the known asymptotic form of Γ(z +a/γ(z +b (e.g., [], and take the limit. We now introduce the constants η of the Laurent expansion ζ (s ζ(s = s p=0 η p (s p, s < 3, (52 with η 0 = γ. These coefficients are important in the theory of the function lnζ(s; hence they are connected with the behaviour of the prime counting function π(x. The alternating binomial sum S 2 (n k = ( ( k η is key in the Li criterion for the RH; its sublinearity in n would imply the latter conecture (e.g., [0]. 4

15 We apply the identity d dsζ(s = ζ (s ζ(s ζ(s, (53 and evaluate the derivatives as s from the right. We obtain Proposition 5. Put a 0 = Γ(k +/2/ π and a q = Γ(k +/2 [ Y q g(s,g (s,...,g (q (s ], q, (54 q! π s= where the function g and its derivatives are given in Eq. (2. We then have for q q+ q qa q+ c k = a c k η q. (55 = Proof. We first recall from the representation (6 that ζ(s = c k q=0 ( d q! ds q ( s P k (s q 2 s= (56a = c k a q (s q. q= (56b For Eq. (56a we have kept in mind that P k (s is a polynomial of degree k in s. For Eq. (56b we have used Lemma together with Eq. (5a. We then reorder the double summation, obtaining and ζ(s = q a q c k (s q, q=0 q= d dsζ(s = q+ (q +a q+ c k (s q. (57a (57b We then apply identity (53, a form of which is ( d s ( s c k ds P k = k P k 2 c 2 s + η p (s p. (58 p=0 5

16 We carry out the necessary multiplication of series on the right side of Eq. (58 and reorder the second term. We then equate coefficients of like powers of s on both sides and Eq. (55 follows. Other summatory relations The authors of the extremely recent Ref. [6] derived the identity c k s k = s ( s k s, ζ(2k+2 Re s < /2 (59 and noted the value ( k c k = 2 k /ζ(2k We first illustrate that Eq. (59 can provide the basis of a family of summatory relations and have Proposition 6. For Re t /2 we have giving Corollary 5. c 0 ln( t+ c k k tk = ( k ( t k, (60 kζ(2k+2 t and c k k c 0 ln(2/3+ c k k 3 k = ( k , (6a kζ(2k+2 2k c 0 ln2+ ( kc k k = , kζ(2k+2 2k (6b c k c 0 ln2+ k 2 = ( k k kζ(2k , (6c [( k + ] = 2 k [( k + ] kζ(2k+2 2 k (6d 6

17 Proof of Proposition 6. We write Eq. (59 in the form c 0 s s + c k s k = s ( s k s ζ(2k+2. (62 We then divide both sides by s, integrate on s from 0 to t, and Eq. (60 obtains. Special cases of t in Eq. (60 yield Corollary 5. Equations (59 and (60 may represent the only so far known series associated with the zeta function where reciprocal zeta values at integer argument occur in the summand. Equation (59 is extended by Proposition 7. For b > and Re a > 0 we have c k (b,as k = s ( s k s, ζ(bk +b,a Re s < /2, (63 where c k (b,a is defined in Eq. (9. In particular, we have Corollary 6. ( k c k (b,a = 2 k ζ(bk,a. (64 Proof. We use the definition (9, reorder a double sum, and apply the binomial expansion: c k (b,as k = = = =0 = =0 =0 s k =0( ( k ( ζ(b +b,a ( s ζ(b +b,a k=s k( k s k ( k + ( s ζ(b +b,a ( s ζ(b +b,a +. (65 7

18 The alternating sum (64 obtains at s =. Similarly, Proposition 6 and Corollary 5 may be extended to include the values c k (b,a. We have Proposition 8. For Re t /2, b >, and Re a > 0 there holds We omit the proof. c 0 (b,aln( t+ c k (b,a t k = k ( k ( t k. (66 kζ(bk +b,a t Third extension of the constants c k Let χ be a Dirichlet character mod k and L(s, χ the corresponding Dirichlet L-function (e.g., [3] L(s,χ = n= χ(n, Re s >. (67 ns We recall that We have L(s,χ = n= χ(nµ(n n s, Re s >. (68 Proposition 9. For b > there holds where L(s,χ = c k (b,χp k (s/b, Re s >, (69 c k (b,χ =0( ( k L(b +b,χ. (70 Proof. We proceed as in Ref. [3]. By simply noting that χ(n the estimates given there ustify the interchange of infinite summations. 8

19 Final remarks The representation of Proposition 9 is expected to extend to automorphic L functions. Accordingly we expect a criterion on the rate of growth of c k (b,χ and its generalization to be equivalent to the extended and generalized Riemann hypothesis, respectively. Based upon a special case of Proposition 3 (or 4 there is an extended Maślanka type representation of Dirichlet L functions. This follows since Dirichlet L functions may be written as a combination of Hurwitz zeta functions. The analog of Stieltes constants and the constants η exist for Dirichlet L functions (e.g., Appendix E of Ref. [0] and our method of Proposition or 5 would equally well apply for relating them to c k (b,χ. If we insert the Euler product for ζ(s into the expression ( for c k we have c k = =0( ( k [ ] p (2+2 p [ ] = δ 0k =0( ( k p ( (7 p In Eq. (7 the product or sum over p is over all primes and δ k is the Kronecker symbol. The first sum in brackets on the right side of this equation may be estimated as p xp (2+2 Ei[ (2 +lnx], where Ei is the exponential integral. An approximate expression for c k for large values of k is given by [6] c k n= µ(n n 2 e k/n2. (72 9

20 We note that alternatively this approximation may be written as a Fourier transform: c k 2 π n= µ(n n 2 0 Summary coskt dt n 2 t 2 +n 2. (73 Our results generalize both the Báez-Duarte representation of the reciprocal of the Riemann zeta function [3] and the Maślanka representation of the zeta function itself [4]. The Maślanka representation has been generalized to the Hurwitz zeta function, the Hurwitz-Lerch zeta function, the multiple zeta function, and other important special functions. We have further extended the Báez-Duarte representation of /ζ to the representation of the reciprocal of Dirichlet L functions and it is anticipated that this may be generalized to automorphic L functions. By way of our generalization of the Báez-Duarte representation in terms of Pochhammer polynomials, we have effectively demonstrated the equivalence of the Riesz [5] and Báez-Duarte criteria for the Riemann hypothesis. We have obtained summatory relations for the coefficients c k of the Báez-Duarte criterion for the Riemann hypothesis and related them to important constants of analytic number theory. In describing the relation of c k to the Stieltes and other constants we have made use of the (exponential complete Bell polynomials Y. Acknowledgement I thank Prof. L. Báez-Duarte for his comments upon reading the manuscript, in particular in clarifying the statement of Conecture. 20

21 Appendix: Interpolating binomial sums from the Faà di Bruno formula A significant source of alternating binomial sums is the Faà di Bruno formula, a generalization of the chain rule. Put D z d/dz and x = x(z. Then we have D n z f(x = n D k x f(x( k =0( ( k x k D n z x. (A. In particular, we have for real a ( a+n n Dz n x a = a n =0( ( n a+ x a Dz n x, (A.2 or equivalently ( ( n a n+a x a Dz n x a = x Dz n =0 n x. (A.3 In connection with developing alternative representations of analytic functions, we point out that Eq. (A.3 can be viewed as an immediate consequence of Lagrange interpolation. This follows from ( a ( n+a n = n k k +a, 0 n. (A.4 k With the Faà di Bruno formula the exponential Bell polynomials again make an appearance []. 2

22 References [] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Washington, National Bureau of Standards (964. [2] G. E. Andrews, R. Askey, and R. Roy, Special functions, Cambridge University Press (999. [3] L. Báez-Duarte, A sequential Riesz-like criterion for the Riemann hypothesis, Int. J. Math. Math. Sci. 2, (2005. [4] L. Báez-Duarte, On Maślanka s representation for the Riemann zeta function, math.nt/ (2003. [5] S. Beltraminelli and D. Merlini, The criteria of Riesz, Hardy-Littlewood et al. for the Riemann hypothesis revisited using similar functions, math.nt/06038 (2006. [6] J. Cislo and M. Wolf, Equivalence of Riesz and Báez-Duarte criterion for the Riemann hypothesis, math.nt/ (2006. [7] M. W. Coffey, A set of identities for a class of alternating binomial sums appearing in computing applications, to appear in Util. Math. (2007. [8] M. W. Coffey, New results on the Stieltes constants: Asymptotic and exact evaluation, J. Math. Anal. Appl. 37, (2006; arxiv:math-ph/

23 [9] M. W. Coffey, New summation relations for the Stieltes constants, Proc. Royal Soc. A 462, (2006. [0] M. W. Coffey, Towards verification of the Riemann hypothesis, Math. Phys., Analysis and Geometry 8, (2005. [] L. Comtet, Advanced Combinatorics, D. Reidel (974. [2] A. Ivić, The Riemann Zeta-Function, Wiley (985. [3] H. Iwaniec and E. Kowalski, Analytic number theory, American Mathematical Society (2004. [4] K. Maślanka, Báez-Duarte criterion for the Riemann hypothesis and Rice s integrals, math.nt/ v2 (2006. [5] M. Riesz, Sur l hypothèse de Riemann, Acta Math. 40, (96. [6] J. Riordan, An introduction to combinatorial analysis, Wiley (958. [7] J. Riordan, Combinatorial identities, Wiley (968. [8] R. Schimming and W. Strampp, Differential polynomial expressions related to the Kadomtsev-Petviashvili and Korteweg-de Vries hierarchies, J. Math. Phys. 40, (999. [9] I. Vardi, Determinants of Laplacians and multiple gamma functions, SIAM J. Math. Analysis 9, (

24 [20] M. Wolf, Evidence in favor of the Báez-Duarte criterion for the Riemann hypothesis, math.nt/ (