PART 1: The Riemann Zeta function ζ() [RZF]
- 1.
-
As defined in literature (Sondow et al, Weisstein, Edwards, Jekel)
- 1.1.
Figure 1.
Riemann Zeta function in R.
Figure 1.
Riemann Zeta function in R.
-
- 1.2.
Euler Product Formula that links
with the distribution of prime numbers
- 1.3.
Integral definition of
:
Where , is the Gamma function
- 1.4.
Analytical continuation of for :
-k<Re(s) [Kopp, Konrad. 1945]:
- 1.5.
Laurent series of
:
where
are the Stieltjes constants.
- 1.6.
Hurwitz function
:
- 1.7.
Generalized Harmonic Function
:
- 1.8.
ζ(s) converges for s>1 to the following values (Sloane):
Table 1.
Values of ζ(s).
s |
ζ(s) |
Known ζ(s) representations over π |
2 |
1.6449179 |
π2/6 |
4 |
1.0823232 |
π4/90 |
6 |
1.0173431 |
π6/945 |
8 |
1.0040774 |
π8/9450 |
- 1.9.
An approximation for the values of for s>1 in R
Based on this expression, for s sufficiently large, one can represent
as a multiple of
:
with a very good approximation given by:
The error between the calculated and actual is very small for s>4.
Some calculated values of
calculated and
actual:
Table 2.
Values of calculated and actual.
Table 2.
Values of calculated and actual.
One can use [2] to propose the following approximation for
:
Table 3.
Comparing and CZ(s).
Table 3.
Comparing and CZ(s).
|
s=3 |
s=4 |
s=10 |
s=14 |
|
1.20206 |
1.0823 |
1.000994 |
1.0000612 |
|
1.18659 |
1.0784 |
1.000988 |
1.0000611 |
Graphically:
Figure 2.
Caceres’ approximation for the Riemann Zeta function in R.
Figure 2.
Caceres’ approximation for the Riemann Zeta function in R.
-
- 1.10.
-
Zeros of (Weisstein)
- 1.10.1.
There are no zeros of . converges for s>1.
- 1.10.2.
is the Harmonic function and diverges, and:
(the Euler-Mascheroni constant)
- 1.10.3.
has trivial zeros at , with natural numbers
has nontrivial zeros at at the critical line (Brentt)
The Riemann Hypothesis [RH] establishes that all nontrivial zeros of zeta have List of imaginary part of the first nontrivial zeros of zeta (Odlyzko):
ß(1) = 14.134725142
ß(3) = 25.010857580
ß(4) = 30.424876126
ß(6) = 37.586178159
ß(7) = 40.918719012
ß(9) = 48.005150881
ß(10) = 49.773832478
Representation of |:
Figure 3.
Modulus of the Riemann Zeta function |
Figure 3.
Modulus of the Riemann Zeta function |
PART 3: A decomposition of ζ(z) based on the C-transformation of . for
The C-values of
from [6] and [7] are equal to the ζ(
) when Re(z)>1, this error |
- ζ(
)| grows significantly in the critical strip
as observed in the following table:
Table 5.
Values of for 0Re(z)<1 for N=500.
Table 5.
Values of for 0Re(z)<1 for N=500.
A |
ß |
|
ζ
(
)
|
ζ
(
)|
|
0.0 |
0 |
CN{f} for N=500
|
-0.5 |
0.5 |
0.2 |
2 |
0.399824505+0.322650799 i |
0.360103 + 0.266246 i |
|
0.7 |
0 |
-2.777900606 |
-2.7783884455 |
|
To understand better the value of the difference - ζ(), one can plot the difference for : (Similar exponential charts occur for all values of for any given value of ß)
Figure 4a.
Where a=Re(z) and b=Im(z).
Figure 4a.
Where a=Re(z) and b=Im(z).
And plot the difference for variable values of : (Similar sine charts occur for all values of for any given value of α)
Figure 4b.
Where a=Re(z) and b=Im(z).
Figure 4b.
Where a=Re(z) and b=Im(z).
These charts lead to the following calculation of the difference
:
With when .
And one can finally write:
and the C-value of for is the Riemann Zeta function ζ(z).
- 2.
-
Decomposition of
One can rewrite [8] and [9] creating the and functions:
where:
and define:
to write:
The following table compared the values of
and
:
Table 6.
ζ(z) compared to X(z) - Y(z).
Table 6.
ζ(z) compared to X(z) - Y(z).
z= 0 +j* 0 and n=500 |
Zeta(z) = -0.5 + i* 0.0 X(z)-Y(z) = -0.5 +i* 0.0 ---> Error = 0.0 +i* 0.0 |
z= 0.2 +j* 2 and n=500 |
Zeta(z) = 0.360102590022591 + i* -0.266246199765574 X(z)-Y(z) = 0.360102741838091 +i* -0.266246128959438 ---> Error= -1.5181550 e-7 +i* -7.080613 e-8 |
z= 0.4 +j* 0 and n=500 |
Zeta(z) = -1.13479778386698 + i* 0.0 X(z)-Y(z) = -1.1347977871726 +i* 0.0 ---> Error= 3.305619 e-9 +i* 0.0 |
The highest error for αn=1000 is .
- 3.
Representation of the function for Re(z)=1/2
Figure 5.
ζ(z) = X(z) - Y(z).
Figure 5.
ζ(z) = X(z) - Y(z).
- 4.
Representation of the function for Re(z)=1/2
Figure 6.
|ζ(z)| = |X(z) - Y(z)|.
Figure 6.
|ζ(z)| = |X(z) - Y(z)|.
- 5.
-
Representation of the function )
The following chart represents for and and .
- 6.
-
Representation of the function
The following chart represents for and and
- 7.
-
Representation of |
Wave representation for . [Figure. 11]
- 8.
-
Representation of |
Polynomial representation for . [Figure. 14]
- 9.
Conclusion PART 3
Using the defined C-transformation, one can write the Riemann Zeta function as the difference of two functions X(z) and Y(z) for . This will provide a new way of analyzing the zeros of the Zeta function, and a new approach to the Riemann Hypothesis.
The decomposition is as follows:
, where:
and:
and:
Observations:
has a wave representation
becomes a parable when z is a nontrivial zero of Riemann Zeta
becomes a line when z is a nontrivial zero of RZF with slope equal
has a polynomial representation
becomes a parable when z is a nontrivial zero of Riemann Zeta
becomes a line when Re(z)=1/2 with slope equal
So, the only common representation for |X(z)| and |Y(z)| occurs when Re(z)=1/2, so
= 0 if and only if Re(z)=1/2
PART 4: Proof of the Riemann Hypothesis using the decomposition.
- 1.
-
Analysis of Absolute Square
This could be observed in Figure. 14, 15, 16.
1.1. is a straight line if and only if
The slope of
with respect to
n is given by:
can only be a line when the slope is constant, which can only happen if and only if:
therefore:
1.2. Summary for for :
- ⇨
the slope is constant if and only if
- ⇨
When α=1/2, =
- ⇨
The slope for is for
- 2.
-
Analysis of Absolute Square
Applying properties of infinite series (Kopp):
One can express the previous expression replacing:
When is represented graphically, one can observe that:
- -
is a wave that converges when
and α>1 (
Figure 17)
- -
is a wave that does not converge when
and α<1 (
Figure 18)
- -
is a wave that collapses to a line when
and α=1/2 and ß=Im(
) (
Figure 19) where z* is a nontrivial zero of RZF.
2.1. converges when and α>1 to
The limit of
outside the critical strip [0,1] can be calculated using [16]:
As one can see in some examples in the following table where z=α+iß:
Table 7.
α |
ß |
|
|
1.0 |
7 |
1.074711506185445 |
1.074756 |
1.0 |
10 |
1.4413521753699579 |
1.441430 |
2.5 |
7 |
1.0093487944300192 |
1.009349 |
2.5 |
10 |
1.0507402208589398 |
1.050740 |
One can observe that the graphs for α=1 do not converge while graphs for α>1 they all converge.
This observation can be used to prove that there are no zero values of ζ(z) for z with Re(z)>1.
2.2. diverges when for α1
diverges when
for α<1 because of [16] and [17]:
2.3. does not collapse to any polynomial function
One can prove it with a reduction to absurd.
Let’s assume that
where C and t integers C>0 and t>0
So, must be 0 which is an absurd.
2.4. collapses to a straight-line
The proposition says that the following limit exists only for Re(z) = 1/2;
2.4.1. For α>1/2, one can see that /= 0:
= 0 because 2α>1 and the series is convergent
2.4.2. For α<1/2, one can see that/=as:
Where the summations are replaced by the number of elements in the matrix (n x n) times the smallest value in each row (1/n) then 1>(2-1-2α)>0 when α<1/2
2.4.3. Limit for α=1/2.
When α=1/2, one can express (
/n
) as:
Using the integral approximation of the infinite series
So, if
exists will be equal to:
if z=1/2+iß
And this limit can only exist when
is monotonous which means that the curve will cross the x-axis only once.
- 3.
-
Calculating the zeros of
Let’s define the function
in R (where z=a+bi) such that:
With the following wave representation for :
Figure 21.
for a=0.4 and variable b.
Figure 21.
for a=0.4 and variable b.
Figure 22.
for a=0.5 and variable b.
Figure 22.
for a=0.5 and variable b.
Figure 23.
for a=0.6 and variable b.
Figure 23.
for a=0.6 and variable b.
As a wave, can have one or more zeros. For to have only one zero, it must cross the axis y=0 only once, which means that the wave collapses to a polynomial line. A numeric method has been created and coded (Python) to find the values of such that =0. The following table shows an example of those calculated values, where x=n, a=Alfa, and b=Beta:
Table 8.
Number of Zeros of for different values of a=Alfa, and b=Beta.
Table 8.
Number of Zeros of for different values of a=Alfa, and b=Beta.
The calculations for and only found single zeros for for values of as shown in the following table that summarizes the single zeros found in those intervals:
Table 9.
List of first Zeros of .
Table 9.
List of first Zeros of .
And the calculated values of
for the values of (a,b) from
Table 9 are:
Table 10.
Limit of
for b in
Table 10 and x->∞.
Table 10.
Limit of
for b in
Table 10 and x->∞.
Table 11.
Comparing “b” calculated with known zeros of ζ(z).
Table 11.
Comparing “b” calculated with known zeros of ζ(z).
Therefore, = C has the following special properties for all (a,b) such that ζ(a+bi)=0.
if S=
Graphically:
Figure 24.
such that ζ(1/2+b*i)=0.
Figure 24.
such that ζ(1/2+b*i)=0.
- 4.
-
Theorem: For Re(z)≥0, if z* is a nontrivial zero of ζ(z), then Re(z*)=1/2
Proof:
- ⮚
From [10], [11], [12]: ζ(z) = X(z) − Y(z) for Re(z)>0, z≠1
- ⮚
From [13]: |Y(z, n)|2 is always a polynomial line.
- ⮚
From [14]: |Y(z, n)|2 is only straight line if and only if Re(z)= ½
= tends to a straight line with slope
- ⮚
From [15]: |X(z, n)|2 is a wave function that has only one polynomial representation in the form of a straight line if and only if Re(z)= ½ [18] and for certain values of Im(z)=ß* calculated using [19]. These values of ß* coincide with the imaginary parts of the nontrivial zeros of Riemann Zeta z*, so:
= tends to a straight line with slope
when Re(z)=1/2 and ß=NTZ of RZF
- ⮚
Therefore and only occur when Re(z*) = ½
- ⮚
As ζ(z)=X(z)-Y(z), therefore all zeros of ζ(z) for z>=0, z≠1 have Re(z)=1/2. [QED]
Figure 25.
for ζ(z)=X(z)-Y(z)=0 -> |X(z)| = |Y(z)| for Re(z) = 1/2,.
Figure 25.
for ζ(z)=X(z)-Y(z)=0 -> |X(z)| = |Y(z)| for Re(z) = 1/2,.
PART 5: On the distribution of the zeroes of the RZF in the critical line
- 1.
-
From [17] one can write:
therefore, the limit of
From [18], if
exists will be equal to:
- 2.
Calculating the nontrivial zeros of ζ(z) using the Harmonic function
From the previous equations, and for any
, a nontrivial zero of Zeta in the critical line α=1/2, one can write:
Where
is the Harmonic function. One can simplify the expression by creating functions
:
From the definition of limit, one can write that for any ε arbitrarily small, there exists and
N such that for any
n>N:
If
, then [20] can be written as:
The following chart shows the representation of H(n), O(n), and P(n) [P(n) is a straight line with slope ]:
Figure 26.
Straight Lines P(n).
Figure 26.
Straight Lines P(n).
The equation [20] can be used to create an algorithm to find the nontrivial zeros of zeta in the critical line without knowing any of them based on their connection to the Harmonic function.
An example of a Python code to calculate the zeros of zeta in the critical line with 1 decimal places accuracy based on [20]:
# __Pedro Caceres__ 2020 Feb 17
#Rough code to find zeros of Riemann Zeta using the Harmonic function
print('Harmonic(',nn,')=', harmo)
a1 = nn/((1-alfa)**2 + b**2)
b1 += (k*j)**(-alfa) * m.cos(b * m.log(k/j))
if abs(h1-harmo) < epsilon:
print('------> Solution beta=',b, ' ... and->', h1-harmo)
This code tends the following results:
Harmonic( 50 )= 4.4792053383294235
--------------> Solution beta= 14.1 ... and error -> 0.0067952158225219605
--------------> Solution beta= 25.0 ... and error -> -0.008460202279115592
--------------> Solution beta= 30.4 ... and error -> 0.0024237587453344034
--------------> Solution beta= 37.6 ... and error -> 0.0012958863904977136
--------------> Solution beta= 40.9 ... and error -> -0.009083573623293262
--------------> Solution beta= 48.0 ... and error -> -0.0027214317425938717
--------------> Solution beta= 49.6 ... and error -> 0.0024275253143217768
These values compared to (Odlyzko):
ß(1) = 14.134725142
ß(3) = 25.010857580
ß(4) = 30.424876126
ß(6) = 37.586178159
ß(7) = 40.918719012
ß(9) = 48.005150881
ß(10) = 49.773832478
Changing the values of “n” and epsilon, one can increase the accuracy of the results.
The fact that the Harmonic function, Hn , can be expressed in an infinite number of ways as a function of any ß=Im(z) imaginary part of a nontrivial solution of ζ(z), provides also an algorithm to calculate all nontrivial zeros from any known zero.
Let’s define the function:
For α=1/2, and ε arbitrarily small, for any two nontrivial zeros of zeta
and
, there exists and N such that for any n>N:
This proposition means that the nontrivial zeros of the Riemann Zeta are not distributed randomly, and they follow a defined structure.
Sample code to show how [21] can be used to find zeros based on a known zero:
# Code to find zeros from any known zero
# __Pedro Caceres__ 2020 Feb 17
nn = 60 #Not really high. Used for a rough calculation
#Calculating H(1/2,zero,n) = a - b
a2 = nn/((1-alfa)**2 + zero**2)
b2 += (k * j) ** (-alfa) * m.cos(zero * m.log(k / j))
h2 = a2-b2 #H2 to compare against
# range to find additional zeros of zeta
for b in range (245000,310000): #adding digits increases accuracy
a1 = nn/((1-alfa)**2 + b**2)
b1 += (k*j)**(-alfa) * m.cos(b * m.log(k/j))
#If error < epsilon, then print potential zero
print('-----------> Solution beta=',b, ' ... and error ->', h1-h2)
Results:
-----> Solution beta= 25.0155 ... and error -> +1.442262027140373 e-05
-----> Solution beta= 30.4385 ... and error -> -1.140533215249206 e-05
These values compared to (Odlyzko):
ß(3) = 25.010857580
ß(4) = 30.424876126
Changing the values of the variable “nn” and epsilon in the code, the accuracy can be increased to more decimal places.
- 3.
Conclusion
The distribution of the nontrivial zeros of the Riemann Zeta function in the critical line is not random. They are located in values of
that verify that for any ß, and ε arbitrarily small, there exists and
N such that for any
n>N:
Nomenclature and conventions
- –
ζ(z)= is the Riemann Zeta Function (RZF).
- –
z* is any nontrivial zero (NTZ) of the RZF verifying that ζ(z*) = 0.
- –
ß*(n) is the nth zero of the Riemann function in the critical line Re(z)=1/2 in C
- –
α=Re(z) is the real part of z
- –
ß=Im(z) is the imaginary part of z
- –
If z=α+iß, Modulus(z)= and SquareAbsolute(z)
References
- B. Riemann, “Ueber die Anzahl der Primzahlen unter einer gegebenen Grӧsse”, Monat. der Kӧnigl. Preuss. Akad. der Wissen. zu Berlin aus der Jahre 1859 (1860), 671–680; also, Gesammelte math. Werke und wissensch. Nachlass, 2. Aufl. 1892, 145–155.
- Brent, R. P. ``On the Zeros of the Riemann Zeta Function in the Critical Strip.'' Math. Comput. 33, 1361-1372, 1979. [CrossRef]
- Edwards , Harold M. Riemann’s Zeta Function. 1974. Dover Publications.
- Encyclopedia Britannica, Editors, 2008, https://www.britannica.com/science/Riemann-zeta-function, accessed March 2018.
- Jekel, David. The Riemann Zeta Function. 1963. University of Washington. https://sites.math.washington.edu/~morrow/336_13/papers/david.pdf.
- Kopp, Konrad. “Theory and Application of Infinite Series.” Blackie & Son Limited, London, 1928.
- Odlyzko, A. “The 10^20 th Zero of the Riemann Zeta Function and 70 Million of Its Neighbors.”.
- Python Software Foundation. Python Language Reference, version 3.7. Available at http://www.python.org.
- Sloane, N. J. A. Sequence A002432. An On-Line Version of the Encyclopedia of Integer Sequences. http://oeis.org/search?q=A002432&language=english&go=Search.
- Sondow, Jonathan and Weisstein, Eric W. "Riemann Zeta Function." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/RiemannZetaFunction.html.
- Weisstein, Eric W. "Riemann Zeta Function Zeros." From MathWorld--A Wolfram Web Resource. https://mathworld.wolfram.com/RiemannZetaFunctionZeros.html.
- Weisstein, Eric, 1999. “Riemann Zeta Function”, https://archive.lib.msu.edu/crcmath/math/math/r/r316.htm.
- Weisstein, Eric, 1999. “Riemann Hypothesis”, https://archive.lib.msu.edu/crcmath/math/math/r/r294.htm.
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).