A Direct Proof of the Prime Number Theorem using Riemann’s Prime-counting Function Zihao Liu International Department, The Affiliated High School of SCNU E-mail: travor [email protected] Abstract. In this paper, we develop a novel analytic method to prove the prime number theorem in de la Vall´ee Poussin’s form: c√log x π(x) = li(x)+ (xe− ) O Instead of performing asymptotic expansion on Chebyshev functions as in conventional analytic methods, this new approach uses contour-integration method to analyze Riemann’s prime counting function J(x), which only differs from π(x) by (√x log log x/ log x). O 1. Introduction The prime number theorem[1][3] has been a popular topic in analytic number theory since the 19th century. Its first proofs were independently given by Hadamard and de la Vall´ee Poussin in 1896 [13] using analytic methods. Since then, mathematicians such as Apostol [1], Levinson [5], Newman [9][7], Selberg [10], Stein [11], and Wright [14] have explored different approaches (e.g. elementary methods, Tauberian theorems, and contour integrations) to prove this theorem. Instead of analyzing π(x) directly, present proofs, regardless elementary or analytic, first attacked either of the Chebyshev functions ϑ(x)= log p (1) p x arXiv:2105.05317v3 [math.NT] 18 Jun 2021 X≤ ψ(x)= Λ(n)= ϑ(x1/k) (2) n x k log x X≤ ≤X2 to derive asymptotic formulae, and estimates for π(x) were later obtained using elementary methods, such as partial summation, to ϑ(x). Although Landau [1][4] derived the prime number theorem by proving 1 lim µ(n) = 0 (3) x x →∞ n x X≤ µ(n) lim = 0 (4) x n →∞ n x X≤ he still justified the equivalence of (3) and (4) to the prime number theorem with Chebyshev functions. Hence, it is an interesting task to investigate a method that avoids the use of (1) and (2). In this paper, we propose a proof of the prime number theorem that avoids the use of Chebyshev functions. In particular, we first introduce Riemann’s prime counting function J(x) to our problem and relate it to π(x). Then, by inverse Mellin transform, an integral representation is obtained for π(x). In the subsequent part of this paper, we evaluate the contour integral and obtain the prime number theorem in the version of de la Vall´ee Poussin. That is, there exists a positive constant c such that c√log x π(x) = li(x)+ (xe− ) (5) O where li(x) is the logarithmic integral: 1 ε x dt li(x) = lim − + (6) ε 0+ log t → Z0 Z1+ε In essence, the main contribution of this paper is a new proof of the prime number theorem. The significance of this proof is that it is an analytic method that avoids the use of Chebyshev functions. 2. Riemann’s prime counting function J(x) Let an be defined by 1/k n = pk,p prime an = (7) (0 otherwise In his paper on π(x), Riemann [2] defined J(x) as the summatory function for an: 1 J(x)= a = (8) n k n x 1<pk x X≤ X≤ This function has a remarkable property that allows us to connect it to the standard prime counting function π(x): Theorem 1. For x 2, we have ≥ J(x)= π(x)+ (√x log log x) (9) O Proof. It follows from (8) that π(x1/k) π(√x) J(x) π(x)= − k ≤ k 2 k log x 2 k log x ≤ X≤ 2 ≤ X≤ 2 √x √x log log x ≤ k ≪ 2 k log x ≤ X≤ 2 By (9), asymptotic formula for J(x) can be ported to π(x) with ease, so we can study π(x) simultaneously when we analyze J(x). 3. Integral representation for J(x) and π(x) Riemann [2] had shown that 1 k+i xs J(x)= ∞ log ζ(s)ds (10) 2πi k i s Z − ∞ for noninteger x > 2 and k > 1 using Fourier analysis. However, because (10) involves infinite integrals, it would be difficult to analyze its asymptotic behavior directly, so we wish to convert (10) into a quantitative form. In other words, we prove a special case of Perron’s formula [13][8]. Lemma 1. For (s) > 1, we have ℜ ∞ a n = log ζ(s) (11) ns n=1 X Proof. Since (s) > 1, the left hand side converges absolutely, we can change the order of summation safely:ℜ ∞ a ∞ 1 1 n = = log ns kpks 1 p s n=1 p p − X X Xk=1 X − 1 = log = log ζ(s) 1 p s p − Y − Lemma 2. For k> 1 and noninteger x> 2, if we set 1 k+iT xs J(x)= log ζ(s)ds + R1(x,k,T ) (12) 2πi k iT s Z − then xk ∞ 1 R (x,k,T ) (13) 1 ≪ T nk log x/n nX=1 | | Proof. If we define 0 0 <y< 1 h(y)= (14) (1 y > 1 then J(x) can be rewritten into ∞ x J(x)= a h (15) n n n=1 X The right hand side converges certainly because the summand vanishes once n>x. Due to Titchmarsh [13], for k> 1 we have 1 k+iT ys yk h(y)= ds + (16) 2πi k iT s O T log y Z − | | Since log ζ(s) converges absolutely whenever (s) = k > 1, we plug (16) and (11) into (15). Finally, a 1 allows us to deduce (13). ℜ | n|≤ Lemma 3. Using the same notation, we have xk x log x x R (x,k,T ) + + (17) 1 ≪ T (k 1) T T x N − | − | where N denotes the nearest integer to x. Proof. To handle the logarithms in the denominator in (13), we partition the sum into xk = + + + Tnk log x/N n n<x/2,n>2x x/2 n N 1 N+1 n 2x X X ≤X≤ − X≤ ≤ | | In the first sum, it is evident that log x/n > log 2, so we have | | xk 1 xk 1 xk ∞ 1 T nk log x/n ≪ T nk ≤ T nk n=1 n<x/X2,n>2x | | n<x/X2,n>2x X xk xk = ζ(k) T ≪ T (k 1) − in which the last inequality follows from the fact that ζ(s) 1/(s 1) when s 1 [13]. When x/2 n 2x, the logarithm function satisfies that∼ − → ≤ ≤ x x n log | − | (18) n ≫ x so we have xk 1 1 1 x 1 T nk log x/n ≪ T log x/n ≪ T x n x/2 n N 1 x/2 n N 1 x/2 n N 1 ≤X≤ − | | ≤X≤ − | | ≤X≤ − − x N 1 dt x x/2 x log x − = log ≪ T x t T x N + 1 ≪ T Zx/2 − − in which the last part follows from the fact that N 1/2 x N + 1/2. Due to symmetry, (18) can be applied to the sum over N + 1 n 2x−to obtain≤ the≤ same estimate. For the last part, using (18), we have ≤ ≤ xk 1 x (19) Tnk log x/N ≪ T log x/N ≪ T x N | | | | | − | Combining all these estimates gives (17). Theorem 2. For half-odd integer x> 2 and k =1+1/ log x, we have 1 k+iT xs π(x)= log ζ(s)ds + R2(x, T ) (20) 2πi k iT s Z − with x log x R (x, T ) + √x log log x (21) 2 ≪ T Proof. Plugging k =1+1/ log x into (17) gives x log x x R (x, 1 + 1/ log x, T ) + (22) 1 ≪ T T x N | − | when x is half-odd, we have x N = 1/2, so the second term of (22) can be removed. Plugging (22) into (12), we get | − | 1 k+iT xs x log x J(x)= log ζ(s)ds + (23) 2πi k iT s O T Z − Finally, we replace J(x) with (23) in (9) to obtain the result. 4. Evaluation of the integral To handle the integral, we consider a rectangular region R with vertices 1 δ iT and k iT . This implies − ± ± 1 k+iT xs 1 xs log ζ(s)ds = log ζ(s)ds (24) 2πi k iT s 2πi ∂R s Z − I 1 1 δ iT 1 δ+iT k+iT + − − + − + (25) 2πi k iT 1 δ iT 1 δ+iT Z − Z − − Z − where ∂R denotes the boundary of R in counterclockwise direction and δ is chosen such that ζ(s) =0in R and on ∂R. This means that the only singularity of the integrand in R is at s = 1. However,6 because s = 1 is a logarithmic singularity, circumvention is needed to evaluate it. Lemma 4. If we define 1 xs 1 f(r)= log ds (26) 2πi s s/r 1 I(r+) − where (r+) denotes any counterclockwise path that encloses only the singularity at s = r, then xr f ′(r)= (27) r for r = 0. 6 Proof. Since ∂ 1 ∂ 1 1 log = [log r log(s r)] = + ∂r s/r 1 ∂r − − r s r − − we have 1 xs 1 1 f ′(r)= + ds 2πi s r s r I(r+) − xs s r xr = lim − + 1 = s r s r r → Lemma 5. Using the same notation, for all η > 0 we have f(1 + iη) = li(x1+iη) iπ (28) − 1 iη f(1 iη) = li(x − )+ iπ (29) − Proof. When η> 0, because (26) implies lima f(a) = 0, we have →∞ 1+iη f(1 + iη)= f(1 + iη) f( + iη)= f ′(r)dr − −∞ +iη Z−∞ 1+iη xr (1+iη) log x eu = dr = du +iη r ( +iη) log x u Z−∞ Z −∞ u=r log x By Cauchy’s integral theorem [12],| we{z can deform} the path of integration of the last integral into a line segment connecting and ε < 0, a clockwise circular arc connecting ε and ε, and finally a line segment connecting−∞ ε and− (1 + ηi) log x.