arxiv:002.47v2 [math.nt] 28 Feb 200 Two arguments that the nontrivial zeros of the Riemann zeta function are irrational Marek Wolf e-mail:mwolf@ift.uni.wroc.pl Abstract We have used the first 2600 nontrivial zeros γ l of the Riemann zeta function calculated with 000 digits accuracy and developed them into the continued fractions. We calculated the geometrical means of the denominators of these continued fractions and for all cases we get values close to the Khinchin s constant, what suggests that γ l are irrational. Next we have calculated the n-th square roots of the denominators q n of the convergents of the continued fractions obtaining values close to the Khinchin Lèvy constant, again supporting the common believe that γ l are irrational. Introduction Bernhard G.F. Riemann has shown [6] that the l.h.s. of the identity valid only for R[s] > : ζ(s) := n = ( ), s = σ + it () s p s n= p=2 can be analytically continued to the whole complex plane devoid of s = by means of the following contour integral: ζ(s) = Γ( s) ( x) s dx (2) 2πi e x x P
where the integration is performed along the path P Now dozens of integrals and series representing the ζ(s) function are known, for collection of such formulas see for example the entry Riemann Zeta Function in [20] and references cited therein. The ζ(s) function has trivial zeros 2, 4, 6,... and infinity of nontrivial complex zeros ρ l = β l + iγ l in the critical strip: β l (0, ). The Riemann Hypothesis (RH) asserts that β l = for all l i.e. all zero lie on the critical line R(s) =. 2 2 Presently it is added that these nontrivial zeros are simple: ζ (ρ l ) 0 many explicit formulas of number theory contain ζ (ρ l ) in the denominators. In 94 G. Hardy [6] proved that infinitely many zeros of ζ(s) lie on the critical line. A. Selberg [8] in 942 has shown that at least a (small) positive proportion of the zeros of ζ(s) lie on the critical line. The first quantitative result was obtained by N. Levinson in 974 [9] who showed that at least one-third of the zeros lie on the critical line. In 989 B. Conrey [2] improved this to two-fifths and quite recently with collaborators [] to over 4%. It was checked computationally [5] that the 0 3 first zeros of the Riemann Zeta function fulfill the condition β l =. A. Odlyzko checked that RH is 2 true in different intervals around 0 20 [0], 0 2 [], 0 22 [4], see also [5] for the two billion zeros from the zero 0 24. There is no hope to obtain the analytical formulas for the imaginary parts γ l of the nontrivial zeros of ζ(s) but the common belief is that they are irrational and perhaps even transcendental [3]. The problem of any linear relations between γ l with integral coefficients appeared for the first time in the paper of A.E. Ingham [7] in connection with the Mertens conjecture. This conjecture specifies the growth of the function M(x) defined by M(x) = n x µ(n), (3) where µ(n) is the Möbius function for n = µ(n) = 0 when p 2 n (4) ( ) r when n = p p 2... p r The Mertens conjecture claims that M(x) < x 2. (5) 2
From this inequality the RH would follow. A. E. Ingham in [7] showed that the validity of the Merten s conjecture requires that the imaginary parts of the nontrivial zeros should fulfill the relations of the form: N c l γ l = 0, (6) l= where c l are integers not all equal to zero. This result raised the doubts in the inequality (5) and indeed in 985 A. Odlyzko and H. te Riele [5] disproved the Merten s conjecture. In this paper we are going to exploit two facts about the continued fractions: the existence of the Khinchin constant and Khinchin Lèvy constant, see e.g. [4,.8], to support the irrationality of γ l. Let r = [a 0 (r); a (r), a 2 (r), a 3 (r),...] = a 0 (r) + a (r) + a 2 (r) + a 3 (r) +... be the continued fraction expansion of the real number r, where a 0 (r) is an integer and all a k (r) with k are positive integers. Khinchin has proved [8], see also [7], that ( lim a (r)... a n (r) ) n = n m= { + (7) } log2 m K 0 2.68545200... (8) m(m + 2) is a constant for almost all real r [4,.8]. The exceptions are of the Lebesgue measure zero and include rational numbers, quadratic irrationals and some irrational numbers too, like for example the Euler constant e = 2.78288285... for which the limit (8) is infinity. The constant K 0 is called the Khinchin constant. If the quantities K(r; n) = ( a (r)a 2 (r)... a n (r) ) n (9) for a given number r are close to K 0 we can regard it as an indication that r is irrational. Let the rational p n /q n be the n-th partial convergent of the continued fraction: p n q n = [a 0 ; a, a 2, a 3,..., a n ]. (0) For almost all real numbers r the denominators of the finite continued fraction approximations fulfill: ( qn (r) ) /n = e π 2 /2 ln 2 L 0 = 3.27582298728... () lim n where L 0 is called the Khinchin Lèvy s constant [4,.8]. Again the set of exceptions to the above limit is of the Lebesgue measure zero and it includes rational numbers, quadratic irrational etc. 3
2 The computer experiments First 00 zeros γ l of ζ(s) accurate to over 000 decimal places we have taken from [2]. Next 2500 zeros of ζ(s) with precision of 000 digits were calculated using the built in Mathematica v.7 procedure ZetaZero[m]. We have checked using PARI/GP [9] that these zeros were accurate within at least 996 places in the sense that in the worst case ζ(ρ l ) < 0 996, l =, 2,..., 2600. PARI has built in function contfrac(r, {nmax}) which creates the row vector a(r) whose components are the denominators a n (r) of the continued fraction expansion of r, i.e. a = [a 0 (r); a (r),..., a n (r)] means that r a 0 (r) + a (r) + a 2 (r) +... a n (r) (2) The parameter nmax limits the number of terms a nmax (r); if it is omitted the expansion stops with a declared precision of representation of real numbers at the last significant partial quotient. By trials we have determined that the precision set to \p 2200 is sufficient in the sense that scripts with larger precision generated exactly the same results: the rows a(γ l ) obtained with accuracy 2200 digits were the same for all l as those obtained for accuracy 2600 and the continued fractions accuracy set to 200 digits gave different denominators a n (γ l ) With the precision set to 2200 digits we have developed the 000 digits values of each γ l, l =, 2,... 2600, into the continued fractions γ l. = [a0 (l); a (l), a 2 (l), a 3 (l),..., a n(l) (l)] a(l) (3) without specifying the parameter nmax, thus the length of the vector a(l) depended on γ l and it turns out that the number of denominators was contained between 788 and 2072. The value of the product a a 2... a n(l) was typically of the order 0 800 0 870. Next for each l we have calculated the geometrical means: n(l) K l (n(l)) = a k (l) k= /n(l). (4) The results are presented in the Fig.. Values of K l (n(l)) are scattered around the red line representing K 0. To gain some insight into the rate of convergence of K l (n(l)) we have plotted in the Fig. 2 the number of sign changes S(l) of K l (m) K 0 for each l when m = 00, 0,... n(l), i.e. S K (l) = number of such m that (K l (m + ) K 0 )(K l (m) K 0 ) < 0. (5) 4
The largest S K (l) was 22 and it occurred for the zero γ 94 and for 38 zeros there were no sign changes at all. In the Fig. 3 we present plots of K l (m) as a function of m for a few zeros γ l. Let the rational p n(l) (γ l )/q n(l) (γ l ) be the n-th partial convergent of the continued fractions (3): p n(l) (γ l ) = a(l). (6) q n(l) (γ l ) For each zero γ l we have calculated the partial convergents p n(l) (γ l )/q n(l) (γ l ). Next from these denominators q n(l) (γ l ) we have calculated the quantities L l (n(l)): L l (n(l)) = ( q n(l) ) /n(l), l =, 2,..., 2600 (7) The obtained values of L l (n(l)) are presented in the Fig.4. These values scatter around the red line representing the Khinchin Lèvy s constant L 0. As in the case of K l (m) the Fig.5 presents the number of sign changes of the difference L l (m) L 0 as a function of the index m of the denominator of the m-th convergent p m /q m S L (l) = number of such m that (L l (m + ) L 0 )(L l (m) L 0 ) < 0. (8) The maximal number of sign changes was 36 and appeared for the zero γ 389 and there were 396 zeros without sign changes. In the Fig. 6 we have plotted the running absolute difference between K l (m) and K 0 averaged over all 2600 zeros: A K (m) = 2600 2600 K l (m) K 0, m = 00, 0,... 788 (9) l= and the similar average for the difference between L l (m) and L 0 : A L (m) = 2600 2600 Ll (m) L 0, m = 00, 0,... 788. (20) l= These two averages very rapidly tend to zero. Although it does not prove nothing, the fact that the curves representing A K (m) and A L (m) almost coincide is very convincing. In the inset the plot on double logarithmic scale reveals that both A L (m) and A k (m) decrease like C K,L / m where C K = 2.3868... and C L = 2.4473.... It is a pure speculation linking the power of these dependencies m /2 to the ordinate of the critical line. 3 Concluding remarks There are generalizations of above quantities K(n) given by (9). It can be shown that the following s-mean values of the denominators a k (r) of the continued fraction 5
for a real number r: M(n, s; r) = ( n n ( ak (r) ) ) /s s (2) are divergent for s and convergent for s < for almost all real r [4,.8]. It can be shown that for s < ( ) lim M(n, s; r) = k s log 2 ( ) /s K n (k + ) 2 s (22) k= where M(n, s; r) for almost all r are the same. The quantities (2) can be computed for imaginary parts of nontrivial zeta zeros M(n, s; γ l ) and compared with values of K s but we leave it for further investigation. The continued fraction were used in the past in the Apery s proof of irrationality of ζ(3). In the paper [3] values of ζ(n) for all n 2 were expressed in terms of rapidly converging continued fractions. These results were analytical, but in case of the nontrivial zeros of the ζ(s) function we are left only with the computer experiments. The results reported in this paper suggests that they are irrational. k= References [] B. Conrey, H. M. Bui, and M. P. Young. More than 4% of the zeros of the zeta function are on the critical line math.nt/002.427, 200. http: //xxx.lanl.gov/abs/math/002.427v. [2] J. B. Conrey. More than two-fifths of the zeros of the Riemann zeta function are on the critical line. J. Reine Angew. Math., 399: 26, 989. [3] D. Cvijovi and J. Klinowski. Continued-fraction expansions for the riemann zeta function and polylogarithms. Proceedings of the American Mathematical Society, 25(9):2543 2550, 997. [4] S. Finch. Mathematical Constants. Cambridge University Press, 2003. [5] X. Gourdon. The 0 3 first zeros of the Riemann Zeta Function, and zeros computation at very large height, Oct. 24, 2004. http://numbers.computation. free.fr/constants/miscellaneous/zetazerose3-e24.pdf. [6] G. H. Hardy. Sur les zéros de la fonction ζ(s) de Riemann. C. R. Acad. Sci. Paris, 58:02 04, 94. [7] A. E. Ingham. On two conjectures in the theory of numbers. Amer. J. Math., 64:33 39, 942. [8] A. Y. Khinchin. Continued Fractions. Dover Publications, New York, 997. 6
[9] N. Levinson. More than one third of zeros of Riemann s zeta-function are on σ = /2. Advances in Math., 3:383 436, 974. [0] A. M. Odlyzko. The 0 20 -th zero of the Riemann zeta function and 75 million of its neighbors. 992 revision of 989 manuscript. [] A. M. Odlyzko. The 0 2 -st zero of the Riemann zeta function. Nov. 998 note for the informal proceedings of the Sept. 998 conference on the zeta function at the Edwin Schroedinger Institute in Vienna. [2] A. M. Odlyzko. Tables of zeros of the Riemann zeta function. http://www. dtc.umn.edu/~odlyzko/zeta_tables/index.html. [3] A. M. Odlyzko. Primes, quantum chaos, and computers. Number Theory, National Research Council, pages 35 46, 990. [4] A. M. Odlyzko. The 0 22 -nd zero of the Riemann zeta function. In M. van Frankenhuysen and M. L. Lapidus, editors, Dynamical, Spectral, and Arithmetic Zeta Functions, number 290 in Amer. Math. Soc., Contemporary Math. series, pages 39 44, 200. [5] A. M. Odlyzko and H. J. J. te Riele. Disproof of the Mertens Conjecture. J. Reine Angew. Math., 357:38 60, 985. [6] B. Riemann. Ueber die anzahl der primzahlen unter einer gegebenen grösse. Monatsberichte der Berliner Akademie, pages 67 680, November 859. [7] C. Ryll-Nardzewski. On the ergodic theorems II (Ergodic theory of continued fractions). Studia Mathematica, 2:7479, 95. [8] A. Selberg. On the zeros of Riemann s zeta-function. Skr. Norske Vid. Akad. Oslo I., 0:59, 942. [9] The PARI Group, Bordeaux. PARI/GP, version 2.3.2, 2008. available from http://pari.math.u-bordeaux.fr/. [20] E. W. Weisstein. CRC Concise Encyclopedia of Mathematics. Chapman & Hall/CRC, 2009. 7
Fig. The plot of K l (n(l)). 8
Fig.2 The number of such m that (K l (m + ) K 0 )(K l (m) K 0 ) < 0 for each l (the initial transient values of m were skipped sign changes were detected for m = 00, 0,... n(l)). 9
Fig.3 For γ 66 there was no sign change of the difference K 66 (m) K 0. For γ 94 there were 22 sign changes of the difference K 94 (m) K 0 it was the largest number of sign changes among all zeros. For γ 434 there were 63 sign changes of the difference K 434 (m) K 0. 0
Fig.4 The plot of L l (n(l)).
Fig.5 The number of such m that (L l (m + ) L 0 )(L l (m) L 0 ) < 0 for each l (the initial transient values of m were skipped sign changes were detected for m = 00, 0,... n(l)). 2
Fig.6 The averaged over all 2600 zeros differences K l (m) K 0 and L l (m) L 0 plotted for m = 00, 0,... 788). In the inset the same curves are plotted on the double logarithmic scale together with fits obtained by the least square method. The equations of the fits are 2.3868/m 0.509 for A K (m) (green line) and 2.4473/m 0.5028 for A L (m) (blue line) 3