Asymptotic harmonic behavior in the prime number distribution
We consider on , where the sum is over all primes . If is bounded on , then the Riemann hypothesis is true or there are infinitely many zeros Re . The first 21 zeros give rise to asymptotic harmonic behavior in defined by the prime numbers up to one trillion.
The Riemann-zeta function is the analytic extension of
where Euler’s identity on the right hand side expresses the relation of the integers to the primes. The zeros of Riemann’s analytic continuation of (1) comprise the negative even integers, , and an infinite number of nontrivial zeros in the strip .
A general approach to find zeros is by continuation (Keller, 1987). If is a starting point of a path with tangent ,
then the endpoint is a zero of , all of which are isolated. All known nontrivial zeros satisfy to within numerical precision, the first three of which are , By the symmetry
Continuation (2) is determined by the prime numbers, since
The poles of at the zeros are therefore expressed by the prime number distribution.
In this paper, we study the distribution of zeros by Fourier analysis of the function
on , where
with summation over all primes. In what follows, we put
The are absolutely summable by Stirling’s formula and the asymptotic distribution of .
Theorem 1.1. In the limit as becomes small, we have the asymptotic behavior
In (9), is evidently unbounded in the limit as approaches zero whenever a finite number of zeros exists off the critical line Re .
Corollary 1.2. If is bounded, then the Riemann hypothesis is true or there are infinitely many zeros Re .
based on the Chebyshev functions
where the sum is over all primes and integers . In (9), has a normalization by according to and is absolutely convergent for all , whereas in (10) is normalized by and the sum is not absolutely convergent. Similar to Corollary 1.2, the left hand side of (10) will be bounded in the limit of large if the Riemann hypothesis is true.
§2 presents some preliminaries on . §3 gives an integral representation of and a discussion on its singularity at . In §4, Cauchy’s integral formula is applied to derive a sum of residues associated with the . The proof Theorem 1.1 follows from a Fourier transform and asymptotic analysis (§5). In §6, we illustrate a direct evaluation of using the primes up to one trillion, showing harmonic behavior arising from by the first few zeros . We summarize our findings in §7.
Riemann obtained an analytic extension of by expressing in terms of ,
Here, satisfies as approaches zero by the identity for the Jacobi function . 111When is an integer, is one-half the surface area of . On Re, it obtains the meromorphic expression (e.g. Borwein et al. 2006)
which gives a maximal analytic continuation of and shows a simple pole at with residue 1.
Riemann further introduced the symmetric form , satisfying , whereby
in terms of the digamma function
in the limit of large .
Lemma 2.1. In the limit of large , the logarithmic derivative of satisfies
Lemma 2.2. Along the line , we have the asymptotic expansion in the limit of large , whereby the are absolutely summable.
Proof. Recall (8) and the asymptotic expansion with a branch cut along the negative real axis. In the limit of large , , and hence , since as becomes large. Hence, the are absolutely summable. Numerically, their sum is small, based on a large number of known zeros .
3 An integral representation of
Following the same steps leading to the Riemann integral for , we have
where absorbs the simple pole in at due to the simple pole in at , leaving analytic at . Following a decomposition ,
and substitution , appears as the Laplace transforms
These integral expressions allow continuations to , respectively, the entire complex plane.
Lemma 3.1. Analytic extension of extends to .
Proof. With , the second term on the right hand side in (5) satisfies
where is analytic at . By (22), as approaches from the right, we have
where is analytic about .
Fig. 2 shows a numerical evaluation of for small evaluated for the 37.6 billion primes up to one trillion, allowing down to () in view of the requirement for an accurate truncation in as defined by (7). The result shows asymptotic harmonic behavior in the limit as becomes small.
If the integral
is absolutely convergent as approaches zero, e.g., when is of one sign in some neighborhood of , as in the numerical evaluation shown in Fig. 2, then has an analytic extension into Re with no singularities, implying the absence of in this region. However, this requires information on the point wise behavior of , which goes beyond the relatively weaker integrability property (23).
To make a step in this direction, we next apply a linear transform to (5) to derive the asymptotic behavior of in terms of the distribution .
4 A sum of residues associated with the non-trivial zeros
and its Fourier transform
Lemma 4.1. has a simple pole at with residue 1 and simple poles at each of the nontrivial zeros of with residue .
Proof. We have (e.g. Borwein et al. 2006)
where is a constant, so that
where denotes the digamma function as before, includes contributions from the logarithmic derivative of the factor to in (31), whose singularities are restricted to the trivial zeros of .
We now consider the Fourier integral over Re as part of contour integration closed over and Re .
Proposition 4.2. The Fourier transform of over Re satisfies
in the limit of large .
Proof. Integration over gives
where we choose to be between two consecutive values of . We have
In the limit as approaches infinity, approaches zero and becomes small by Lemma 2.2., whence
Next, integration over with a small semicircle around obtains an result in the limit of large by application of Lemma 2.1-3 and the Riemann-Lebesque Lemma. The result now follows in the limit as approaches infinity, taking into account the residue sum associated with the and absolute summability of the .
5 Proof of Theorem 1.1
Multiplying (5) by , we have
We thus consider
which ab initio is defined on Re by Euler’s identity with Fourier transform
Turning to the right hand side of (40), we consider the coefficients
Here, since . In particular, and has a well defined limit and in the limit as becomes arbitrarily large.
Lemma 5.1. The sum is well-defined on Re .
Proof. The result follows from the case . By the Prime Number Theorem, , whereby summation over the tails satisfy
whenever . Hence, for , whenever . It follows that
on Re .
Lemma 5.2. For any , the Fourier transform of over Re satisfies
Proof. The Fourier integral can be obtained in a contour integration with closure over and the edges for large . In the notation (42), it obtains a residue at , since , whence
The integral (46) exists by virtue of a removable singularity of at . It asymptotically decays to zero for large when by the Riemann-Lebesque Lemma.
Lemma 5.3. For , the Fourier transform
in the limit of large .
Proof. Since is analytic in Re , we are at liberty to consider the transform on . The result follows from the Riemann-Lebesque Lemma.
Proof of Theorem 1.1. The Fourier transform of (47) is
By Proposition 4.2 and Lemmas 5.1-5.2, we have
With , Theorem 1.1 now follows.
6 Numerical illustration of asymptotic harmonic behavior
The harmonic behavior emerges in
To search for higher harmonics associated with the zeros in , we compare the spectrum of by taking a Fast Fourier Transform with respect to ,
and compare the results with an analytic expression for the Fourier coefficients of the ,
where denotes the Bessel function of the first of order . Fig. 3 shows the first 21 harmonics in our evaluation of , which is about the maximum that can be calculated by direct summation in quad precision.
The zeros of the Riemann-zeta function are endpoints of continuation, defined by an expressed by a regularized sum over the prime numbers defined by (6).
The zeros of introduce asymptotic harmonic behavior in as a function of , defined by the sum of residues of the , shown in Figs. 2-3. Primes up to 4 billion are needed to identify the first 4 harmonics, up to 70 billion for the 10 and up to 1 trillion for the first 21. It appears that the prime number range scales approximately exponentially with the number of harmonics it contains.
Theorem 1.1 describes a correlation between the distribution of the primes and the distribution of the nontrivial zeros . Suppose there are a finite number of zeros in Re . We may then consider for which gives rise to dominant exponential growth in in the limit as becomes large. This observation leads to Corollary 1.2. can remain bounded in only if the Riemann hypothesis is true, or if remains fortuitously bounded as an infinite sum over with no maximum in . Conversely, Riemann hypothesis implies
According to (9) and our numerical calculation shown in Fig. 3, the number of primes relevant to the observed asymptotic harmonic behavior scales approximately exponentially with the number of zeros . The zeros explored to large values by existing numerical experiments hereby constrain the distribution of an exponentially large number of primes.
Acknowledgment. The author gratefully acknowledges stimulating discussions with Fabian Ziltener and Anton F.P. van Putten. Some of the manuscript was prepared at the Korea Institute for Advanced Study, Dongdaemun-Gu, Seoul. This research was supported in part by the National Science Foundation through TeraGrid resources provided by Purdue University under grant number TG-DMS100033. We specifically acknowledge the assistance of Vicki Halberstadt, Rich Raymond and Kimberly Dillman. The computations have been carried out using Lahey Fortran 95.
- Borwein et al. (2006) Borwein, P., Choi, S., Rooney, B., & Weirathmueller, A., 2006, The Riemann Hypothesis (Berlin: Springer Verlag)
- Cheng (1999) Cheng, Y., 1999, An explicit upper bound for the Riemann Zeta function near the line , Rocky Mountain J. Math. 29, 115 - 140
- Dusart (1999) Dusart, P.. 1999, Inégalités explicites pour , , et les nombres premiers., C. R. Math. Rep. Acad. Sci. Canada, 21, 53 - 59
- Ford (2002) Ford, K., 2002, Zero-free regions for the Riemann zeta function, in Number Theory for the Millenium (Urbana, IL, 2000), 2, 25 - 26
- Hadamard (1893) Hadamard, J., 1893, Etude sur les propriétés des fonctions entires et en particulier d’une fonction, J. Math. Pures Appl., 171 - 216
- Keiper (1992) Keiper, J.B., Power Series Expansions of Riemann’s Function, Math. Comput., 58, 765 - 773
- Keller (1987) Keller, H.B., 1987, Numerical Methods in Bifurcation Problems (Berlin: Springer Verlag/Tata Institute for Fundamental Research)
- Lehmer (1988) Lehmer, D.H., 1988, The Sum of Like Powers of the Zeros of the Riemann Zeta Function, Math. Comput., 50, 265 - 273
- Littlewood (1922) Littlewood, J.E., 1922, Researches in the theory of the Riemann function, Proc. London Math. Soc., ser. 2, 20, xxii - xxvii
- Littlewood (1924) Littlewood, J.E., 1926, On the Riemann zeta-function, Proc. London Math. Soc., ser. 2, 24, 175 - 201
- Littlewood (1927) Littlewood, J.E., 1928, Mathematical Notes (5): On the function , Proc. London Math. Soc., ser. 2, 27, 349 - 357
- Richert (1967) Richert, H.E., 1967, Zur Abschätzung der Riemannschen Zetafunktion in der Nähe der Vertikalen , Math. Ann. 169, 97 - 101
- von Mangoldt (1895) von Mangoldt, H., Zu Riemann’s Abhandlung ‘Über die Anzahl der Priemzahlen unter einer gegebenen Grösse’, , J. Reine Angew. Math. 114, 255 - 305
- Titchmarsh (1986) Titchmarsh, E.C., 1986, The theory of the Riemann Zeta-function (Oxford, 2nd ed., revised by D. R. Heath-Brown)
- Wintner (1941) Wintner, A., 1941, On the asymptotic behavior of the Riemann zeta-Function on the Line , Amer. J. Math., 63, 575 - 580