A Proof of the Riemann Hypothesis Based on a New Expression of ξ(s)

1. Introduction

The Riemann zeta function is originally defined in the half-plane $\Re (s)>1$ by the absolutely convergent series

$$\zeta (s)=\sum _{n=1}^{\infty }{\displaystyle \frac{1}{n^s}},\Re (s)>1$$

(1)

The connection between the above-defined Riemann zeta function and prime numbers was discovered by Euler, i.e., the famous Euler product

$$\zeta (s)=\sum _{n=1}^{\infty }{\displaystyle \frac{1}{n^s}}=\prod _p{(1-p^{-s})}^{-1},\Re (s)>1$$

(2)

where p runs over the prime numbers.

Riemann showed in his paper in 1859 how to extend the zeta function to the whole complex plane $\mathbb{C}$ by analytic continuation [1]

$$\zeta (s)={\displaystyle \frac{\Gamma (1-s)}{2\pi i}}{\int }_{\infty }^{\infty }{\displaystyle \frac{{(-x)}^s}{e^x-1}}·{\displaystyle \frac{dx}{x}}$$

(3a)

where $"{\int }_{\infty }^{\infty }"$ is the symbol adopted by Riemann to represent the contour integral from $+\infty$ to $+\infty$ around a domain which includes the value 0 but no other point of discontinuity of the integrand in its interior.

Or equivalently,

$$\zeta (s)={\displaystyle \frac{{\pi }^{s/2}}{\Gamma (s/2)}}\{{\displaystyle \frac{1}{s(s-1)}}+{\int }_1^{\infty }(x^{{\displaystyle \frac{s}{2}}-1}+x^{-{\displaystyle \frac{s}{2}}-{\displaystyle \frac{1}{2}}})·({\displaystyle \frac{\theta (x)-1}{2}})dx\}$$

(3b)

where $\theta (x)={\sum }_{-\infty }^{\infty }{e}^{-n^2\pi x}$ is the Jaccobi theta function, $\Gamma$ is the Gamma function in the following Weierstrass expression

$${\displaystyle \frac{1}{\Gamma (s)}}=s·e^{\gamma s}\prod _{n=1}^{\infty }(1+{\displaystyle \frac{s}{n}})e^{-s/n}$$

(4)

where $\gamma$ is the Euler-Mascheroni constant.

As shown by Riemann, $\zeta (s)$ extends to $\mathbb{C}$ as a meromorphic function with only a simple pole at $s=1$, with residue 1, and satisfies the following functional equation

$${\pi }^{-{\displaystyle \frac{s}{2}}}\Gamma ({\displaystyle \frac{s}{2}})\zeta (s)={\pi }^{-{\displaystyle \frac{1-s}{2}}}\Gamma ({\displaystyle \frac{1-s}{2}})\zeta (1-s)$$

(5)

The Riemann zeta function $\zeta (s)$ has zeros at the negative even integers: $-2$, $-4$, $-6$, $-8$, … and one refers to them as the trivial zeros. The other zeros of $\zeta (s)$ are complex numbers, i.e., non-trivial zeros.

In 1896, Hadamard [2] and Poussin [3] independently proved that no zeros could lie on the line $\Re (s)=1$, together with the functional equation $\xi (s)=\xi (1-s)$ and the fact that there are no zeros with real part greater than 1, this showed that all non-trivial zeros must lie in the interior of the critical strip$0<\Re (s)<1$. Later on, Hardy (1914) [4], Hardy and Littlewood (1921) [5] showed that there are infinitely many zeros on the critical line $\Re (s)={\displaystyle \frac{1}{2}}$.

To give a summary of the related research works on the RH, we have the following results on the properties of the non-trivial zeros of $\zeta (s)$ [2–7].

Lemma 1: Non-trivial zeroes of $\zeta (s)$, noted as $\rho =\alpha +j\beta$, have the following properties

1) The number of non-trivial zeroes is infinity;

2) $\beta \ne 0$;

3) $0<\alpha <1$;

4) $\rho ,\overline{\rho },1-\overline{\rho },1-\rho$ are all non-trivial zeroes.

As further study, the completed zeta function $\xi (s)$ is proposed, i.e.

$$\xi (s)={\displaystyle \frac{1}{2}}s(s-1){\pi }^{-{\displaystyle \frac{s}{2}}}\Gamma ({\displaystyle \frac{s}{2}})\zeta (s)$$

(6)

It is well-known that $\xi (s)$ is an entire function of order 1. This implies $\xi (s)$ is analytic, and can be expressed as infinite product of polynomial factors, in the whole complex plane $\mathbb{C}$. In addition, replacing s with $1-s$ in Eq.(6), and combining Eq.(5), we obtain the following functional equation

$$\xi (s)=\xi (1-s)$$

(7)

According to the definition of $\xi (s)$, and recalling Eq.(4), the trivial zeros of $\zeta (s)$ are canceled by the poles of $\Gamma ({\displaystyle \frac{s}{2}})$. The zero of $s-1$ and the pole of $\zeta (s)$ cancel; the zero $s=0$ and the pole of $\Gamma ({\displaystyle \frac{s}{2}})$ cancel [7–9]. Thus, all the zeros of $\xi (s)$ are exactly the nontrivial zeros of $\zeta (s)$. Then we have the following Lemma 2.

Lemma 2: The zeros of $\xi (s)$ coincide with the non-trivial zeros of $\zeta (s)$.

Consequently, the following two statements are equivalent.

Statement 1: All the non-trivial zeros of $\zeta (s)$ have real part equal to ${\displaystyle \frac{1}{2}}$.

Statement 2: All zeros of $\xi (s)$ have real part equal to ${\displaystyle \frac{1}{2}}$.

To prove the RH, a natural thinking is to estimate the numbers of non-trivial zeros of $\zeta (s)$ inside or outside some certain areas according to Argument Principle. Along this train of thought, there are many research works. Let $N(T)$ denote the number of non-trivial zeros of $\zeta (s)$ inside the rectangle: $0<\alpha <1,0<\beta \le T$, and let $N_0(T)$ denote the number of non-trivial zeros of $\zeta (s)$ on the line $\alpha ={\displaystyle \frac{1}{2}},0<\beta \le T$. Selberg proved that there exist positive constants c and $T_0$, such that $N_0(T)>cN(T),(T>T_0)$ [10], later on, Levinson proved that $c\ge {\displaystyle \frac{1}{3}}$ [11], Lou and Yao proved that $c\ge 0.3484$ [12], Conrey proved that $c\ge {\displaystyle \frac{2}{5}}$ [13], Bui, Conrey and Young proved that $c\ge 0.41$ [14], Feng proved that $c\ge 0.4128$ [15], Wu proved that $c\ge 0.4172$ [16].

On the other hand, many non-trivial zeros have been calculated by hand or by computer programs. Among others, Riemann found the first three non-trivial zeros [17]. Gram found the first 15 zeros based on Euler-Maclaurin summation [18]. Titchmarsh calculated the 138^th to 195^th zeros using the Riemann-Siegel formula [19,20]. Here are the first three (pairs of) non-trivial zeros: ${\displaystyle \frac{1}{2}}\pm j14.1347251;{\displaystyle \frac{1}{2}}\pm j21.0220396;{\displaystyle \frac{1}{2}}\pm j25.0108575$.

The idea of this paper originates from Euler’s work on proving the famous equality

$$1+{\displaystyle \frac{1}{2^2}}+{\displaystyle \frac{1}{3^2}}+{\displaystyle \frac{1}{4^2}}+{\displaystyle \frac{1}{5^2}}+\dots ={\displaystyle \frac{{\pi }^2}{6}}$$

(8)

This result was deduced by comparing the coefficients of two infinite expressions of ${\displaystyle \frac{sinx}{x}}$: one as a power series and the other as an infinite product,

$$\begin{array}{cc}\hfill {\displaystyle \frac{sinx}{x}}& =1-{\displaystyle \frac{x^2}{3!}}+{\displaystyle \frac{x^4}{5!}}-{\displaystyle \frac{x^6}{7!}}+\dots =(1-{\displaystyle \frac{x^2}{{\pi }^2}})(1-{\displaystyle \frac{x^2}{4{\pi }^2}})(1-{\displaystyle \frac{x^2}{9{\pi }^2}})\dots \hfill \end{array}$$

(9)

Motivated by this approach, we conjecture that $\xi (s)$ can be factored into the form $(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})$, which is verified by pairing ${\rho }_i$ and ${\overline{\rho }}_i$ in the Hadamard product representation of $\xi (s)$, i.e. $(1-{\displaystyle \frac{s}{{\rho }_i}})(1-{\displaystyle \frac{s}{{\overline{\rho }}_i}})={\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})$.

The Hadamard product expansion of $\xi (s)$, first proposed by Riemann and later rigorously justified by Hadamard [21], is given by

$$\xi (s)=\xi (0)\prod _{\rho }(1-{\displaystyle \frac{s}{\rho }}),s\in \mathbb{C}$$

(10)

where $\xi (0)={\displaystyle \frac{1}{2}}$, $\rho$ runs over all zeros of $\xi (s)$.

Hadamard showed that to ensure the absolute convergence of this infinite product expansion, $\rho$ and $1-\rho$ must be paired. Later in Section 4, we will demonstrate that pairing $\rho$ with its complex conjugate $\overline{\rho }$ can also be used to ensure the absolute convergence.

2. Preliminary Lemmas

This section provides preliminary lemmas supporting the proof of the key lemma - Lemma 8 in the next section.

The key point of this section is to extend, by employing the divisibility concept of entire functions, both the transitivity of polynomial divisibility and a property of irreducible polynomials (Lemma 3) to the context of infinite products of polynomial factors.

We begin with the ring of real polynomials $\mathbb{R}\left[x\right]$, defined as

$$\mathbb{R}\left[x\right]=\{\sum _{i=0}^{\infty }{a}_i{x}^i|a_i\in \mathbb{R},a_i\ne 0\mathrm{for}\mathrm{all}\mathrm{but}\mathrm{a}\mathrm{finite}\mathrm{number}\mathrm{of}i\}$$

and equipped with the operations + (addition) and · (multiplication).

The ring of real polynomials is a subset of the ring of entire functions, which is defined as the set of all holomorphic functions on the whole complex plane $\mathbb{C}$, together with the operations of addition and multiplication, denoted as $\mathbb{H}(\mathbb{C})$ [22,23].

Both rings possess properties of divisibility, coprimality, and the greatest common divisor, denoted as "gcd". There are also differences between these two rings. Among others, the polynomial ring is a unique factorization domain (UFD), while the ring of entire functions is not a UFD. For entire functions, their divisibility, coprimality and common factors are determined by the relationships between their zero sets [23,24].

To facilitate the subsequent discussions (particularly the proof of Lemma 8), we provide the following definition.

Definition 2.1: Let $f(x)={\prod }_{i=1}^{\infty }{p}_i(x),p_i(x)\in \mathbb{R}\left[x\right]$, be an entire function, and $h(x)\in \mathbb{R}\left[x\right]$. We say $h(x)$ divides $f(x)$, denoted as $h(x)\mid f(x)$, if there exists an entire function $g(x)={\prod }_{i=1}^{\infty }{q}_i(x),q_i(x)\in \mathbb{R}\left[x\right]$, such that $f(x)=h(x)·g(x)$.

To support the proof of the key lemma - Lemma 8 in next section. We also need the following lemmas.

Lemma 3: Let $m(x),g_1(x),...,g_n(x)\in \mathbb{R}\left[x\right],n\ge 2$. If $m(x)$ is irreducible (prime) and divides the product $g_1(x)\dots g_n(x)$, then $m(x)$ divides one of the polynomials $g_1(x),...,g_n(x)$.

Lemma 4: Let $f(x),m(x)\in \mathbb{R}\left[x\right]$. If $m(x)$ is irreducible and $f(x)$ is any polynomial, then either $m(x)$ divides $f(x)$ or $gcd(m(x),f(x))=1$.

Lemma 5: Let $f(x)$ be an entire function on $\mathbb{C}$. Suppose $f(x)={\prod }_{i=1}^{\infty }{g}_i(x)$ converges absolutely on $\mathbb{C}$ and uniformly on every compact subset of $\mathbb{C}$, where each $g_i(x)\in \mathbb{R}\left[x\right]$ is irreducible in $\mathbb{R}\left[x\right]$ of degree $d\in \{1,2\}$. If $m(x)\in \mathbb{R}\left[x\right]$ is irreducible in $\mathbb{R}\left[x\right]$ of degree d and divides $f(x)$, then $m(x)$ divides $g_i(x)$ for some $i\ge 1$.

Remark: The contents of Lemma 3 and Lemma 4 can be found in many textbooks of linear algebra, modern algebra, or abstract algebra, see for example Refs.[25-27].

Below we give the proof of Lemma 5.

Proof: Let $\alpha$ be a root of $m(x)$, i.e., $m(\alpha )=0$. Since $m(x)\mid f(x)$, we have $f(\alpha )=0$. By uniform convergence of ${\prod }_{i=1}^{\infty }{g}_i(x)$, there exists at least one index $i\in \mathbb{N}$ such that $g_i(\alpha )=0$, otherwise $f(\alpha )\ne 0$. See Theorem 2 in Ref.[8] on page 178 for more details (also listed in Appendix for convenience), with $f_n(x)$ therein corresponding to ${\prod }_{i=1}^n{g}_i(x)$.

As $g_i(x)$ and $m(x)$ are irreducible over $\mathbb{R}$ with $deg(g_i(x))=deg(m(x))=d$, they share the root $\alpha$. Thus:

• If $d=1$, then $g_i(x)=a(x-\alpha )$ and $m(x)=b(x-\alpha )$ for $a,b\ne 0$, so $m(x)\mid g_i(x)$.
• If $d=2$, then both $g_i(x)$ and $m(x)$ have roots $\{\alpha ,\overline{\alpha }\}$, so $g_i(x)=c·m(x)$ for $c\ne 0$, hence $m(x)\mid g_i(x)$.

In both cases, $m(x)$ divides $g_i(x)$.

That completes the proof of Lemma 5.

Remark: In the proof of Lemma 5, the absolute convergence of ${\prod }_{i=1}^{\infty }{g}_i(x)$ is not directly used; it is only required to justify that the infinite product defines an entire function.

Additionally, we also need the following results on properties of zeros of entire function for understanding the multiplicity of zeros of $\xi (s)$.

Lemma 6: Let $f(s)$ be a non-zero entire function, and let $s_0$ be a zero of $f(s)$. Then the multiplicity of $s_0$ is a finite positive integer.

Proof: Let $f(s)\neg \equiv 0,s\in \mathbb{C}$, be an entire function, which means it is holomorphic on the whole complex plane. Suppose $f(s)$ has a zero at $s_0\in \mathbb{C}$ of multiplicity m, then $f(s)={(s-s_0)}^{m}g(s)$, where $g(s)$ is also an entire function and $g(s_0)\ne 0$.

Assume for contradiction that m is infinite, which implies there exists an accumulation point of zeros in the neighbor of $s_0$. Then, by Identity Theorem for holomorphic functions, and considering "0" is also an entire function, we have $f(s)\equiv 0,s\in \mathbb{C}$, which contradicts the given condition that $f(s)\neg \equiv 0,s\in \mathbb{C}$. Thus, the assumption is false, i.e., m must be a finite positive integer.

That completes the proof of Lemma 6.

Remark: Statements similar to Lemma 6 can be found in Ref.[28] and other related textbooks/monographs.

Lemma 7: Let $f(s)$ be a non-zero entire function, and let $s_0$ be a zero of $f(s)$. Then the multiplicity of $s_0$ is unique.

Proof: Let $f(s)\neg \equiv 0,s\in \mathbb{C}$, be an entire function, which has a multiple zero at $s_0\in \mathbb{C}$ of multiplicity m. We can write: $f(s)={(s-s_0)}^{m}g(s)$, where $g(s)$ is also an entire function and $g(s_0)\ne 0$.

Assume for contradiction that there exists another integer $n\ne m$ such that n is also a multiplicity of the zero $s_0$. This means we can also write: $f(s)={(s-s_0)}^{n}h(s)$, where $h(s)$ is an entire function and $h(s_0)\ne 0$.

Since both expressions for $f(s)$ must be equal, we then obtain ${(s-s_0)}^{m}g(s)={(s-s_0)}^{n}h(s)$. Without loss of generality, consider $m>n$, then we have: ${(s-s_0)}^{m-n}g(s)=h(s)\Rightarrow h(s_0)=0$, which is a contradiction to $h(s_0)\ne 0$. Thus, the assumption is false, i.e., the multiplicity of a zero of any non-zero entire function is unique.

That completes the proof of Lemma 7.

3. Key Lemma

In this section, we prove the key lemma - Lemma 8, which is substantial for the proof of the RH.

Lemma 8: Given that the entire function

$$f(s)=\prod _{i=1}^{\infty }{\left(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)}^{m_i},s\in \mathbb{C}$$

(11)

converges absolutely on $\mathbb{C}$ and uniformly on every compact subset of $\mathbb{C}$, where $0<{\alpha }_i<1$ and ${\beta }_i\ne 0$ are real numbers with $|{\beta }_1|\le |{\beta }_2|\le \dots$, $m_i\ge 1$ is the multiplicity of zero conjugates ${\rho }_i={\alpha }_i+j{\beta }_i,{\overline{\rho }}_i={\alpha }_i-j{\beta }_i$, ${\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{{\beta }_i^2}}<\infty$.

Under the substitution $s\mapsto 1-s$, we have

$$f(1-s)=\prod _{i=1}^{\infty }{\left(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)}^{m_i},s\in \mathbb{C}$$

(12)

Then

$$f(s)=f(1-s)⟺{\alpha }_i={\displaystyle \frac{1}{2}};i=1,2,3,\dots ;0<|{\beta }_1|<|{\beta }_2|<\dots$$

(13)

Remark: It should be noted that $m_i$ is actually the multiplicity of quadruplets of zeros $({\rho }_i,{\overline{\rho }}_i,1-{\rho }_i,1-{\overline{\rho }}_i)$ under the assumption $f(s)=f(1-s)$.

Proof: Considering

$$(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}{)}^{m_i}=\underset{m_i\mathrm{times}}{\underbrace{(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})\dots }}$$

(14)

we have from $f(s)=f(1-s)$ and Eq.(11)

$$\prod _{i=1}^{\infty }\left(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)=\prod _{i=1}^{\infty }\left(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)$$

(15)

where the $i^{th}$ factor appears $m_i$ times in any order.

Due to absolute convergence, both products ${\prod }_{i=1}^{\infty }\left(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)$ and ${\prod }_{i=1}^{\infty }\left(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)$ can be rearranged into any single quadratic polynomial factor times the remaining product, which stays absolutely convergent and hence defines an entire function. Then by the definition of divisibility for entire functions [23,24], more specifically Definition 2.1, Eq.(15) implies that every polynomial factor on one side of the equation divides the infinite product on the other side, i.e.,

$$\left\{\begin{array}{cc}& (1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})\mid \prod _{i=1}^{\infty }\left(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)\hfill \\ \hfill & (1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}})\mid \prod _{i=1}^{\infty }\left(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}\right)\hfill \end{array}\right.$$

(16)

where "∣" is the divisible sign, $i\in \mathbb{N}=\{1,2,3,\dots \}$.

Both polynomials $1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}$ and $1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}$ have discriminant $\Delta =-4·{\displaystyle \frac{1}{{\beta }_i^2}}<0$. Hence they are irreducible in $\mathbb{R}\left[x\right]$. By Lemma 5, Eq.(16) yields:

$$\left\{\begin{array}{cc}& (1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})\mid (1+{\displaystyle \frac{{(1-s-{\alpha }_l)}^2}{{\beta }_l^2}})\hfill \\ \hfill & (1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}})\mid (1+{\displaystyle \frac{{(s-{\alpha }_l)}^2}{{\beta }_l^2}})\hfill \\ \hfill & i,l\in \mathbb{N}\hfill \end{array}\right.$$

(17)

For the special kind of polynomials in Eq.(17), "divisible" means "equal", which can be verified by comparing the like terms in the following equation to get $k=1$.

$$(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})=k(1+{\displaystyle \frac{{(1-s-{\alpha }_l)}^2}{{\beta }_l^2}}),k\ne 0\in \mathbb{R}$$

(18)

Further, due to the uniqueness of the multiplicity $m_i$, the only solution to Eq.(17) is $i=l$, otherwise, duplicated zeros (in quadruplets) with ${\alpha }_i+{\alpha }_l=1,{\beta }_i^2={\beta }_l^2,l\ne i$ would be generated to change $m_i$. Therefore we have from Eq.(17):

$$(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})=(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}),i\in \mathbb{N}$$

(19)

By comparing the like terms in Eq.(19), we obtain ${\alpha }_i={\displaystyle \frac{1}{2}}(\forall i)$. Further, to ensure the uniqueness of $m_i$ while ${\alpha }_i={\displaystyle \frac{1}{2}}$, we need limit the ${\beta }_i$ values to be distinct, i.e., $0<|{\beta }_1|<|{\beta }_2|<|{\beta }_3|<\dots$. Thus, we obtain

$$f(s)=f(1-s)⟹{\alpha }_i={\displaystyle \frac{1}{2}};i=1,2,3,\dots ;0<|{\beta }_1|<|{\beta }_2|<\dots$$

(20)

On the other hand, we have the following fact.

$$\begin{array}{cc}\hfill & {\alpha }_i={\displaystyle \frac{1}{2}},i=1,2,3,\dots ,0<|{\beta }_1|<|{\beta }_2|<\dots \hfill \\ \hfill & \Rightarrow (\mathrm{considering}{\beta }_i\ne 0)\hfill \\ \hfill & (1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})=(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}})\hfill \\ \hfill & \Rightarrow (\mathrm{considering}{m}_i\ge 1)\hfill \\ \hfill & (1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}{)}^{m_i}=(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}{)}^{m_i}\hfill \\ \hfill & \Rightarrow (\mathrm{taking}\mathrm{infinite}\mathrm{products}\mathrm{on}\mathrm{both}\mathrm{sides}\mathrm{of}\mathrm{the}\mathrm{above}\mathrm{equations})\hfill \\ \hfill & \prod _{i=1}^{\infty }(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}}{)}^{m_i}=\prod _{i=1}^{\infty }(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}}{)}^{m_i}\hfill \end{array}$$

(21)

From Eqs.(20) and (21), we finally get:

$$f(s)=f(1-s)⟺{\alpha }_i={\displaystyle \frac{1}{2}};i=1,2,3,\dots ;0<|{\beta }_1|<|{\beta }_2|<\dots$$

That completes the proof of Lemma 8.

In addition, Lemma 9 will also be used in the proof of the RH in the next section.

Lemma 9: The infinite product ${\prod }_{i=1}^{\infty }({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}{)}^{m_i}$ converges to a non-zero constant, given the conditions: $0<{\alpha }_i<1,{\beta }_i\ne 0,{\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{{\beta }_i^2}}<\infty$, and $m_i\ge 1$ is the multiplicity of zero ${\alpha }_i+j{\beta }_i$.

Proof: First of all, we know that

$$\prod _{i=1}^{\infty }({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}{)}^{m_i}=\prod _{i=1}^{\infty }{\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}$$

where in the right side expression, $i^{th}$ factor appears $m_i$ times.

Let $a_i={\displaystyle \frac{{\alpha }_i^2}{{\alpha }_i^2+{\beta }_i^2}}$, then ${\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}=1-{\displaystyle \frac{{\alpha }_i^2}{{\alpha }_i^2+{\beta }_i^2}}=1-a_i$.

Since $0<{\alpha }_i<1$ and ${\beta }_i\ne 0$, we have: $0<a_i<{\displaystyle \frac{1}{{\beta }_i^2}}$. Then ${\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{{\beta }_i^2}}<\infty$ (given condition) implies ${\sum }_{i=1}^{\infty }\left|a_i\right|={\sum }_{i=1}^{\infty }{a}_i<\infty$ (absolute convergence).

Further, the absolute convergence of ${\sum }_{i=1}^{\infty }{a}_i$ guarantees that the product ${\prod }_{i=1}^{\infty }(1-a_i)={\prod }_{i=1}^{\infty }{\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}={\prod }_{i=1}^{\infty }({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}{)}^{m_i}$ converges to a non-zero constant.

That completes the proof of Lemma 9.

4. A Proof of the RH

This section presents a proof of the Riemann Hypothesis. We first prove that Statement 2 of the RH is true, and then by Lemma 2, Statement 1 of the RH is also true. To be brief, to prove the Riemann Hypothesis, it suffices to show that ${\alpha }_i={\displaystyle \frac{1}{2}},i=1,2,3,\dots$ in the new expression of $\xi (s)$ as shown in Eq.(25).

Proof of the RH: The proof proceeds in three steps.

Step 1:

It is well-known that zeros of $\xi (s)$ always come in complex conjugate pairs. Then by pairing ${\rho }_i={\alpha }_i+j{\beta }_i$ and ${\overline{\rho }}_i={\alpha }_i-j{\beta }_i$ in the Hadamard product as shown in Eq.(10), we have

$$\begin{array}{cc}\hfill \xi (s)& =\xi (0)\prod _{\rho }(1-{\displaystyle \frac{s}{\rho }})=\xi (0)\prod _{i=1}^{\infty }(1-{\displaystyle \frac{s}{{\rho }_i}})(1-{\displaystyle \frac{s}{{\overline{\rho }}_i}})\hfill \\ \hfill & =\xi (0)\prod _{i=1}^{\infty }(1-{\displaystyle \frac{s}{{\alpha }_i+j{\beta }_i}})(1-{\displaystyle \frac{s}{{\alpha }_i-j{\beta }_i}})=\xi (0)\prod _{i=1}^{\infty }({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}})\hfill \end{array}$$

(22)

where $0<{\alpha }_i<1,{\beta }_i\ne 0$ (according to Lemma 1) are real numbers with $|{\beta }_1|\le |{\beta }_2|\le |{\beta }_3|\le \dots$.

The absolute convergence of the infinite product in Eq.(22) in the form

$$\xi (s)=\xi (0)\prod _{i=1}^{\infty }(1-{\displaystyle \frac{s}{{\rho }_i}})(1-{\displaystyle \frac{s}{{\overline{\rho }}_i}})=\xi (0)\prod _{i=1}^{\infty }(1-{\displaystyle \frac{s(2{\alpha }_i-s)}{|{\rho }_i{|}^2}})$$

(23)

depends on the convergence of infinite series ${\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{|{\rho }_i{|}^2}}$ (since $\left|s\right|<\infty$ and $0<{\alpha }_i<1\Rightarrow \left|s(2{\alpha }_i-s)\right|<\infty$), which is a well-known fact according to Theorem 2 in Section 2, Chapter IV of Ref.[9] (also listed in Appendix for convenience). Thus, the infinite products as shown in Eq.(23) and Eq.(22) are absolutely convergent on the whole complex plane $\mathbb{C}$. According to Weierstrass M test, these two infinite products are also uniformly convergent on any compact subset of $\mathbb{C}$ (Let $K\subset \mathbb{C}$ be an arbitrary compact set, then there exists constant $M_K>0$ such that $\forall s\in K,{\displaystyle \frac{|s(2{\alpha }_i-s)|}{|{\rho }_i{|}^2}}\le {\displaystyle \frac{M_K}{|{\rho }_i{|}^2}},{\sum }_{i=1}^{\infty }{\displaystyle \frac{M_K}{|{\rho }_i{|}^2}}=M_K{\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{|{\rho }_i{|}^2}}<\infty$).

Further, considering the absolute convergence of

$$\xi (s)=\xi (0)\prod _{i=1}^{\infty }(1-{\displaystyle \frac{s(2{\alpha }_i-s)}{|{\rho }_i{|}^2}})=\xi (0)\prod _{i=1}^{\infty }({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}})$$

(24)

we have the following new expression of $\xi (s)$ by putting all the possible multiple factors (zeros) together:

$$\xi (s)=\xi (0)\prod _{i=1}^{\infty }({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}}{)}^{m_i}$$

(25)

where $m_i\ge 1$ is the multiplicity of ${\rho }_i/{\overline{\rho }}_i$, $i=1,2,3,\dots$.

Step 2: Replacing s with $1-s$ in Eq.(25), we obtain the infinite product expression of $\xi (1-s)$, i.e.,

$$\xi (1-s)=\xi (0)\prod _{i=1}^{\infty }{({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}})}^{m_i}$$

(26)

where $m_i\ge 1$ is the multiplicity of $1-{\rho }_i/1-{\overline{\rho }}_i$, $i=1,2,3,\dots$.

Step 3: According to the functional equation $\xi (s)=\xi (1-s)$, and considering Eq.(25) and Eq.(26), we have

$$\xi (0)\prod _{i=1}^{\infty }{({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}})}^{m_i}=\xi (0)\prod _{i=1}^{\infty }{({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}})}^{m_i}$$

(27)

which is equivalent to

$$\prod _{i=1}^{\infty }{(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})}^{m_i}=\prod _{i=1}^{\infty }{(1+{\displaystyle \frac{{(1-s-{\alpha }_i)}^2}{{\beta }_i^2}})}^{m_i}$$

(28)

We know that $0<{\alpha }_i<1,|{\rho }_i{|}^2\to \infty (\sin \mathrm{ce}{\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{|{\rho }_i{|}^2}}<\infty ,\mathrm{then}{\displaystyle \frac{1}{|{\rho }_i{|}^2}}\to 0)\Rightarrow {\left|{\beta }_i\right|}^2\to \infty$.

Then we have ${lim}_{i\to \infty }{\displaystyle \frac{{\beta }_i^2}{|{\rho }_i{|}^2}}={lim}_{i\to \infty }{\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}=1$, that means ${\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{{\beta }_i^2}}$ and ${\sum }_{i=1}^{\infty }{\displaystyle \frac{1}{|{\rho }_i{|}^2}}$ have the same convergence. Thus, both sides of Eq.(28) are absolutely and uniformly convergent, with reference to the analysis of Eqs. (22) and (23) in Step 1.

Furthermore, both sides of Eq.(28) converge to entire functions, because they differ with the entire functions $\xi (s)$ and $\xi (1-s)$ respectively by a non-zero multiplicative constant, i.e.

$$\begin{array}{cc}\hfill \xi (s)=& \xi (0)\prod _{i=1}^{\infty }{({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}}+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\alpha }_i^2+{\beta }_i^2}})}^{m_i}\hfill \\ \hfill =& \xi (0)\prod _{i=1}^{\infty }{({\displaystyle \frac{{\beta }_i^2}{{\alpha }_i^2+{\beta }_i^2}})}^{m_i}\prod _{i=1}^{\infty }{(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})}^{m_i}\hfill \\ \hfill =& c·\prod _{i=1}^{\infty }{(1+{\displaystyle \frac{{(s-{\alpha }_i)}^2}{{\beta }_i^2}})}^{m_i}\hfill \end{array}$$

(29)

where c is a non-zero constant, see Lemma 9 for details.

Finally, according to Lemma 8, Eq.(28) is equivalent to

$${\alpha }_i={\displaystyle \frac{1}{2}};0<|{\beta }_1|<|{\beta }_2|<|{\beta }_3|<\dots ;i=1,2,3,\dots$$

(30)

Thus, we conclude that all the zeros of the completed zeta function $\xi (s)$ have real part equal to ${\displaystyle \frac{1}{2}}$, i.e., Statement 2 of the RH is true. According to Lemma 2, Statement 1 of the RH is also true, i.e., all the non-trivial zeros of the Riemann zeta function $\zeta (s)$ have real part equal to ${\displaystyle \frac{1}{2}}$.

That completes the proof of the RH.