AMApplied Mathematics2152-7385Scientific Research Publishing10.4236/am.2021.128050AM-111639ArticlesPhysics&Mathematics The Proof of Riemann Hypothesis, the Key to the Door Is the Periodicity JinliangWang1*Research Institute of ESMD Method and Its Applications, College of Science, Qingdao University of Technology, Qingdao, China1208202112087237652, August 202128, August 2021 31, August 2021© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The Riemann hypothesis is a well-known mathematical problem that has been in suspense for 162 years. Its difficulty lies in the fact that it is involved in an infinite integral which includes infinite series with complex variables. To detour this is in vain, since all the messages are hid in it. To unscramble them, there is a totally new idea, that is, the “periodicity”! By investigating the numerical approximate values of zero points, an explicit distribution law on the critical line was found. To accord with this, a periodic form for the real part of Xi function was constructed and rigidly proved. The Riemann hypothesis can be divided into three progressive propositions. The first proposition (the number of zero points in the critical strip satisfies a certain estimation) had been proved in 1905. The second proposition (the number of zero points on the critical line satisfies the same estimation as in the critical strip) is ever in suspense. It can be solved perfectly with the newly found “periodicity”. The third proposition (all the nontrivial zero points are on the critical line), that is, the Riemann hypothesis, is also true. The proof is a combination of the symmetry, monotonicity, periodicity of the Xi function and the extremum principle of the harmonic functions. It is the moment to draw full stop for this suspending problem.

Riemann Hypothesis Riemann Zeta Function Distribution Law of Zero Point Periodicity Monotonicity Extremum Principle
1. Introduction

The Riemann hypothesis is a well-known mathematical problem. It had ever been a global hot topic when British mathematician Michael F. Atiyah (1929- 2019) reported his proof in Heidelberg Laureate Forum on Sep. 24, 2018. Unfortunately, his approach does not work and now this problem is still in suspense. To crack it, just as reported by New Scientist in , “unless with a totally new idea.” To read through the monographs in    and the popular readings in  , one can find the difficulty of Riemann hypothesis lies in the complexity of the Zeta function. Essentially, it is involved in an infinite integral which includes infinite series with complex variables. To detour this infinite integral is in vain, since all the messages are hid in it. To unscramble them, in my brain many ideas appeared and died during the past three years. Among them there is only one which survives and could be seen as “a totally new idea”. This is the “periodicity”! The proof will be a combination of the symmetry, monotonicity, periodicity and extremum principle. Notice that the present article may become the terminator of this problem, a whole story will be told.

This problem was left by German mathematician Bernhard Riemann (1826- 1866) in 1859   . He found that the distribution of prime numbers among all natural numbers is very closely related to the analytic continuation of the following infinite series on complex plane ℂ :

ζ ( s ) = ∑ n = 1 ∞ 1 n s = 1 1 s + 1 2 s + 1 3 s + ⋯ (1)

with s = σ + i t , which is usually called Riemann Zeta function. This series converges only for the case σ > 1 , so all the mentioned ζ ( s ) below should be understood in the sense of analytic continuation. If there is a s 0 which satisfies ζ ( s 0 ) = 0 , then we call it one zero point of ζ ( s ) . This function has real zero points at the negative even integers − 2, − 4, − 6 , ⋯ and one refers to them as the trivial zero points. Relatively, the other complex zero points of it are called the nontrivial zero points.

Riemann hypothesis: All the nontrivial zero points of ζ ( s ) have real part σ = 1 / 2 .

Is it true? In the official millennium problem description , E. Bombieri had reviewed that, in 1986 the first 1.5 × 109 nontrivial zero points of ζ ( s ) (arranged by increasing positive imaginary part) had been checked by J. Lune et al with numerical approach in , and the result showed that they are simple and all possess real part σ = 1 / 2 . More zero points had been checked by the followers, the same thing occurred. Till now the known record was set by X. Gourdon in 2004 , 1013 zero points. So the Riemann hypothesis is very likely true. What lacks is the theoretical proof.

The analytic continuation for ∑ n = 1 ∞     n − s (defined on σ > 1 ) is not unique. Among them, the most popular one is given by:

ζ ( s ) = π s / 2 ξ ( s ) ( s − 1 ) Γ ( s / 2 + 1 ) , (2)

which is meromorphic on ℂ with a unique pole at s = 1 (its residual is 1)   . Here

Γ ( s / 2 + 1 ) = ∫ 0 ∞     x s / 2 e − x d x (3)

with a default requirement σ / 2 + 1 > 0 . Its analytic continuation to the negative direction is done by Γ ( s / 2 ) = Γ ( s / 2 + 1 ) / ( s / 2 ) step by step. In this sense, − 2, − 4, − 6 , ⋯ are the zero points of 1 / Γ ( s / 2 + 1 ) , and this is why they are the trivial zero points of ζ ( s ) . Naturally, the function ξ ( s ) satisfies the inverse relation which was written down by Riemann:

ξ ( s ) = ( s − 1 ) π − s 2 Γ ( s / 2 + 1 ) ζ ( s ) = s ( s − 1 ) 2 π − s 2 Γ ( s / 2 ) ζ ( s ) = s ( s − 1 ) 2 [ 1 s ( s − 1 ) + ∫ 1 ∞     ψ ( x ) ( x s 2 − 1 + x − s + 1 2 ) d x ] = 1 2 − s ( 1 − s ) 2 ∫ 1 ∞     ψ ( x ) ( x s 2 − 1 + x − s + 1 2 ) d x (4)

with ψ ( x ) = ∑ n = 1 ∞     e − π n 2 x , which satisfies a peculiar relation bellow (see the deduction in   ):

2 ψ ( x ) + 1 = x − 1 / 2 [ 2 ψ ( x − 1 ) + 1 ] . (5)

Notice that ψ ( x ) is positive and exponentially decreasing on the interval [ 1, ∞ ) , the integral in Equation (4) is finite and hence ξ ( s ) is analytic on ℂ . Particularly, the zero points of ξ ( s ) coincide with the nontrivial zero points of ζ ( s ) , and the exploring of Riemann hypothesis can be converted to considering ξ ( s ) . To substitutes by 1 − s the expression maintains unchanged, so ξ ( s ) = ξ ( 1 − s ) . This symmetric property implies the most important relation between the function ζ ( s ) and the critical line σ = 1 / 2 . Correspondingly, the complex region with 0 ≤ σ ≤ 1 is called the critical strip which includes all the zero points of ξ ( s ) .

In fact, ξ ( s ) has no zero points outside the critical strip. With the aid of Euler formula:

ζ ( s ) = ∑ n = 1 ∞     n − s = ∏ p ( 1 − p − s ) − 1 (6)

(where the last expression is a product respect to all the prime numbers p = 2,3,5, ⋯ ), one can easily checked that | ζ ( s ) | ≠ 0 , that is | ξ ( s ) | ≠ 0 , for the case σ > 1 (see  ). Furthermore, the relation ξ ( s ) = ξ ( 1 − s ) indicates that the same thing is true for σ < 0 .

We note that, as the critical strip concerned, the two lines σ = 1 and σ = 0 can be gotten rid of. As reviewed in  , the proof given by de la Vallée Poussion in 1899 for this is very complex. Yet, as our approach concerned, it doesn’t matter. Any chosen strip with 1 / 2 − a ≤ σ ≤ 1 / 2 + a and a ≥ 1 / 2 can fulfil the request.

To denote ξ ( σ + i t ) = U ( σ , t ) + i V ( σ , t ) , then it follows from the well-known Reflection Principle ξ ( s ) ¯ = ξ ( s ¯ ) of complex conjugate that U ( σ , t ) − i V ( σ , t ) = U ( σ , − t ) + i V ( σ , − t ) . Meanwhile, the relation ξ ( s ) = ξ ( 1 − s ) yields U ( σ , t ) + i V ( σ , t ) = U ( 1 − σ , − t ) + i V ( 1 − σ , − t ) . Based on these two relations, the symmetric properties of ξ ( s ) are clarified as:

Lemma 1. To denote ξ ( σ + i t ) = U ( σ , t ) + i V ( σ , t ) , then its real part and imaginary part separately satisfy the symmetries and anti-symmetries bellow:

U ( σ , − t ) = U ( σ , t ) ,   V ( σ , − t ) = − V ( σ , t ) ,

U ( 1 − σ , t ) = U ( σ , t ) ,   V ( 1 − σ , t ) = − V ( σ , t ) .

For the particular cases with t = 0 and σ = 1 / 2 , the relations for V read V ( σ ,0 ) = − V ( σ ,0 ) and V ( 1 / 2 , t ) = − V ( 1 / 2 , t ) . So V ( σ , 0 ) = V ( 1 / 2 , t ) = 0 . This indicates that the values of ξ ( σ + i t ) are real on the lines t = 0 and σ = 1 / 2 . This lemma can be understood as: U and V are symmetric and anti-sym- metric about the two lines t = 0 and σ = 1 / 2 , respectively. That is why most of the arguments are about the upper quarter strip with 1 / 2 ≤ σ ≤ 1 and t > 0 . Particularly, due to the direct relationship with the Riemann hypothesis, the line σ = 1 / 2 has drawn much attention.

It is enlightening to divide Riemann hypothesis into three progressive propositions (see , pages 22-23):

Proposition 1. In the region bounded by 0 ≤ σ ≤ 1 and 0 < t < T , the number of zero points of ξ ( s ) is about T / 2 π log ( T / 2 π e ) .

Proposition 2. On the critical line σ = 1 / 2 with 0 < t < T , the number of zero points of ξ ( s ) is also about T / 2 π log ( T / 2 π e ) .

Proposition 3. All the zero points of ξ ( s ) are on the critical line σ = 1 / 2 .

Till now only the first proposition had been proved. It was finished by German mathematician Hans von Mangoldt in 1905   . His estimation for the zero-points number of ξ ( s ) is

N ( T ) = T 2 π l o g T 2 π e + O ( l o g ( T ) ) , (7)

which is understood as taking the integer part. We note that, O ( l o g ( T ) ) does not imply the existence of a term C l o g ( T ) for some constant C. Its meaning is that | N ( T ) − ( T / 2 π ) log ( T / 2 π e ) | ≤ C l o g ( T ) . In fact, the previously mentioned checks are done according to this theoretical result, and the numerical approximate values for the first four zero points are

1 / 2 + 14.1347251 i ,   1 / 2 + 21.0220396 i ,

1 / 2 + 25.0108575 i ,   1 / 2 + 30.4248761 i .

On the aspect of theoretical study, in 1914 Hardy firstly proved that, there are infinitely many zero points of ξ ( s ) on the critical line σ = 1 / 2 (see the review in  ). Due to the efforts of Selberg , Levinson  and Conrey , the ratio was lifted step by step. Now the known optimal estimation is that more than 40% of zero points of ξ ( s ) are on the critical line σ = 1 / 2 (they are also simple ones). These indicate that Proposition 2 is far from settled. Our approach is a direct attack to it. In addition, as reviewed in , “All the zero points of ξ ( s ) are simple ones” already became an accompanying conjecture to Riemann hypothesis. Yet, under our approach, this is not a thing.

The above survey is about the known knowledge of Riemann hypothesis, and the next is about our new findings. We begin with investigating the distribution law of zero points on the critical line σ = 1 / 2 in Section 2; The third section is about the main content, which includes 3 subsections. In the last subsection the periodicity of U ( σ , t ) is detailed illustrated. During the argument processes of periodicity, the monotonicity of the related functions become the preconditions. With the aid of explicit expression of U ( σ , t ) , the distribution law of zero points is verified. This indicates Proposition 2 given by Riemann holds true. The proof of Proposition 3, that is the Riemann hypothesis, is given in Section 4. The conclusions and related remarks are given in the last.

2. Distribution Law for the Zero Points on the Critical Line

On the critical line σ = 1 / 2 , the values of ξ ( s ) are real (ensured by Lemma 1), that is, only the real part U ( 1 / 2 , t ) of ξ ( 1 / 2 + i t ) is left. To recall the expression of ξ ( s ) in (4), it follows from   that

U ( 1 / 2 , t ) = ( s − 1 ) π − s 2 Γ ( s 2 + 1 ) ζ ( s ) | s = 1 / 2 + i t = − 1 + 4 t 2 8 π − s 2 Γ ( s 2 ) ζ ( s ) | s = 1 / 2 + i t = − a ( t ) e i ϑ ( t ) ζ ( 1 2 + i t ) , (8)

where ( ℜ and ℑ mean the taking of real and imaginary parts, respectively)

a ( t ) = 1 + 4 t 2 8 π − 1 / 4 e x p [ ℜ l o g Γ ( 1 + 2 t i 4 ) ] , (9)

ϑ ( t ) = ℑ l o g Γ ( 1 + 2 t i 4 ) − t 2 l o g π = t 2 l o g t 2 π e − π 8 + 1 48 t + 7 5760 t 3 + ⋯ . (10)

Notice that a ( t ) > 0 , the zero points on the critical line are only determined by the sign shift of Z ( t ) = e i ϑ ( t ) ζ ( 1 / 2 + i t ) . Furthermore, as reviewed in , the inverse expression can be employed to make the estimation [the zero points of ζ ( 1 / 2 + i t ) coincide with that of ξ ( 1 / 2 + i t ) ]:

ζ ( 1 / 2 + i t ) = e − i ϑ ( t ) Z ( t ) = Z ( t ) c o s ϑ ( t ) − i Z ( t ) s i n ϑ ( t ) . (11)

According to this, for the first time, Gram calculated the first 15 zero points in 1903. His method is an estimation approach based on s i n ϑ ( t ) = 0 , which yields the known Gram points at ϑ ( t k ) = k π with an approximation:

t k 2 log t k 2 π e − π 8 = k π ,   k = 1 , 2 , ⋯ . (12)

Certainly, to estimate the location of a zero point, the change of Z ( t ) should be also considered. This adds the complexity. By the way, 29 years later, the approximate approach for Z ( t ) was improved. That is the known Riemann-Siegel formula, which was discovered by Siegel in 1932 among Riemann’s private papers (see  ).

Though the method given by Gram is so rough, it is enlightening. The distribution of zero points on the critical line σ = 1 / 2 may possess some kind of periodicity. In fact, for some given t = T , the argument ϑ ( T ) obtained on the critical line and the zero-points number N ( T ) in the region 0 ≤ σ ≤ 1 and 0 < t < T have close relationship:

ϑ ( T ) = T 2 l o g T 2 π e − π 8 + ⋯ ,

N ( T ) = T 2 π l o g T 2 π e + 7 8 + ⋯ .

Exactly, as reviewed in  (page 132-134), in 1918 Backlund had revealed that, for all T ≥ 2 ,

| N ( T ) − π − 1 ϑ ( T ) − 1 | < 4.350 + 0.137 l o g T + 0.443 log log T .

Enlightened by this, we have a new idea. Proposition 2 holds true provided that U ( 1 / 2 , t ) = ξ ( 1 / 2 + i t ) can be expressed in the following form:

U ( 1 / 2 , t ) = A ( t ) s i n ( t 2 l o g t 2 π e + α + ε ( t ) ) , (13)

where A ( t ) , α and ε ( t ) are the amplitude (>0), initial phase angle and small perturbation function. Certainly, for the case with t = 2 π we have ( t / 2 ) l o g ( t / 2 π e ) = − π . Hence, a suitable phase angle α should be chosen to accord with the fact that ξ ( 1 / 2 + 2 π i ) = U ( 1 / 2 ,2 π ) > 0 (since ξ ( 0 ) = − ζ ( 0 ) = 1 / 2 and there is no zero points for 0 ≤ t ≤ 14 , see , page 31). If U ( 1 / 2 , t ) is indeed in such a periodic form, the thing becomes very simple. All the zero points on the critical line σ = 1 / 2 can be solved one by one with the following formula:

t k 2 log t k 2 π e + α + ε ( t k ) = k π ,   k = 1 , 2 , ⋯ . (14)

By the way, to count the number of zero points it can be used in a reverse way, that is,

N = t N 2 π l o g t N 2 π e + α π + ε ( t N ) π ,

which accords with Proposition 2. Relative to Gram's approximation in Equation (12), this formula is closer to the truth. In the following we give some numerical evidences for this. The theoretical proof will be given in the next sections.

To neglect the perturbation term ε ( t k ) and make fitting with the known numerical approximate values t k * to the zero points (provided by Odlyzko in  ), it leads to the results in Figure 1. The first 50 zero points fitting yields an optimal phase angle α = 4 π / 3 , that is, the k -th zero point on the critical line is very close to the solution of

t k 2 log t k 2 π e + 4 3 π = k π ,   k = 1 , 2 , ⋯ . (15)

We call this newly found formula as the distribution law of zero points.

The check in Figure 2 with 10,000 zero points shows that the above distribution

law is also obeyed. The error is small. Though the value of t k * is up to 9877, except the 14-th and 64-th zero points, all the corresponding errors satisfy | t k − t k * | < 1 .

In addition to the perturbation term ε ( t k ) , these errors may also partly ascribe to the numerical scheme for t k * , though its calculation accuracy is up to 109. Not forget that, essentially, the calculation of ξ ( s ) is involved in an infinite integral which includes infinite series with complex variables. At least, the approximate values for the 14-th and 64-th zero points are questionable. To take the 14-th zero as an example. t 14 * = 60.831778525 , t 15 * = 65.112544048 and hence

t 15 * 2 l o g t 15 * 2 π e − t 14 * 2 l o g t 14 * 2 π e ≈ 43.5681 − 38.6353 = 4.9328 > π . (16)

Yet, it follows from N ( T ) ≈ π − 1 ( T / 2 ) l o g ( T / 2 π e ) that, to add one π means to add a zero point. This indicates that 4.9328 is a too big difference between two adjacent zero points, it should be understood as a calculation error. With this understanding, the departure from the newly-found distribution law caused by the perturbation term ε ( t k ) may be smaller than the difference of | t k * − t k | depicted here.

Further check is given in Figure 3. This law is also obeyed for the zero points numbered from 99,900 to 100,000.

3. To Prove the Periodicity of U ( σ , t )

Can U ( 1 / 2 , t ) = ξ ( 1 / 2 + i t ) be expressed in the periodic form of Equation (13)? It needs a theoretical proof. Our theme is not limited by this, and in order to solve the Riemann hypothesis the periodicity of U ( σ , t ) with 1 / 2 ≤ σ ≤ 1 and t > 0 will be considered in a unified way.

Firstly, the 2-terms expression for ξ ( s ) in Equation (4) is awkward for use. In the following we adopt the derived form (see , pages 16-17):

ξ ( s ) = 1 2 − s ( 1 − s ) 2 ∫ 1 ∞     ψ ( x ) ( x s 2 − 1 + x − s + 1 2 ) d x = 4 ∫ 1 ∞ [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 cosh [ 1 2 ( s − 1 2 ) log x ] d x = 4 ∫ 1 ∞ [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 ⋅ 1 2 [ e 1 2 ( σ + i t − 1 2 ) log x + e − 1 2 ( σ + i t − 1 2 ) log x ] d x = 2 ∫ 1 ∞ [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 [ x 1 2 ( σ − 1 2 ) + x − 1 2 ( σ − 1 2 ) ] cos ( t 2 log x ) d x       + 2 i ∫ 1 ∞ [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 [ x 1 2 ( σ − 1 2 ) − x − 1 2 ( σ − 1 2 ) ] sin ( t 2 log x ) d x , (17)

here the superscript “ ' ” denotes the taking of conventional derivative. Hence,

U ( σ , t ) = 2 ∫ 1 ∞ [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 [ x 1 2 ( σ − 1 2 ) + x − 1 2 ( σ − 1 2 ) ] c o s ( t 2 l o g x ) d x ,

V ( σ , t ) = 2 ∫ 1 ∞ [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 [ x 1 2 ( σ − 1 2 ) − x − 1 2 ( σ − 1 2 ) ] s i n ( t 2 l o g x ) d x . (18)

In the following we mainly concern the variation of U ( σ , t ) .

To denote

f ( σ , x ) = 2 [ x 3 2 ψ ′ ( x ) ] ′ x − 1 4 [ x 1 2 ( σ − 1 2 ) + x − 1 2 ( σ − 1 2 ) ] ,

θ ( t ) = t 4 l o g t 2 π e + π 6 , (19)

with a substitution of

x = ( 2 π e t ) 1 / 2 e − π / 3 t y ,

it follows from the first formula in Equation (18) that

U ( σ , t ) = ∫ 1 ∞     f ( σ , x ) c o s ( t 2 l o g x ) d x = ∫ ( t 2 π e ) 1 / 2 e π / 3 t ∞ f [ σ , ( 2 π e t ) 1 / 2 e − π / 3 t y ] ( 2 π e t ) 1 / 2 e − π / 3 t     ⋅ cos { t 2 log [ ( 2 π e t ) 1 / 2 e − π / 3 t y ] } d y = ∫ ( t 2 π e ) 1 / 2 e π / 3 t ∞ f [ σ , ( 2 π e t ) 1 / 2 e − π / 3 t y ] ( 2 π e t ) 1 / 2 e − π / 3 t     ⋅ cos [ t 2 log y − θ ( t ) ] d y = 2 t [ I 1 ( σ , t ) cos θ ( t ) + I 2 ( σ , t ) sin θ ( t ) ] , (20)

where, under another substitution y = e 2 z / t ,

I 1 ( σ , t ) = t 2 ∫ ( t 2 π e ) 1 / 2 e π / 3 t ∞     f [ σ , ( 2 π e t ) 1 / 2 e − π / 3 t e 2 z / t ] ( 2 π e t ) 1 / 2 e − π / 3 t cos ( t 2 log y ) d y = ∫ θ ( t ) ∞     f [ σ , ( 2 π e t ) 1 / 2 e − π / 3 t e 2 z / t ] ( 2 π e t ) 1 / 2 e − π / 3 t e 2 z / t cos z d z = ∫ θ ( t ) ∞     F ( σ , x ( z , t ) ) cos z d z , (21)

I 2 ( σ , t ) = t 2 ∫ ( t 2 π e ) 1 / 2 e π / 3 t ∞     f [ σ , ( 2 π e t ) 1 / 2 e − π / 3 t y ] ( 2 π e t ) 1 / 2 e − π / 3 t s i n ( t 2 l o g y ) d y = ∫ θ ( t ) ∞     F ( σ , x ( z , t ) ) s i n z d z , (22)

in which F ( σ , x ) = f ( σ , x ) x with

x ( z , t ) = ( 2 π e t ) 1 / 2 e − π / 3 t y = ( 2 π e t ) 1 / 2 e − π / 3 t ⋅ e 2 z / t .

3.1. To Prove the Monotonicity of F ( σ , x ) about x

To recall that ψ ( x ) = ∑ n = 1 ∞     e − π n 2 x , with a denotation δ = ( σ − 1 / 2 ) / 2 , the definition of f ( σ , x ) in (19) leads to

F ( σ , x ) = 2 [ x 3 / 2 ψ ′ ( x ) ] ′ x − 1 / 4 ( x δ + x − δ ) ⋅ x = 2 [ 3 2 x 1 / 2 ψ ′ ( x ) + x 3 / 2 ψ ″ ( x ) ] x 3 / 4 ( x δ + x − δ ) = 2 [ 3 2 ψ ′ ( x ) + x ψ ″ ( x ) ] x 5 / 4 ( x δ + x − δ ) = 2 [ 3 2 ∑ n = 1 ∞ ( − π n 2 ) e − π n 2 x + x ∑ n = 1 ∞ ( − π n 2 ) 2 e − π n 2 x ] x 5 / 4 ( x δ + x − δ ) = 2 x 5 / 4 ( x δ + x − δ ) ∑ n = 1 ∞     a n ( a n x − 3 2 ) e − a n x (23)

with a n = π n 2 . Since π n 2 x − 3 / 2 > 0 for n ≥ 1 and x ≥ 1 , we have F ( σ , x ) > 0 . In addition, to recall that δ = ( σ − 1 / 2 ) / 2 , we have

∂ ∂ σ ( x δ + x − δ ) = l o g x 2 ( x δ − x − δ ) ≥ 0

for x ≥ 1 , and hence F ( σ , x ) increases along with the increasing of σ on [ 1 / 2 ,1 ] . By the way, f ( σ , x ) > 0 also holds and it follows from Equation (18) that

ξ ( 1 / 2 ) = U ( 1 / 2 ,0 ) = ∫ 1 ∞     f ( 1 / 2 , x ) d x > 0,

which accords with the documented result ξ ( 1 / 2 ) = − ζ ( 1 / 2 ) > 0 in  (see page 122).

To take the partial derivative of F ( σ , x ) with respect to x , it yields

F x ( σ , x ) = 2 ∑ n = 1 ∞     a n ∂ ∂ x [ ( x δ + x − δ ) ( a n x 9 / 4 − 3 2 x 5 / 4 ) e − a n x ] = 2 ∑ n = 1 ∞       a n δ ( x δ − 1 − x − δ − 1 ) ( a n x 9 / 4 − 3 2 x 5 / 4 ) e − a n x     + 2 ∑ n = 1 ∞     a n ( x δ + x − δ ) ( 9 4 a n x 5 / 4 − 15 8 x 1 / 4 ) e − a n x     − 2 ∑ n = 1 ∞     a n 2 ( x δ + x − δ ) ( a n x 9 / 4 − 3 2 x 5 / 4 ) e − a n x

= 2 δ ( x δ − x − δ ) x 1 / 4 ∑ n = 1 ∞     a n ( a n x − 3 2 ) e − a n x     + 2 ( x δ + x − δ ) x 1 / 4 ∑ n = 1 ∞     a n ( − a n 2 x 2 + 15 4 a n x − 15 8 ) e − a n x = 2 ( x δ + x − δ ) x 1 / 4 ∑ n = 1 ∞     a n [ − a n 2 x 2 + 15 4 a n x − 15 8 + δ x 2 δ − 1 x 2 δ + 1 ( a n x − 3 2 ) ] e − a n x

≤ 2 ( x δ + x − δ ) x 1 / 4 ∑ n = 1 ∞     a n [ − a n 2 x 2 + 15 4 a n x − 15 8     + 1 4 x − 1 1 + 1 ( a n x − 3 2 ) ] e − a n x = − 2 ( x δ + x − δ ) x 1 / 4 ∑ n = 1 ∞     a n [ a n ( a n − 1 8 ) x 2     − ( 29 8 a n − 3 16 ) x + 27 16 ] e − a n x (24)

due to the fact that for x ≥ 1 , n ≥ 1 and 0 ≤ δ ≤ 1 / 4 it holds a n x − 3 / 2 = π n 2 x − 3 / 2 > 0 and 1 ≤ x 2 δ ≤ x .

Let

h ( n , x ) = a n ( a n − 1 8 ) x 2 − ( 29 8 a n − 3 16 ) x + 27 16 ,

to solve the equation h ( n , x ) = 0 it yields two roots

x 1,2 = [ ( 29 / 8 ) a n − 3 / 16 ] ± Δ 2 a n ( a n − 1 / 8 ) ,

in which x 2 stands for the bigger one. Here

Δ = ( 29 8 a n − 3 16 ) 2 − 4 a n ( a n − 1 8 ) ⋅ 27 16 = 1 64 [ a n ( 409 a n − 33 ) + 9 4 ] > 0 ,

due to the fact that 409 a n − 33 = 409 π n 2 − 33 > 0 for n ≥ 1 .

For n = 1,2,3, ⋯ , the values of x 2 are about 1.0047,0.2464,0.1091, ⋯ . In view of x ≥ 1 , we have h ( n , x ) > 0 for n ≥ 2 . In case n = 1 , the inequality h ( 1, x ) > 0 also holds for x > x 2 ≈ 1.0047 . On [ 1, x 2 ] , h ( 1, x ) ≤ 0 and to ensure F x ( σ , x ) < 0 it suffices to consider the variation of the summation on this interval. In fact, on this interval,

∑ n = 1 ∞     a n [ a n ( a n − 1 8 ) x 2 − ( 29 8 a n − 3 16 ) x + 27 16 ] e − a n x > ∑ n = 1 3     a n [ a n ( a n − 1 8 ) x 2 − ( 29 8 a n − 3 16 ) x + 27 16 ] e − a n x

whose value increases from 8.1862 × 1017 to 0.0040. Hence, for all σ ∈ [ 1 / 2 ,1 ] , F x ( σ , x ) < 0 always holds and F ( σ , x ) is monotone decreasing about x on [ 1, ∞ ) .

3.2. To Prove I 1 2 ( σ , t ) + I 2 2 ( σ , t ) > 0 and Estimate Its Upper Bound

In order to combine the two terms I 1 ( σ , t ) c o s θ ( t ) + I 2 ( σ , t ) s i n θ ( t ) into a single periodic function, it requires I 1 2 ( σ , t ) + I 2 2 ( σ , t ) > 0 , that is I 1 2 ( σ , t ) + I 2 2 ( σ , t ) > 0 . To split z into z = z * + θ ( t ) , then

x ( z , t ) = ( 2 π e t ) 1 / 2 e − π / 3 t e 2 z / t = ( 2 π e t ) 1 / 2 e − π / 3 t e 2 ( z * + θ ( t ) ) / t = e 2 z * / t . (25)

It follows from Equation (21) and Equation (22) that

I 1 2 ( σ , t ) + I 2 2 ( σ , t ) = ∫ θ ( t ) ∞     F ( σ , x ( z , t ) ) cos z d z ⋅ ∫ θ ( t ) ∞     F ( σ , x ( w , t ) ) cos w d w       + ∫ θ ( t ) ∞     F ( σ , x ( z , t ) ) sin z d z ⋅ ∫ θ ( t ) ∞     F ( σ , x ( w , t ) ) sin w d w = ∫ θ ( t ) ∞   ∫ θ ( t ) ∞     F ( σ , x ( z , t ) ) F ( σ , x ( w , t ) ) cos ( z − w ) d w d z = ∫ 0 ∞   ∫ 0 ∞     F ( σ , x ( z * + θ ( t ) , t ) ) F ( σ , x ( w * + θ ( t ) , t ) ) cos ( z * − w * ) d w * d z *

= ∫ 0 ∞   ∫ 0 ∞ F ( σ , e 2 z * / t ) F ( σ , e 2 w * / t ) cos ( z * − w * ) d w * d z * = ∫ 0 ∞     F ( σ , e 2 z * / t ) cos ( z * ) d z * ⋅ ∫ 0 ∞     F ( σ , e 2 w * / t ) cos ( w * ) d w *     + ∫ 0 ∞     F ( σ , e 2 z * / t ) sin ( z * ) d z * ⋅ ∫ 0 ∞     F ( σ , e 2 w * / t ) sin ( w * ) d w * = [ ∫ 0 ∞     F ( σ , e 2 z / t ) cos ( z ) d z ] 2 + [ ∫ 0 ∞     F ( σ , e 2 z / t ) sin ( z ) d z ] 2 . (26)

For a fixed t ( > 0 ) , the variations of F ( σ , e 2 z / t ) , F ( σ , e 2 z / t ) c o s ( z ) and F ( σ , e 2 z / t ) s i n ( z ) for the case with σ = 1 / 2 are simulated in Figure 4.

Notice that F x ( σ , x ) < 0 for x ≥ 1 we have

∂ ∂ z F ( σ , x ( z , t ) ) ) = F x ( σ , x ) ∂ x ∂ z = F x ( σ , x ) 2 t e 2 z / t < 0

for the case with x = e 2 z / t . Hence, for fixed t > 0 , F ( σ , e 2 z / t ) decreases along with the increasing ofz (see Figure 4). In view of this,

∫ 0 ∞     F ( σ , e 2 z / t ) sin ( z ) d z = ∑ k = 0 ∞   ∫ 2 k π 2 ( k + 1 ) π     F ( σ , e 2 z / t ) sin ( z ) d z = ∑ k = 0 ∞   ∫ 0 2 π     F ( σ , e 2 ( z * + 2 k π ) / t ) sin ( z * ) d z * = ∑ k = 0 ∞   ∫ 0 π [ F ( σ , e 2 ( z * + 2 k π ) / t ) − F ( σ , e 2 ( z * + 2 k π + π ) / t ) ] sin ( z * ) d z * = ∑ k = 0 ∞ [ F ( σ , e 2 ( z ′ k + 2 k π ) / t ) − F ( σ , e 2 ( z ′ k + 2 k π + π ) / t ) ] ∫ 0 π sin ( z * ) d z * = 2 ∑ k = 0 ∞ [ F ( σ , e 2 ( z ′ k + 2 k π ) / t ) − F ( σ , e 2 ( z ′ k + 2 k π + π ) / t ) ] > 0 (27)

with 0 ≤ z ′ k ≤ π , here the mean-value theorem of integrals is used.

To combine Equation (26) and Equation (27) it leads to I 1 2 ( σ , t ) + I 2 2 ( σ , t ) > 0 . The assertion is proved. In addition, I 1 2 ( σ , t ) + I 2 2 ( σ , t ) is also upper bounded.

In fact, on the one hand, it follows from Equation (27) that

∫ 0 ∞     F ( σ , e 2 z / t ) sin ( z ) d z = 2 F ( σ , e 2 z ′ 0 / t ) − 2 ∑ k = 1 ∞ [ F ( σ , e 2 ( z ′ k + 2 k π − π ) / t ) − F ( σ , e 2 ( z ′ k + 2 k π ) / t ) ] < 2 F ( σ , e 2 z ′ 0 / t ) ≤ 2 F ( σ , 1 ) = 2 F 0 ( 1 ) (28)

with

F 0 ( 1 ) = 4 ∑ n = 1 ∞     a n ( a n − 3 2 ) e − a n , (29)

which does not rely on σ [see the expression in Equation (23)]. On the other hand,

∫ 0 ∞     F ( σ , e 2 z / t ) cos ( z ) d z = ∫ 0 π / 2     F ( σ , e 2 z / t ) cos ( z ) d z + ∑ k = 0 ∞   ∫ ( 4 k + 1 ) π / 2 ( 4 k + 5 ) π / 2     F ( σ , e 2 z / t ) cos ( z ) d z = F ( σ , e 2 z 0 / t ) − ∑ k = 0 ∞   ∫ 0 2 π F ( σ , e [ 2 z * + ( 4 k + 1 ) π ] / t ) sin ( z * ) d z * = F ( σ , e 2 z 0 / t ) − ∑ k = 0 ∞   ∫ 0 π [ F ( σ , e [ 2 z * + ( 4 k + 1 ) π ] / t )         − F ( σ , e [ 2 z * + ( 4 k + 3 ) π ] / t ) ] sin ( z * ) d z * = F ( σ , e 2 z 0 / t ) − 2 ∑ k = 0 ∞ [ F ( σ , e [ 2 z ′ k + ( 4 k + 1 ) π ] / t ) − F ( σ , e [ 2 z ′ k + ( 4 k + 3 ) π ] / t ) ] (30)

with 0 ≤ z 0 ≤ π / 2 and 0 ≤ z ′ k ≤ π . Its upper bound is F ( σ , e 2 z 0 / t ) ≤ F 0 ( 1 ) , and lower bound is F ( σ , e 2 z 0 / t ) − 2 F ( σ , e ( 2 z ′ 0 + π ) / t ) ≥ F ( σ , e π / t ) − 2 F ( σ , e π / t ) = − F ( σ , e π / t ) > − F 0 ( 1 ) .

Based on the above estimations, we have the final result:

0 < I 1 2 ( σ , t ) + I 2 2 ( σ , t ) < F 0 2 ( 1 ) + [ 2 F 0 ( 1 ) ] 2 = 5 F 0 ( 1 ) . (31)

3.3. Detailed Arguments on the Periodicity of U ( σ , t )

Since I 1 2 ( σ , t ) + I 2 2 ( σ , t ) , the real part of ξ ( σ + i t ) can be combined into a single periodic one:

U ( σ , t ) = 2 t [ I 1 ( σ , t ) cos θ ( t ) + I 2 ( σ , t ) sin θ ( t ) ] = 2 t A ( σ , t ) [ sin ϕ ( σ , t ) cos θ ( t ) + cos ϕ ( σ , t ) sin θ ( t ) ] = 2 t A ( σ , t ) sin [ ϕ ( σ , t ) + θ ( t ) ] , (32)

in which A ( σ , t ) = I 1 2 ( σ , t ) + I 2 2 ( σ , t ) , θ ( t ) = ( t / 4 ) l o g ( t / 2 π e ) + π / 6 and

s i n ϕ ( σ , t ) = I 1 ( σ , t ) A ( σ , t ) ,   c o s ϕ ( σ , t ) = I 2 ( σ , t ) A ( σ , t ) . (33)

Notice that A ( σ , t ) > 0 , their signs are determined by I 1 ( σ , t ) and I 2 ( σ , t ) , respectively. In the following we take the first one as the research object.

Firstly, we make certain the variation of θ ( t ) = ( t / 4 ) l o g ( t / 2 π e ) + π / 6 . As the function ( t / 4 ) l o g ( t / 2 π e ) concerned, in case 0 < t < 2 π e it is negative; in case t > 2 π e it is positive. To take the derivative of θ ( t ) , it leads to a unique extreme point at t = 2 π . Notice that the first zero point of U ( 1 / 2 , t ) = ξ ( 1 / 2 + i t ) is at about t = 14.1347251 , we only consider the case with t ∈ [ 2 π , ∞ ) on which θ ( t ) is increasing with value ≥-π/3.

There is a fact that, along with the increasing of θ ( t ) , for every given σ ∈ [ 1 / 2 ,1 ] , the function I 1 ( σ , t ) oscillates in a periodic manner about t. Explicitly, for the case with θ ( t ˜ m ) = 2 m π + π / 2 ( m = 0,1,2, ⋯ ), it follows from Equation (21), Equation (26) and Equation (27) that

I 1 ( σ , t ˜ m ) = ∫ θ ( t ˜ m ) ∞     F ( σ , x ( z , t ˜ m ) ) cos z d z = ∫ 0 ∞     F ( σ , e 2 z * / t ˜ m ) cos [ z * + θ ( t ˜ m ) ] d z * = ∫ 0 ∞     F ( σ , e 2 z * / t ˜ m ) cos ( z * + 2 m π + π / 2 ) d z * = − ∫ 0 ∞ ( σ , e 2 z * / t ˜ m ) sin z * d z * < 0. (34)

Similarly, we have I 1 ( σ , t ^ m ) > 0 for the case with θ ( t ^ m ) = θ ( t ˜ m ) + π = 2 m π + 3 π / 2 ( m = 0,1,2, ⋯ ).

The above results accord with the fact that s i n ϕ ( σ , t ) shifts its sign when ϕ ( σ , t ) is added by π. This consistency implies that θ ( t ) contributes to the variation of s i n ϕ ( σ , t ) , and is possibly the main part of ϕ ( σ , t ) . In the following we verify this by considering the detailed variation of I 1 ( σ , t ) on the interval [ t ˜ m , t ^ m ] (together with another half period [ t ^ m , t ˜ m + 1 ] ).

For any given t 0 ( ≥ 2 π ) , in case t > t 0 we have

I 1 ( σ , t ) − I 1 ( σ , t 0 ) = ∫ θ ( t ) ∞     F ( σ , x ( z , t ) ) cos z d z − ∫ θ ( t 0 ) ∞     F ( σ , x ( z , t 0 ) ) cos z d z = − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) cos z d z + ∫ θ ( t ) ∞ [ F ( σ , x ( z , t ) ) − F ( σ , x ( z , t 0 ) ) ] cos z d z . (35)

Relative to θ ( t 0 ) = 2 m π + π / 2 (for a given m ), to verify that θ ( t ) is the main part of ϕ ( σ , t ) , it suffices to show that the monotone increasing term − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) cos z d z dominates the relative change of I 1 ( σ , t ) with respect to θ ( t ) ∈ [ 2 m π + π / 2 ,2 m π + 3 π / 2 ] [that is, for t ∈ ( t ˜ m , t ^ m ) ]. In other word, the contribution from the additional term is relatively small. In the following we estimate them in detail.

For the first term, in view of the monotone decreasing of F ( σ , x ( z , t ) ) about z, with a denotation Δ θ = θ ( t ) − θ ( t 0 ) and a substitution z = z * + θ ( t 0 ) , we have

− ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) cos z d z = − ∫ 0 Δ θ     F ( σ , e 2 z * / t 0 ) cos ( z * + 2 m π + π / 2 ) d z * = ∫ 0 Δ θ     F ( σ , e 2 z * / t 0 ) sin z * d z * . (36)

3.3.1. To Estimate the Additional Term

The second term is an additional one which is very complex. We estimate it step by step. In the neighbourhood of t 0 , to make the Taylor expansion of F ( σ , x ( z , t ) ) it leads to

F ( σ , x ( z , t ) ) = F ( σ , x ( z , t 0 ) ) + F t ( σ , x ( z , t 0 ) ) Δ t + ⋯ ,

where Δ t = t − t 0 , the subscript “t” means the taking of partial derivative with respect to t.

Recall that x ( z , t ) = ( 2 π e / t ) 1 / 2 e − π / 3 t e 2 z / t , to replace z by z * + θ ( t ) = z * + ( t / 4 ) l o g ( t / 2 π e ) + π / 6 , it leads to

∂ x ∂ t | t = t 0 = − 1 t [ 1 2 + 2 ( z − π / 6 ) t ] ( 2 π e t ) 1 / 2 e − π / 3 t e 2 z / t | t = t 0 = − 1 t [ 1 2 + 2 ( z * + θ ( t ) − π / 6 ) t ] ( 2 π e t ) 1 / 2 e − π / 3 t e 2 ( z * + θ ( t ) ) / t | t = t 0 = − 1 t 0 ( 2 z * t 0 + 1 2 log t 0 2 π ) e 2 z * / t 0 . (37)

Furthermore, notice that x ( z , t 0 ) = x ( z * + θ ( t 0 ) , t 0 ) = e 2 z * / t 0 , to recall the expression of F x ( σ , x ) in Equation (24), we get

F t ( σ , x ( z , t 0 ) ) = F x ( σ , x ) | z = z * + θ ( t 0 ) ⋅ ∂ x ∂ t | t = t 0 = ∑ n = 1 ∞     a n [ − a n 2 x 2 + 15 4 a n x − 15 8 + δ x 2 δ − 1 x 2 δ + 1 ( a n x − 3 2 ) ] e − a n x     ⋅ 2 ( x δ + x − δ ) x 1 / 4 ⋅ [ − 1 t 0 ( 2 z * t 0 + 1 2 log t 0 2 π ) x ]

= 2 t 0 ( log x + 1 2 log t 0 2 π ) ( x δ + x − δ ) x 5 / 4     ⋅ ∑ n = 1 ∞     a n [ a n 2 x 2 − 15 4 a n x + 15 8 − δ x 2 δ − 1 x 2 δ + 1 ( a n x − 3 2 ) ] e − a n x = : 1 t 0 G 1 ( σ , x ) + 1 2 t 0 log t 0 2 π G 2 ( σ , x ) (38)

with x ( z * , t 0 ) = e 2 z * / t 0 , a n = π n 2 and δ = ( σ − 1 / 2 ) / 2 .

It follows from the proofs in Section 3.1 that G 1 ( σ , e 2 z * / t 0 ) > 0 and G 2 ( σ , e 2 z * / t 0 ) > 0 for z * ≥ 0 , t 0 ≥ 2 π and 1 / 2 ≤ σ ≤ 1 . With respect to the single variable w = 2 z * / t 0 , each of them has a unique extreme point (maximum one), say w 1 and w 2 respectively, which can be estimated by their dominate terms with n = 1 . Due to the existence of strong decaying factor e − a n e w , the effect from σ is very small, see Figure 5. In the following we mainly consider the case with σ = 1 / 2 .

The approximate values for the maximum extreme points of G 1 ( 1 / 2 , e w ) and

G 2 ( 1 / 2 , e w ) are at about ( 0.6484,0.6957 ) and ( 0.4746,1.2304 ) , respectively. Explicitly, w 1 ≈ 0.6484 and w 2 ≈ 0.4746 satisfy separately the equations ∂ G 1 ( 1 / 2 , e w ) / ∂ w = 0 and ∂ G 2 ( 1 / 2 , e w ) / ∂ w = 0 (for each of them only the first term is concerned), which read:

( y 1 3 − 7 y 1 2 + 165 16 y 1 − 75 32 ) l o g ( y 1 π ) = y 1 2 − 15 4 y 1 + 15 8 ,

y 2 3 − 7 y 2 2 + 165 16 y 2 − 75 32 = 0

with y 1 = π e w 1 and y 2 = π e w 2 . By the way, the maximum extreme point of F ( 1 / 2 , e w ) is at about ( 0,0.8934 ) (to ensure the monotonicity, the first 3 terms are considered), which will be used in the related estimation.

When σ varies in [ 1 / 2 ,1 ] , the locations of the maximum extreme points of G 1 ( σ , e w ) and G 2 ( σ , e w ) have small shifts. Explicitly, when σ = 1 / 2 the location of G 1 ( σ , e w ) is at about ( 0.6484,0.6957 ) , and along with the increasing of σ its value increases a little. At the right end σ = 1 , the location is shifted to about ( 0.6538,0.6967 ) . Similarly, that of G 2 ( σ , e w ) is shifted from about ( 0.4746,1.2304 ) to about ( 0.4820,1.2237 ) . These indicate that the effect of σ to the maximum values of G 1 and G 2 can be neglected, yet the effects of σ to the horizontal shifting can not be neglected, since z * = w 1 t 0 / 2 and z * = w 2 t 0 / 2 may magnify the small differences along with the increasing of t 0 and result in considerable impacts to the corresponding integrals. By the way, we note that w 1 and w 2 depend on the variable σ , and should be in the forms w 1 ( σ ) and w 2 ( σ ) , to neglect the variable is just for simplicity.

Based on the above analysis, we estimate the second term of Equation (35) in the following:

∫ θ ( t ) ∞ [ F ( σ , x ( z , t ) ) − F ( σ , x ( z , t 0 ) ) ] cos z d z ≈ ∫ θ ( t ) ∞     F t ( σ , x ( z , t 0 ) ) Δ t cos z d z = Δ t ∫ Δ θ ∞     F t ( σ , x ( z * + θ ( t 0 ) , t 0 ) ) cos ( z * + θ ( t 0 ) ) d z * = Δ t ∫ Δ θ ∞ [ 1 t 0 G 1 ( σ , e 2 z * / t 0 ) + 1 2 t 0 log t 0 2 π G 2 ( σ , e 2 z * / t 0 ) ] cos ( z * + 2 m π + π / 2 ) d z * = − Δ t t 0 ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * − Δ t 2 t 0 log t 0 2 π ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) sin z * d z * . (39)

In the following we give some illustrations for the case σ = 1 / 2 . In view of Figure 5, G 1 ( 1 / 2 , e 2 z * / t 0 ) has maximum value at 2 z * / t 0 = w 1 , that is, at about z * = t 0 w 1 / 2 ≈ 0.3242 t 0 . Similarly, G 2 ( 1 / 2 , e 2 z * / t 0 ) has maximum value at about z * = t 0 w 2 / 2 ≈ 0.2373 t 0 . For example, for the cases with m = 0 , 1 and 2, to solve

θ ( t 0 ) = t 0 4 l o g ( t 0 2 π e ) + π 6 = 2 m π + π / 2 ,

it yields t 0 = 20.8747,37.4043 and 50.3595, respectively. The corresponding maximum extreme points of G 1 ( 1 / 2 , e 2 z * / t 0 ) are shifted to ( 6.7676,0.6957 ) , ( 12.1265,0.6957 ) and ( 16.3265,0.6957 ) , respectively. Similarly, those for G 2 ( 1 / 2 , e 2 z * / t 0 ) are shifted to ( 4.9536,1.2304 ) , ( 8.8760,1.2304 ) and ( 11.9503 , 1.2304 ) , respectively (see Figure 6). What should be mentioned is that, in the neighborhoods of these two maximum extreme points, G 1 ( 1 / 2 , e 2 z * / t 0 ) s i n z * and G 2 ( 1 / 2 , e 2 z * / t 0 ) s i n z * may possess their peak absolute values separately, which dominate the integrals. Yet, since t 0 w 1 / 2 ≈ 0.3242 t 0 and t 0 w 2 / 2 ≈ 0.2373 t 0 , the changes of G 1 ( 1 / 2 , e 2 z * / t 0 ) and G 2 ( 1 / 2 , e 2 z * / t 0 ) are not synchronous with that of θ ( t ) . This adds complexity to the additional term.

Notice that w = w 1 is the unique extreme point of G 1 ( σ , e w ) at which it possesses the maximum value, we get an increasing interval [ 0 , t 0 w 1 / 2 ] and a decreasing interval [ t 0 w 1 / 2 , ∞ ) for G 1 ( σ , e 2 z * / t 0 ) with respect to the variable z * . For m ≥ 0 , there must be an integer k 0 ≥ 1 such that z * = t 0 w 1 / 2 ∈ [ k 0 π , ( k 0 + 1 ) π ) (see Figure 6). Furthermore, in case k 0 = 2 j 0 with j 0 ≥ 1 , we have

∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * = − ∫ 0 Δ θ     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * + ∑ k = 0 ∞ ( − 1 ) k ∫ 0 π     G 1 ( σ , e 2 ( z + k π ) / t 0 ) sin z d z = − ∫ 0 Δ θ     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * + ∫ 0 π     G 1 ( σ , e 2 ( z + 2 j 0 π ) / t 0 ) sin z d z     − ∑ j = 0 j 0 − 1   ∫ 0 π [ G 1 ( σ , e 2 ( z + ( 2 j + 1 ) π ) / t 0 ) − G 1 ( σ , e 2 ( z + 2 j π ) / t 0 ) ] sin z d z     − ∑ j = j 0 + 1 ∞   ∫ 0 π [ G 1 ( σ , e 2 ( z + ( 2 j − 1 ) π ) / t 0 ) − G 1 ( σ , e 2 ( z + 2 j π ) / t 0 ) ] sin z d z < ∫ 0 π     G 1 ( σ , e 2 ( z + 2 j 0 π ) / t 0 ) sin z d z < 2 G 1 ( σ , e w 1 ) ≤ 2 G 1 ( 1 , e w 1 ) . (40)

Here the approximation 0.6957 ≈ G 1 ( 1 / 2 , e w 1 ) ≤ G 1 ( σ , e w 1 ) ≤ G 1 ( 1 , e w 1 ) ≈ 0.6967 is used.

Similarly, in case k 0 = 2 j 0 − 1 with j 0 > 1 , we have

∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * = ∫ Δ θ π     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * + ∑ k = 1 ∞ ( − 1 ) k ∫ 0 π     G 1 ( σ , e 2 ( z + k π ) / t 0 ) sin z d z = ∫ Δ θ π     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * − ∫ 0 π     G 1 ( σ , e 2 ( z + ( 2 j 0 − 1 ) π ) / t 0 ) sin z d z     + ∑ j = 1 j 0 − 1   ∫ 0 π [ G 1 ( σ , e 2 ( z + 2 j π ) / t 0 ) − G 1 ( σ , e 2 ( z + ( 2 j − 1 ) π ) / t 0 ) ] sin z d z     + ∑ j = j 0 ∞   ∫ 0 π [ G 1 ( σ , e 2 ( z + 2 j π ) / t 0 ) − G 1 ( σ , e 2 ( z + ( 2 j + 1 ) π ) / t 0 ) ] sin z d z > − ∫ 0 π     G 1 ( σ , e 2 ( z + ( 2 j 0 − 1 ) π ) / t 0 ) sin z d z > − 2 G 1 ( σ , e w 1 ) ≥ − 2 G 1 ( 1 , e w 1 ) . (41)

In a word, with respect to θ ( t 0 ) = 2 m π + π / 2 , for all m ≥ 0 it holds the estimation bellow:

| ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) s i n z * d z * | < 2 G 1 ( 1 , e w 1 ) ≈ 2 × 0.6967 = 1.3934. (42)

When G 1 ( σ , e 2 z * / t 0 ) is substituted by G 2 ( σ , e 2 z * / t 0 ) , with the same approach we get a similar estimation:

| ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) s i n z * d z * | < 2 G 2 ( 1 / 2 , e w 2 ) ≈ 2 × 1.2304 = 2.4608, (43)

since for this case 1.2237 ≈ G 2 ( 1, e w 2 ) ≤ G 2 ( σ , e w 2 ) ≤ G 2 ( 1 / 2 , e w 2 ) ≈ 1.2304 .

3.3.2. To Estimate the Change of I 1 ( σ , t ) with Respect to θ ( t )

With the help of the above results, in the following we estimate the relative change I 1 ( σ , t ) − I 1 ( σ , t 0 ) with respect to θ ( t ) ∈ ( 2 m π + π / 2 ,2 m π + 3 π / 2 ) ( m ≥ 0 ).

Relative to θ ( t 0 ) = 2 m π + π / 2 , there is a Taylor expansion for θ ( t ) :

Δ θ = θ ( t ) − θ ( t 0 ) = θ ′ ( t ) Δ t + θ ″ ( t ) ( Δ t ) 2 2 ! + ⋯ = 1 4 l o g t 0 2 π Δ t + 1 8 t 0 ( Δ t ) 2 + ⋯ . (44)

Hence, Δ t ≈ 4 Δ θ / l o g ( t 0 / 2 π ) . It follows from Equations (35), (36), (39), (42) and (43) that

I 1 ( σ , t ) − I 1 ( σ , t 0 ) = − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) cos z d z + ∫ θ ( t ) ∞ [ F ( σ , x ( z , t ) ) − F ( σ , x ( z , t 0 ) ) ] cos z d z ≈ ∫ 0 Δ θ     F ( σ , e 2 z * / t 0 ) sin z * d z * − Δ t t 0 ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z *     − Δ t 2 t 0 log t 0 2 π ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) sin z * d z *

≈ ∫ 0 Δ θ     F ( σ , e 2 z * / t 0 ) sin z * d z * − 4 Δ θ t 0 log ( t 0 / 2 π ) ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z *     − 2 Δ θ t 0 ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) sin z * d z * = : R ( σ , Δ θ ) . (45)

Furthermore, in view of 0 < Δ θ < π , to take the partial derivative of R ( σ , Δ θ ) with respect to Δ θ , it reads

∂ R ∂ Δ θ = F ( σ , e 2 Δ θ / t 0 ) sin Δ θ − 4 t 0 log ( t 0 / 2 π ) ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z *                     − 2 t 0 ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) sin z * d z *                     + 2 Δ θ t 0 [ 2 log ( t 0 / 2 π ) G 1 ( σ , e 2 Δ θ / t 0 ) + G 2 ( σ , e 2 Δ θ / t 0 ) ] sin Δ θ

> F ( σ , e 2 π / t 0 ) sin Δ θ − 4 t 0 log ( t 0 / 2 π ) | ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z * |     − 2 t 0 | ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) sin z * d z * | > F ( 1 / 2 , e 2 π / t 0 ) sin Δ θ − 4 t 0 [ 2 G 1 ( 1 , e w 1 ) log ( t 0 / 2 π ) + G 2 ( 1 / 2 , e w 2 ) ] ≈ F ( 1 / 2 , e 2 π / t 0 ) Δ θ − 4 t 0 [ 2 G 1 ( 1 , e w 1 ) log ( t 0 / 2 π ) + G 2 ( 1 / 2 , e w 2 ) ] > F ( 1 / 2 , e 2 π / t 0 ) Δ θ − 8 G 2 ( 1 / 2 , e w 2 ) t 0 > 0 , (46)

provided that

Δ θ > 8 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) .

Here the inequality 2 G 1 ( 1, e w 1 ) / log ( t 0 / 2 π ) < G 2 ( 1 / 2 , e w 2 ) is also used. In fact, to solve θ ( t 0 ) = π / 2 it yields t 0 ≈ 20.8747 . To recall that G 1 ( 1, e w 1 ) ≈ 0.6967 and G 2 ( 1 / 2 , e w 2 ) ≈ 1.2304 , we get

2 G 1 ( 1, e w 1 ) l o g ( t 0 / 2 π ) G 2 ( 1 / 2 , e w 2 ) ≈ 0.9432 < 1.

Along with the increasing of m , that is t 0 , the ratio gets smaller and smaller, and the inequality always holds.

For the case with m = 0 and t 0 ≈ 20.8747 , the approximate value of F ( 1 / 2 , e 2 π / t 0 ) is about 0.7205. At this time,

8 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) ≈ 8 × 1.2304 20.8747 × 0.7205 ≈ 0.2083 π .

This indicates that, for the case with θ ( t ) ∈ ( π / 2 , 3 π / 2 ) , no matter the additional term

∫ θ ( t ) ∞ [ F ( σ , x ( z , t ) ) − F ( σ , x ( z , t 0 ) ) ] c o s z d z

is positive or not, the increasing term − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) c o s z d z dominates the relative change of I 1 ( σ , t ) for all σ ∈ [ 1 / 2 ,1 ] , at least, with respect to Δ θ on the interval ( 0.2083 π , π ) . In another word, I 1 ( σ , t ) increases along with the increasing of θ ( t ) on ( 0.7083 π , 3 π / 2 ) (ensured by ∂ R / ∂ Δ θ > 0 ). For the other cases with m ≥ 1 , the interval may become wider than this. To summarize the results, we get a theorem below:

Theorem 1. For the case with

θ ( t ) = t 4 l o g ( t 2 π e ) + π 6 ∈ [ 2 m π + π / 2 ,2 m π + 3 π / 2 )

( m is a non-negative integer), let θ ( t 0 ) = 2 m π + π / 2 ,then for all σ ∈ [ 1 / 2 ,1 ] , I 1 ( σ , t ) increases along with the increasing of θ ( t ) , at least on the interval defined by

2 m π + π / 2 + 8 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) < θ ( t ) < 2 m π + 3 π / 2 . (47)

In the following, we begin with assuming θ ( t 0 ) = 2 m π − π / 2 and consider the case θ ( t ) ∈ ( 2 m π − π / 2 ,2 m π + π / 2 ) with m ≥ 1 . To repeat the previous deduction process it results in

I 1 ( σ , t ) − I 1 ( σ , t 0 ) = − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) cos z d z + ∫ θ ( t ) ∞ [ F ( σ , x ( z , t ) ) − F ( σ , x ( z , t 0 ) ) ] cos z d z ≈ − ∫ 0 Δ θ     F ( σ , e 2 z * / t 0 ) sin z * d z * + 4 Δ θ t 0 log ( t 0 / 2 π ) ∫ Δ θ ∞     G 1 ( σ , e 2 z * / t 0 ) sin z * d z *       + 2 Δ θ t 0 ∫ Δ θ ∞     G 2 ( σ , e 2 z * / t 0 ) sin z * d z * . (48)

Furthermore, to denote it by R ( σ , Δ θ ) , then we get ∂ R / ∂ Δ θ < 0 on condition that

Δ θ > 8 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) .

Hence, there is a similar result bellow:

Theorem 2. For the case with

θ ( t ) = t 4 l o g ( t 2 π e ) + π 6 ∈ [ 2 m π − π / 2 ,2 m π + π / 2 )

( m is a positive integer),let θ ( t 0 ) = 2 m π − π / 2 ,then for all σ ∈ [ 1 / 2 ,1 ] , I 1 ( σ , t ) decreases along with the decreasing of θ ( t ) ,at least on the interval defined by

2 m π − π / 2 + 8 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) < θ ( t ) < 2 m π + π / 2 .

These two theorems have clarified the periodic behaviors of I 1 ( σ , t ) . It follows from Theorem 1 that, due to the existence of decay factor 1 / t 0 in Equation (47) ( 1 / F ( 1 / 2 , e 2 π / t 0 ) also decays), the dominating interval of the increasing term − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) c o s z d z becomes wider and wider along with the increasing of t 0 . In another word, the possible interval for I 1 ( σ , t ) < I 1 ( σ , t 0 ) becomes narrower and narrower. This implies that the minimum extreme point of I 1 ( σ , t ) get closer and closer to the left end of the interval [ 2 m π + π / 2 ,2 m π + 3 π / 2 ) with respect to θ ( t ) as m → ∞ ( m ≥ 0 ). Certainly, this only makes sense when the additional term is negative. In case it is positive, the thing becomes simple. The minimum value of I 1 ( σ , t ) is attained at the left boundary t = t 0 defined by θ ( t 0 ) = 2 m π + π / 2 . Furthermore, from Theorem 2 we see, in the left neighborhood of this boundary it is monotone decreasing, and the minimum extreme point is exactly at t = t 0 which accords with θ ( t 0 ) = 2 m π + π / 2 . Inversely, the maximum extreme point of I 1 ( σ , t ) get closer and closer to the left end of the interval [ 2 m π − π / 2 ,2 m π + π / 2 ) with respect to θ ( t ) as m → ∞ . It only makes sense when the additional term is positive, and in case this term is negative the maximum extreme point is exactly at t = t 0 with θ ( t 0 ) = 2 m π − π / 2 ( m ≥ 1 ).

Based on the above analysis we get a corollary below:

Corollary 1. For the case with

θ ( t ) = t 4 l o g ( t 2 π e ) + π 6 ∈ [ 2 m π ,2 ( m + 1 ) π ]

( m is a non-negative integer), for every given σ ∈ [ 1 / 2 ,1 ] # there is a minimum extreme point and a maximum extreme point of I 1 ( σ , t ) ,which lie separately in the intervals defined by

2 m π + π / 2 ≤ θ ( t ) ≤ 2 m π + π / 2 + 8 G 2 ( 1 / 2 , e w 2 ) t ˜ m F ( 1 / 2 , e 2 π / t ˜ m ) ,

2 m π + 3 π / 2 ≤ θ ( t ) ≤ 2 m π + 3 π / 2 + 8 G 2 ( 1 / 2 , e w 2 ) t ^ m F ( 1 / 2 , e 2 π / t ^ m )

with θ ( t ˜ m ) = 2 m π + π / 2 and θ ( t ^ m ) = 2 m π + 3 π / 2 . As m → ∞ ,these two extreme points in the moving intervals tend separately to their left boundaries t ˜ m and t ^ m .

3.3.3. To Estimate the Change of I 1 ( σ , t ) / A ( σ , t ) with Respect to θ ( t )

According to the definition,

s i n ϕ ( σ , t ) = I 1 ( σ , t ) A ( σ , t ) .

Hence, besides the estimation on the change of I 1 ( σ , t ) , it also requires an estimation on that of A ( σ , t ) .

It follows from Equation (26) that A ( σ , t ) = H 2 ( σ , t ) + S 2 ( σ , t ) with

H ( σ , t ) = ∫ 0 ∞     F ( σ , e 2 z / t ) c o s ( z ) d z ,

S ( σ , t ) = ∫ 0 ∞     F ( σ , e 2 z / t ) s i n ( z ) d z .

To take the partial derivative of A ( σ , t ) about t, it leads to

A t ( σ , t ) = H ( σ , t ) H t ( σ , t ) + S ( σ , t ) S t ( σ , t ) H 2 ( σ , t ) + S 2 ( σ , t ) . (49)

Let x ( z , t ) = e 2 z / t , then ∂ x / ∂ t = − ( 2 z / t 2 ) e 2 z / t . It follows from Equation (38) that

H t ( σ , t ) = ∫ 0 ∞     F t ( σ , x ( z , t ) ) cos ( z ) d z , = ∫ 0 ∞     F x ( σ , x ) ∂ x ∂ t cos ( z ) d z , = 1 t ∫ 0 ∞     G 1 ( σ , e 2 z / t ) cos ( z ) d z .

S t ( σ , t ) = 1 t ∫ 0 ∞     G 1 ( σ , e 2 z / t ) s i n ( z ) d z .

To repeat the deduction processes in Equation (40) and Equation (41) we get | H t ( σ , t ) | ≤ 2 G 1 ( 1, e w 1 ) / t and | S t ( σ , t ) | ≤ 2 G 1 ( 1, e w 1 ) / t .

For the reference point defined by θ ( t 0 ) = 2 m π + π / 2 with m ≥ 0 , it holds Δ t ≈ 4 Δ θ / l o g ( t 0 / 2 π ) [see Equation (44)], and hence

| Δ A | = | A ( σ , t ) − A ( σ , t 0 ) | ≈ | A t ( σ , t 0 ) | Δ t ≈ | H ( σ , t 0 ) H t ( σ , t 0 ) + S ( σ , t 0 ) S t ( σ , t 0 ) | H 2 ( σ , t 0 ) + S 2 ( σ , t 0 ) ⋅ 4 Δ θ log ( t 0 / 2 π ) ≤ | H ( σ , t 0 ) | ⋅ | H t ( σ , t 0 ) | + | S ( σ , t 0 ) | ⋅ | S t ( σ , t 0 ) | H 2 ( σ , t 0 ) + S 2 ( σ , t 0 ) ⋅ 4 Δ θ log ( t 0 / 2 π ) < 2 G 1 ( 1, e w 1 ) t 0 | H ( σ , t 0 ) | + | S ( σ , t 0 ) | H 2 ( σ , t 0 ) + S 2 ( σ , t 0 ) ⋅ 4 Δ θ log ( t 0 / 2 π ) ≤ 2 G 1 ( 1, e w 1 ) t 0 2 ⋅ 4 Δ θ log ( t 0 / 2 π ) ≈ 8 2 G 1 ( 1, e w 1 ) t 0 log ( t 0 / 2 π ) Δ θ . (50)

Due to the existence of the decay factor 1 / t 0 l o g ( t 0 / 2 π ) , this estimation makes sense. For example, for the cases with m = 1,10 and 100, the approximate values are about 0.1572 Δ θ , 0.0213 Δ θ and 0.0025 Δ θ , respectively.

For an arbitrary given reference point t 0 ∈ ( 2 π , ∞ ) and σ ∈ [ 1 / 2 ,1 ] we have

sin ϕ ( σ , t ) = I 1 ( σ , t ) A ( σ , t ) = I 1 ( σ , t ) A ( σ , t 0 ) ⋅ A ( σ , t 0 ) A ( σ , t ) . (51)

Since A ( σ , t 0 ) > 0 , the change characteristics of I 1 ( σ , t ) are maintained by I 1 ( σ , t ) / A ( σ , t 0 ) . The positive factor A ( σ , t 0 ) / A ( σ , t ) does not influence the sign change of sin ϕ ( σ , t ) . As the function I 1 ( σ , t ) / A ( σ , t ) concerned, for θ ( t 0 ) = 2 m π + π / 2 , I 1 ( σ , t 0 ) / A ( σ , t 0 ) < 0 , and for θ ( t 0 ) = 2 m π + 3 π / 2 , I 1 ( σ , t 0 ) / A ( σ , t 0 ) > 0 . That is to say, the oscillation characteristic of I 1 ( σ , t ) is maintained by I 1 ( σ , t ) / A ( σ , t ) . But the variation of the amplitude does make a difference, and a further estimation is needed.

Our task is to verify that the change of I 1 ( σ , t ) / A ( σ , t ) is mainly associated with that of θ ( t ) , that is, the term − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) c o s z d z also dominate. Since this term increases on [ 2 m π + π / 2 ,2 m π + 3 π / 2 ] ( m ≥ 0 ) and decreases on [ 2 m π − π / 2 ,2 m π + π / 2 ] ( m ≥ 1 ) with respect to θ ( t ) , in the following we still consider it in a separate way.

For a fixed σ , relative to t = t 0 defined by θ ( t 0 ) = 2 m π + π / 2 , to denote A ( σ , t ) = A ( σ , t 0 ) + Δ A and I 1 ( σ , t ) = I 1 ( σ , t 0 ) + Δ I 1 , then

s i n ϕ ( σ , t ) − s i n ϕ ( σ , t 0 ) = I 1 ( σ , t ) A ( σ , t ) − I 1 ( σ , t 0 ) A ( σ , t 0 ) = A ( σ , t 0 ) Δ I 1 − I 1 ( σ , t 0 ) Δ A [ A ( σ , t 0 ) + Δ A ] A ( σ , t 0 ) . (52)

Since the unchangeable factor 1 / A ( σ , t 0 ) does not impact the relative change of this increment, in the following arguments we omit it for simplicity.

Since the inequality A ( σ , t ) ≥ | I 1 ( σ , t ) | holds for all t > 2 π , there always be A ( σ , t 0 ) + Δ A ≥ | I 1 ( σ , t 0 ) + Δ I 1 | , and it follows from Equation (45), Equation (46) and Equation (50) that

∂ ∂ Δ θ [ A ( σ , t 0 ) Δ I 1 − I 1 ( σ , t 0 ) Δ A A ( σ , t 0 ) + Δ A ] = A A + Δ A ⋅ ∂ Δ I 1 ∂ Δ θ − I 1 A + Δ A ⋅ ∂ Δ A ∂ Δ θ − A Δ I 1 − I 1 Δ A [ A + Δ A ] 2 ⋅ ∂ Δ A ∂ Δ θ = A A + Δ A [ ∂ Δ I 1 ∂ Δ θ − I 1 + Δ I 1 A + Δ A ⋅ ∂ Δ A ∂ Δ θ ] ≥ A A + Δ A [ ∂ Δ I 1 ∂ Δ θ − | ∂ Δ A ∂ Δ θ | ]

> A A + Δ A { F ( 1 / 2 , e 2 π / t 0 ) Δ θ − 4 t 0 [ 2 G 1 ( 1, e w 1 ) log ( t 0 / 2 π ) + G 2 ( 1 / 2 , e w 2 ) ]         − 8 2 G 1 ( 1, e w 1 ) t 0 log ( t 0 / 2 π ) }

= A A + Δ A { F ( 1 / 2 , e 2 π / t 0 ) Δ θ − 4 t 0 [ 2 ( 2 + 1 ) G 1 ( 1, e w 1 ) log ( t 0 / 2 π ) + G 2 ( 1 / 2 , e w 2 ) ] } > A A + Δ A [ F ( 1 / 2 , e 2 π / t 0 ) Δ θ − 14 t 0 G 2 ( 1 / 2 , e w 2 ) ] > 0 (53)

provided that

Δ θ > 14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) . (54)

Here the inequality 2 ( 2 + 1 ) G 1 ( 1, e w 1 ) / log ( t 0 / 2 π ) < ( 5 / 2 ) G 2 ( 1 / 2 , e w 2 ) is used. In fact, for the case with m = 0 , to recall that G 1 ( 1, e w 1 ) ≈ 0.6967 , G 2 ( 1 / 2 , e w 2 ) ≈ 1.2304 and t 0 ≈ 20.8747 , we get

2 ( 2 + 1 ) G 1 ( 1, e w 1 ) l o g ( t 0 / 2 π ) G 2 ( 1 / 2 , e w 2 ) ≈ 2 ( 2 + 1 ) × 0.6967 l o g ( 20.8747 / 2 π ) × 1.2304 ≈ 2.2771 < 5 / 2 .

Along with the increasing of m , t 0 increases accordingly, and this proportion will become smaller and smaller. So, the inequality always holds true.

For the case with m = 0 and t 0 ≈ 20.8747 , the approximate value of F ( 1 / 2 , e 2 π / t 0 ) is about 0.7205. At this time,

14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) ≈ 14 × 1.2304 20.8747 × 0.7205 ≈ 0.3646 π < 2 5 π .

This indicates that, for the case with θ ( t ) ∈ ( π / 2 , 3 π / 2 ) , no matter the additional term

∫ θ ( t ) ∞ [ F ( σ , x ( z , t ) ) − F ( σ , x ( z , t 0 ) ) ] c o s z d z

of I 1 ( σ , t ) is positive or not, the increasing term − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) c o s z d z

of I 1 ( σ , t ) dominates the relative change of I 1 ( σ , t ) / A ( σ , t ) for all σ ∈ [ 1 / 2 ,1 ] , at least, with respect to Δ θ = θ ( t ) − θ ( t 0 ) on the interval ( 0.3646 π , π ) . In an- other word, s i n ϕ ( σ , t ) increases along with the increasing of θ ( t ) on ( 0.8646 π , 3 π / 2 ) . For the other cases with m ≥ 1 , the interval becomes wider than this.

In summery, for the case with θ ( t ) ∈ [ 2 m π + π / 2 ,2 m π + 3 π / 2 ) ( m ≥ 0 ), the results for I 1 ( σ , t ) / A ( σ , t ) are similar to that of I 1 ( σ , t ) in Theorem 1. The mimic theorem is as follows:

Theorem 3. For the case with

θ ( t ) = t 4 l o g ( t 2 π e ) + π 6 ∈ [ 2 m π + π / 2 ,2 m π + 3 π / 2 )

( m is a non-negative integer), let θ ( t 0 ) = 2 m π + π / 2 ,then for all σ ∈ [ 1 / 2 ,1 ] , s i n ϕ ( σ , t ) increases along with the increasing of θ ( t ) ,at least on the interval defined by

2 m π + π / 2 + 14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) < θ ( t ) < 2 m π + 3 π / 2 . (55)

For another case with θ ( t 0 ) = 2 m π − π / 2 and θ ( t ) ∈ [ 2 m π − π / 2 ,2 m π + π / 2 ) ( m ≥ 1 ), the term − ∫ θ ( t 0 ) θ ( t )     F ( σ , x ( z , t 0 ) ) c o s z d z of I 1 ( σ , t ) decreases and dominates, and there is a considerable dominating range. In fact, the above approach holds and it follows from Equation (48) and Equation (50) that:

∂ ∂ Δ θ [ A ( σ , t 0 ) Δ I 1 − I 1 ( σ , t 0 ) Δ A A ( σ , t 0 ) + Δ A ] = A A + Δ A [ ∂ Δ I 1 ∂ Δ θ − I 1 + Δ I 1 A + Δ A ⋅ ∂ Δ A ∂ Δ θ ] ≤ A A + Δ A { − F ( 1 / 2 , e 2 π / t 0 ) Δ θ + 4 t 0 [ 2 G 1 ( 1, e w 1 ) log ( t 0 / 2 π ) + G 2 ( 1 / 2 , e w 2 ) ]       + 8 2 G 1 ( 1, e w 1 ) t 0 log ( t 0 / 2 π ) } < A A + Δ A [ − F ( 1 / 2 , e 2 π / t 0 ) Δ θ + 14 t 0 G 2 ( 1 / 2 , e w 2 ) ] < 0 (56)

provided that

Δ θ > 14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) .

For this case, the result is as follows:

Theorem 4. For the case with

θ ( t ) = t 4 l o g ( t 2 π e ) + π 6 ∈ [ 2 m π − π / 2 ,2 m π + π / 2 )

( m is a positive integer), let θ ( t 0 ) = 2 m π − π / 2 , then for all σ ∈ [ 1 / 2 ,1 ] , s i n ϕ ( σ , t ) decreases along with the increasing of θ ( t ) ,at least on the interval defined by

2 m π − π / 2 + 14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) < θ ( t ) < 2 m π + π / 2 .

In the following we consider the relative change of s i n ϕ ( σ , t ) on the remainder interval [ ( k − 1 / 2 ) π , ( k − 1 / 2 ) π + p ( t 0 ) ] ( k is a positive integer) with

p ( t 0 ) = 14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) ,

where t 0 is defined by θ ( t 0 ) = ( k − 1 / 2 ) π . We firstly study the odd case with k = 2 m + 1 ( m ≥ 0 ). At this time, the inequality s i n ϕ ( σ , t 0 ) < 0 holds for all σ ∈ [ 1 / 2 ,1 ] . According to Equation (52),

s i n ϕ ( σ , t ) − s i n ϕ ( σ , t 0 ) = 1 A ( σ , t 0 ) + Δ A [ Δ I 1 − I 1 ( σ , t 0 ) A ( σ , t 0 ) Δ A ] ,

whose sign is determined by the factor in the bracket. In the following we consider this part. With the denotations

P ( σ , e 2 z * / t 0 ) = 2 l o g ( t 0 / 2 π ) G 1 ( σ , e 2 z * / t 0 ) + G 2 ( σ , e 2 z * / t 0 ) ,

Q ( σ , t 0 ) = H ( σ , t 0 ) H t ( σ , t 0 ) t 0 + S ( σ , t 0 ) S t ( σ , t 0 ) t 0 H 2 ( σ , t 0 ) + S 2 ( σ , t 0 ) , (57)

to recall the approximations in Equation (45) and Equation (50), we have

Δ I 1 − I 1 ( σ , t 0 ) A ( σ , t 0 ) Δ A ≈ ∫ 0 Δ θ     F ( σ , e 2 z * / t 0 ) sin z * d z * − 2 Δ θ t 0 ∫ Δ θ ∞     P ( σ , e 2 z * / t 0 ) sin z * d z *     − I 1 ( σ , t 0 ) A ( σ , t 0 ) ⋅ 4 Δ θ t 0 log ( t 0 / 2 π ) Q ( σ , t 0 ) = ∫ 0 Δ θ [ F ( σ , e 2 z * / t 0 ) + 2 Δ θ t 0 P ( σ , e 2 z * / t 0 ) ] sin z * d z *     − 2 Δ θ t 0 [ ∫ 0 ∞     P ( σ , e 2 z * / t 0 ) sin z * d z * + 2 sin ϕ ( σ , t 0 ) log ( t 0 / 2 π ) Q ( σ , t 0 ) ] . (58)

Notice that for 0 < Δ θ ≤ p ( t 0 ) < π / 2 the first term is positive, and its sign relies on the second term. Since for a fixed σ , the function

Φ ( σ , t 0 ) = ∫ 0 ∞     P ( σ , e 2 z * / t 0 ) s i n z * d z * + 2 s i n ϕ ( σ , t 0 ) l o g ( t 0 / 2 π ) Q ( σ , t 0 ) (59)

is only associated with t 0 , it can be used to make judgement. In case Φ ( σ , t 0 ) ≤ 0 , Equation (58) implies that s i n ϕ ( σ , t ) > s i n ϕ ( σ , t 0 ) . The thing becomes simple, and θ ( t 0 ) = ( 2 m + 1 / 2 ) π accords with the minimum extreme point of s i n ϕ ( σ , t ) . In case Φ ( σ , t 0 ) > 0 the previous estimations make sense. To take the partial derivative about Δ θ , it reads

∂ ∂ Δ θ [ Δ I 1 − I 1 ( σ , t 0 ) A ( σ , t 0 ) Δ A ] ≈ [ F ( σ , e 2 Δ θ / t 0 ) + 2 Δ θ t 0 P ( σ , e 2 Δ θ / t 0 ) ] s i n Δ θ       + 2 t 0 ∫ 0 Δ θ     P ( σ , e 2 z * / t 0 ) s i n z * d z * − 2 t 0 Φ ( σ , t 0 ) . (60)

Corresponding to the case with Δ θ = 0 , its value is − ( 2 / t 0 ) Φ ( σ , t 0 ) < 0 . On the one hand, from the previous sections we know F , G 1 and G 2 are all positive functions. On the other hand, for fixed σ ∈ [ 1 / 2 ,1 ] , G 1 ( σ , e w ) and G 2 ( σ , e w ) increase along with the increasing of w on [ 0, w 1 ] and [ 0, w 2 ] (with w 1 ≈ 0.6484 , w 2 ≈ 0.4746 ), respectively. Notice that 0 ≤ Δ θ < π / 2 and t 0 ≈ 20.8747 for the initial case with m = 0 , we have w = 2 Δ θ / t 0 < 2 × ( π / 2 ) / 20.8747 ≈ 0.1505 < w 1 , w 2 . So P ( σ , e 2 Δ θ / t 0 ) increases along with the increasing of Δ θ . In addition, though F ( σ , e 2 Δ θ / t 0 ) decreases along with the increasing of Δ θ , the product of s i n Δ θ and it increases. So the function in Equation (60) increases along with the increasing of Δ θ on [ 0, p ( t 0 ) ] for fixed σ and t 0 . This indicates that Δ I 1 − I 1 Δ A / A decreases first and possibly increases later on this small interval. After all, Δ θ = p ( t 0 ) accords with the left end of the increasing interval of s i n ϕ ( σ , t ) − s i n ϕ ( σ , t 0 ) .

So, on the remainder interval [ ( 2 m + 1 / 2 ) π , ( 2 m + 1 / 2 ) π + p ( t 0 ) ] ( m ≥ 0 ), the variation of θ ( t ) also makes significant influence on s i n ϕ ( σ , t ) , and the minimum extreme point shifts to the right or not is determined by the sign of Φ ( σ , t 0 ) . Furthermore, from Equation (57) and Equation (59) we see, due to the existence of decay factor 2 / l o g ( t 0 / 2 π ) in the other two terms,

∫ 0 ∞     G 2 ( σ , e 2 z * / t 0 ) s i n ( z * ) d z *

dominates the sign of Φ ( σ , t 0 ) , and influences the location-shifting of the minimum extreme points. For the even case with k = 2 m ( m ≥ 1 ), to repeat the above deduction process we get a symmetric result for the maximum extreme points.

The correspondingly corollary is as follows:

Corollary 2. For the case with θ ( t ) ∈ [ 2 m π ,2 ( m + 1 ) π ] ( m is a positive integer), let θ ( t ˜ m ) = 2 m π + π / 2 and θ ( t ^ m ) = 2 m π + 3 π / 2 , then there is a minimum extreme point and a maximum extreme point of s i n ϕ ( σ , t ) , which lie separately in the intervals [ 2 m π + π / 2 ,2 m π + π / 2 + p ( t ˜ m ) ] and [ 2 m π + 3 π / 2 ,2 m π + 3 π / 2 + p ( t ^ m ) ] , where the function p is defined by

p ( t ) = 14 G 2 ( 1 / 2 , e w 2 ) t F ( 1 / 2 , e 2 π / t ) . (61)

As m → ∞ , p ( t ˜ m ) , p ( t ^ m ) → 0 , and these two extreme points in the moving intervals tend to their left boundaries t ˜ m and t ^ m .

3.3.4. The Explicit Periodic Form for U ( σ , t )

It follows from Theorem 3 and Theorem 4 that, for a fixed σ ∈ [ 1 / 2 ,1 ] , s i n ϕ ( σ , t ) = I 1 ( σ , t ) / A ( σ , t ) increases or decreases along with the increasing of θ ( t ) on most of the variation range [ ( k − 1 / 2 ) π , ( k + 1 / 2 ) π ) ( k is a positive integer). In case k is odd it increases, and in case k is even it decreases. This indicates that θ ( t ) is the main part of ϕ ( σ , t ) . It follows from Corollary 2 that, As m → ∞ , θ ( t ) = 2 m π + π / 2 accords with the minimum extreme point of s i n ϕ ( σ , t ) , and the phase angle needs to be adjusted by adding one π. In addition to this, in ϕ ( σ , t ) there should be also a perturbation term, say − ε ( σ , t ) . It reflects the contributions from other factors. We note that the negative sign here is to accord with the fact that all the location-shifts of the extreme points are to the right. Based on this, for every given σ ∈ [ 1 / 2 ,1 ] and every t ≥ t 0 with θ ( t 0 ) = π / 2 ( t 0 ≈ 20.8747 ), we have an explicit formula bellow:

s i n ϕ ( σ , t ) = s i n [ θ ( t ) + π − ε ( σ , t ) ] . (62)

To recall the relation in Equation (32), we express the real part of ξ ( σ + i t ) in an explicit periodic form:

U ( σ , t ) = 2 t A ( σ , t ) s i n [ ϕ ( σ , t ) + θ ( t ) ] = 2 t A ( σ , t ) s i n [ 2 θ ( t ) + π − ε ( σ , t ) ] = 2 t A ( σ , t ) s i n { 2 [ t 4 l o g ( t 2 π e ) + π 6 ] + π − ε ( σ , t ) } = A * ( σ , t ) s i n [ t 2 l o g ( t 2 π e ) + 4 3 π − ε ( σ , t ) ] (63)

with A * ( σ , t ) = 2 A ( σ , t ) / t , which matches exactly with the anticipation in Section 2 for the particular case with σ = 1 / 2 .

In the following we give some estimations to the amplitude function A * ( σ , t ) and the perturbation function ε ( σ , t ) .

3.3.5. To Refine the Estimation of Amplitude Function A ∗ ( σ , t )

From Equation (31) we see 0 < A ( σ , t ) < 5 F 0 ( 1 ) . In fact, the lower bound can be refined. To denote x ( z , t ) = e 2 z / t , then it follows from Equation (24) and Equation (38) that

F z ( σ , x ( z , t ) ) = F x ( σ , x ) ⋅ 2 t e 2 z / t = 2 t F x ( σ , x ) x = − 2 t G 2 ( σ , e 2 z / t ) < 0.

Based on this, for a fixed σ ∈ [ 1 / 2 ,1 ] we refine the estimation in Equation (27) as follows:

∫ 0 ∞     F ( σ , e 2 z / t ) sin ( z ) d z = 2 ∑ k = 0 ∞ [ F ( σ , e 2 ( z ′ k + 2 k π ) / t ) − F ( σ , e 2 ( z ′ k + 2 k π + π ) / t ) ] = − 2 ∑ k = 0 ∞   ∫ z ′ k + 2 k π z ′ k + ( 2 k + 1 ) π     F z ( σ , x ( z , t ) ) d z = 4 t ∑ k = 0 ∞   ∫ z ′ k + 2 k π z ′ k + ( 2 k + 1 ) π     G 2 ( σ , e 2 z / t ) d z = 4 π t ∑ k = 0 ∞     G 2 ( σ , e w ^ k ) , (64)

in which 0 ≤ z ′ k ≤ π and

4 k π t ≤ 2 t ( z ′ k + 2 k π ) ≤ w ^ k ≤ 2 t [ z ′ k + ( 2 k + 1 ) π ] ≤ 4 ( k + 1 ) π t . (65)

In view of Figure 5, for a given σ ∈ [ 1 / 2 ,1 ] , G 2 ( σ , e w ) has a unique maximum extreme point at w 2 . On [ 0, w 2 ] it increases and on [ w 2 , ∞ ) it decreases. Along with the increasing of σ , the location of the maximum extreme point is shifted from about ( 0.4746,1.2304 ) to about ( 0.4820,1.2237 ) . Notice that t ≥ t 0 ≈ 20.8747 , 4 π / t ≤ 4 π / t 0 ≈ 0.6020 > w 2 , it follows from Equation (26) and Equation (64) that:

A ( σ , t ) = [ ∫ 0 ∞     F ( σ , e 2 z / t ) c o s ( z ) d z ] 2 + [ ∫ 0 ∞     F ( σ , e 2 z / t ) s i n ( z ) d z ] 2 ≥ ∫ 0 ∞     F ( σ , e 2 z / t ) s i n ( z ) d z = 4 π t ∑ k = 0 ∞     G 2 ( σ , e w ^ k ) > 4 π t G 2 ( σ , e w ^ 1 ) ≥ 4 π t G 2 ( 1 / 2 , e 8 π / t 0 ) ≈ 4 π t ⋅ 0.115 > 2 π 5 t .

To combine this result with that in Equation (31) we get 2 π / 5 t < A ( σ , t ) < 5 F 0 ( 1 ) and hence

4 π 5 t 2 < A * ( σ , t ) < 2 5 F 0 ( 1 ) t , (66)

which has clarified the decaying characteristics of U ( σ , t ) with respect to t.

3.3.6. To Estimate the Perturbation Function ε ( σ , t )

In the following we explore the perturbation term ε ( σ , t ) .

Firstly, for every given σ ∈ [ 1 / 2 ,1 ] , we affirm that the case with

l i m t → ∞ | ε ( σ , t ) | = ∞

is impossible. On the contrary, if it is true, there must be a large enough t * and an integer M such that ε ( σ , t * ) = ( 2 M + 1 ) π . At this time, θ ( t * ) lies either in [ 2 m π + π / 2 ,2 m π + 3 π / 2 ) or in [ 2 m π + 3 π / 2 ,2 m π + 5 π / 2 ) for some m ≥ 0 . In the following we consider the first half period. In case θ ( t * ) lies in the increasing interval, that is, 2 m π + π / 2 + p ( t 0 ) < θ ( t * ) < 2 m π + 3 π / 2 [here p ( t 0 ) is defined in Equation (61) with θ ( t 0 ) = 2 m π + π / 2 ], s i n ϕ ( σ , t ) must increases in the neighborhood of t = t * . Yet at this time,

s i n ϕ ( σ , t * ) = s i n [ θ ( t * ) + π − ( 2 M + 1 ) π ] = s i n θ ( t * )

decreases. This leads to a contradiction. In case θ ( t * ) lies in the remainder interval 2 m π + π / 2 ≤ θ ( t * ) ≤ 2 m π + π / 2 + p ( t 0 ) , there should be s i n ϕ ( σ , t * ) = s i n θ ( t * ) > 0 . Yet, according to the previous arguments,

s i n ϕ ( σ , t 0 ) = I 1 ( σ , t 0 ) / A ( σ , t 0 ) < 0.

As t 0 → ∞ , p ( t 0 ) → 0 and in case t * is large enough, it is sufficiently close to t 0 , and it follows from the continuousness of the function that s i n ϕ ( σ , t * ) < 0 . This is also a contradiction. On another half period with 2 m π + 3 π / 2 ≤ θ ( t * ) < 2 m π + 5 π / 2 , a contradiction is also met. This indicates that ε ( σ , t ) is bounded.

Secondly, on every period there exists an interval on which | ε ( σ , t ) | < π / 2 . As the case θ ( t 0 ) = 2 m π + π / 2 concerned, it follows from the previous arguments that, either t 0 is exactly the minimum extreme point of s i n ϕ ( σ , t ) , or there is a right-shifting of the minimum extreme point. For any given σ ∈ [ 1 / 2 ,1 ] , for the later case the extreme point accords with

ϕ ( σ , t 1 ) = θ ( t 1 ) + π − ε ( σ , t 1 ) = 2 m π + 3 π / 2

for some t = t 1   ( > t 0 ) . According to Corollary 2, θ ( t 1 ) ∈ ( 2 m π + π / 2 ,2 m π + π / 2 + p ( t 0 ) ] . Hence,

0 < ε ( σ , t 1 ) = θ ( t 1 ) − θ ( t 0 ) ≤ p ( t 0 ) = 14 G 2 ( 1 / 2 , e w 2 ) t 0 F ( 1 / 2 , e 2 π / t 0 ) < 2 5 π . (67)

Due to the existence of the decay factor 1 / t 0 , its value should become smaller and smaller as m → ∞ . This indicates that, in any given half period [ 2 m π + π / 2 ,2 m π + 3 π / 2 ] , there indeed exists some region such that | ε ( σ , t ) | < π / 2 . This result also holds for another half period.

Thirdly, the inequality | ε ( σ , t ) | < π / 2 holds for big enough t.

For every given σ ∈ [ 1 / 2 ,1 ] , the extreme point of s i n ϕ ( σ , t ) is unique with respect to ϕ ( σ , t ) ∈ ( k π , ( k + 1 ) π ) ( k is a nonnegative integer), if and only if ϕ ( σ , t ) is increasing along with the increasing of t, that is,

ϕ t ( σ , t ) = 1 4 l o g ( t 2 π ) − ε t ( σ , t ) > 0. (68)

Since ε ( σ , t ) is bounded, ε t ( σ , t ) should be also bounded, say | ε t ( σ , t ) | ≤ C 0 , where C 0 is a positive constant. Certainly, the above inequality holds for the case with ε t ( σ , t ) ≤ 0 . For another case with ε t ( σ , t ) > 0 , the estimation makes sense, and it holds provided that

t > T = 2 π e 4 C 0 . (69)

Under the above assumption, s i n ϕ ( σ , t ) has a unique negative minimum extreme point and a unique positive maximum extreme point on [ 2 m π ,2 ( m + 1 ) π ) with respect to ϕ ( σ , t ) for any big enough m . As the minimum extreme point concerned, according to Corollary 2, it lies in the interval [ 2 m π + π / 2 ,2 m π + π / 2 + p ( t 0 ) ] with respect to θ ( t ) , where t 0 is defined by θ ( t 0 ) = 2 m π + π / 2 . Similarly, the maximum extreme point lies in the interval [ 2 m π + 3 π / 2 ,2 m π + 3 π / 2 + p ( t 0 ) ] with respect to θ ( t ) , where t 0 is defined by θ ( t 0 ) = 2 m π + 3 π / 2 .

If the assertion is false, there must be a t * > T and a σ * ∈ [ 1 / 2 ,1 ] such that ε ( σ * , t * ) = π / 2 or − π / 2 . Let t 1 ,   t 2 ,   t 3 ,   t 4 and t 5 be the points which accord with ϕ ( σ * , t ) = ( 2 m + 1 ) π , ( 2 m + 1 ) π + π / 2 , 2 ( m + 1 ) π , 2 ( m + 1 ) π + π / 2 and ( 2 m + 3 ) π , respectively. Firstly, t * ≠ t 2 ,   t 4 , since they are the extreme points with 0 < ε ( σ , t * ) < 2 π / 5 for all σ ∈ [ 1 / 2 ,1 ] . Secondly, t * = t 1 ,   t 3 ,   t 5 are also impossible. The reason is that, for the case with t * = t 3 , s i n ϕ ( σ * , t * ) = s i n [ 2 ( m + 1 ) π ] = 0 , yet θ ( t 3 ) = ϕ ( σ * , t 3 ) − π + ε ( σ * , t 3 ) = ( 2 m + 1 ) π ± π / 2 accords with s i n ϕ ( σ * , t 3 ) ≠ 0 (see Figure 7). This is a contradiction. The same thing occurs for t * = t 1 and t 5 .

In the following we firstly consider the case with ε ( σ * , t * ) = π / 2 . In case t * lies in the increasing interval ( t 2 , t 3 ) of s i n ϕ ( σ * , t ) , then

s i n ϕ ( σ * , t * ) = s i n [ θ ( t * ) + π − ε ( σ * , t * ) ] = s i n [ ( θ ( t * ) − π / 2 ) + π ] = s i n [ θ ( t ˜ ) + π ]

with t ˜ ∈ ( t 1 , t 2 ) . Yet from Theorem 4 we see on this interval s i n ϕ ( σ * , t ) decreases along with the increasing of θ ( t ) , that is, along with the increasing of t. This is a contradiction. In case t * ∈ ( t 3 , t 4 ) , we have s i n ϕ ( σ * , t ) > 0 . Yet at this time,

s i n ϕ ( σ * , t * ) = s i n [ ( θ ( t * ) − π / 2 ) + π ] = s i n [ θ ( t ˜ ) + π ] < 0

with t ˜ ∈ ( t 2 , t 3 ) . This is also a contradiction. In case t * lies in the decreasing interval ( t 4 , t 5 ) of s i n ϕ ( σ * , t ) , then it follows from Theorem 3 that, there is also a corresponding point t ˜ ∈ ( t 3 , t 4 ) , at which s i n ϕ ( σ * , t ) increases. A contradiction is met. In case t * ∈ ( t 1 , t 2 ) , we have s i n ϕ ( σ * , t * ) < 0 , and there should be a corresponding point t ˜ ∈ ( t 0 , t 1 ) with ϕ ( σ * , t 1 ) = ( 2 m + 1 ) π − π / 2 , such that s i n ϕ ( σ * , t * ) > 0 , and get a contradiction.

For the case with ε ( σ * , t * ) = − π / 2 , to repeat the above processes we get the corresponding contradictions. In summary, for all σ ∈ [ 1 / 2 ,1 ] and all t > T = 2 π e 4 C 0 , it holds the estimation | ε ( σ , t ) | < π / 2 .

In addition, along with the increasing of θ ( t ) on ( π / 2 + p ( t 0 ) , 3 π / 2 ) , s i n ϕ ( σ , t ) increases. This indicates that θ ( t ) dominates and the inequality in Equation (68) holds true. Notice that p ( t 0 ) < 2 π / 5 , to solve ( t * / 4 ) l o g ( t * / 2 π e ) + π / 6 = π / 2 + 2 π / 5 , it yields t * ≈ 26.2694 . Hence, on the inverse manner, from Equation (68) we have

C 0 < 1 4 log ( t * 2 π ) ≈ 0.3576 < 2 5 .

This partly indicates that, the requirement for T is not too high.

In the following we refine the estimation on ε ( σ , t ) with another approach.

Firstly, as the extreme points of s i n ϕ ( σ , t ) concerned, for a given σ ∈ [ 1 / 2 ,1 ] , there must be a positive integer k and a t k * such that θ ( t k * ) + π − ε ( σ , t k * ) = k π + 3 π / 2 , from which we get [same to Equation (67)]

0 ≤ ε ( σ , t k * ) ≤ p ( t ¯ k ) = 14 G 2 ( 1 / 2 , e w 2 ) t ¯ k F ( 1 / 2 , e 2 π / t ¯ k ) ≤ C ¯ t ¯ k (70)

with C ¯ = 14 G 2 ( 1 / 2 , e w 2 ) / F ( 1 / 2 , e 2 π / t ¯ 1 ) , where t ¯ k is defined by θ ( t ¯ k ) = k π + π / 2 . Since this inequality holds for all k , we expect the estimation

| ε ( σ , t ) | ≤ C 1 t (71)

holds for some finite positive constant C 1 which does not rely on σ .

On the contrary, if for any given C > 0 there is always a σ = σ 0 and a t = T such that | ε ( σ 0 , T ) | > C / T > 0 , then there is an integer K > 0 such that θ ( t K * ) , θ ( T ) ∈ [ K π + π / 2 ,   K π + 3 π / 2 ) . In the following we seek for a contradiction with case by case discussion.

Case a: In case ε ( σ 0 , t K + 1 * ) > 0 , Since both | ε ( σ 0 , T ) | and ε ( σ 0 , t K + 1 * ) are positive, there must be a finite proportional coefficient λ > 0 such that | ε ( σ 0 , T ) | = λ ε ( σ 0 , t K + 1 * ) . In view of t ¯ K + 1 > T we have

C t ¯ K + 1 < C T < | ε ( σ 0 , T ) | = λ ε ( σ 0 , t K + 1 * ) ≤ λ p ( t ¯ K + 1 ) ≤ λ C ¯ t ¯ K + 1 .

A contradiction is met only if we choose C > λ C ¯ .

Case b: In case ε ( σ 0 , t K + 1 * ) = 0 , if there is an integer n   ( > 1 ) such that ε ( σ 0 , t K + n * ) > 0 , then the previous contradiction is also met by choosing C > λ C ¯ . At this time, the inequality bellow

C t ¯ K + n < C T < | ε ( σ 0 , T ) | = λ ε ( σ 0 , t K + n * ) ≤ λ p ( t ¯ K + n ) ≤ λ C ¯ t ¯ K + n .

Case c: If for all k > K , the equality ε ( σ 0 , t k * ) = 0 always holds, furthermore, if there is a σ 1 ∈ [ 1 / 2 ,1 ] and an integer n   ( ≥ 1 ) such that ε ( σ 1 , t K + n * ) > 0 , then | ε ( σ 0 , T ) | = λ ε ( σ 1 , t K + n * ) also holds for some λ > 0 , and the previous contradiction is also met. We note that, due to the continuousness of ε ( σ , t ) , the occurrence of this situation is very natural.

Case d: For all σ ∈ [ 1 / 2 ,1 ] and all k > K , the equality ε ( σ , t k * ) = 0 always holds, then the particular point t ¯ k with θ ( t ¯ k ) = k π + π / 2 is exactly the extreme point of s i n ϕ ( σ , t ) . If that happens, the refined estimation in Equation (71) is invalid. However, the estimation | ε ( σ , t ) | < π / 2 is still valid for this case.

We note that, the ideal situation in the last case hardly occurs, since it requires the inequality Φ ( σ , t ¯ k ) ≤ 0 always holds. This, in turn, needs the dominating factor

∫ 0 ∞     G 2 ( σ , e 2 z * / t ¯ k ) s i n ( z * ) d z * ≤ 0.

Yet, take the case with σ = 1 / 2 as an example, G 2 ( 1 / 2 , e 2 z * / t ¯ k ) has peak value at about z * = t ¯ k w 2 / 2 ≈ 0.2373 t ¯ k , where the ratio 0.2373 is not an integer or half integer. Furthermore, since this peak value may dominate the corresponding integral, as a whole, the variation of Φ ( σ , t ¯ k ) is not synchronous with that of θ ( t ¯ k ) , and it prevents the occurrence of this ideal situation. Unfortunately, a rigid proof is unobtainable, and this case can not be excluded.

We note that, for Cases a-c, the requirement t > T = 2 π e 4 C 0 for the uniqueness of extreme point can be loosened. Since it holds | ε ( σ , t ) | ≤ C 1 / t , to express ε ( σ , t ) = h ( σ , t ) / t , then | h ( σ , t ) | ≤ C 1 and there is also a finite positive constant C 2 such that | h t ( σ , t ) | ≤ C 2 . Hence,

| ε t ( σ , t ) | = | ∂ ∂ t [ h ( σ , t ) t ] | = | h t ( σ , t ) t − h ( σ , t ) ⋅ 1 t 2 | = | h t ( σ , t ) − ε ( σ , t ) | t ≤ C 2 + C 1 t .

For this case, the inequality in Equation (68) holds on condition that t > T with

T 4 l o g ( T 2 π ) = C 1 + C 2 .

Certainly, there is another understanding. For Cases a-c, the estimation | ε ( σ , t ) | ≤ C 1 / t can be used in a direct way. To ensure | ε ( σ , t ) | < π / 2 , it only requires t > T = 2 C 1 / π .

3.3.7. The Finial Results for U and the Verification of Proposition 2

As stated in the previous section, for every given σ ∈ [ 1 / 2 ,1 ] , only ift is big enough, say t > T = 2 π e 4 C 0 , there is ϕ t ( σ , t ) = θ ′ ( t ) − ε t ( σ , t ) > 0 , and on every period the uniqueness of the minimum extreme point and maximum extreme point of s i n ϕ ( σ , t ) is ensured. To denote

ϑ ( σ , t ) = 2 θ ( t ) + 4 3 π − ε ( σ , t ) = t 2 log ( t 2 π e ) + 4 3 π − ε ( σ , t ) , (72)

then the inequality ϑ t ( σ , t ) = θ ′ ( t ) + ϕ t ( σ , t ) > 0 holds for t > T = 2 π e 4 C 0 , which is bigger than 2π.

For every permitted integer k , with respect to ϑ ( σ , t ) ∈ ( k π , ( k + 1 ) π ) , we verify the uniqueness of the extreme point of U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) . We note that, since U ( σ , t ) is the real part of the analytic function ξ ( σ + i t ) , both A * ( σ , t ) and ε ( σ , t ) are smooth functions which have arbitrary order derivatives.

For every given σ ∈ [ 1 / 2 ,1 ] , the partial derivative of U ( σ , t ) about t reads:

U t ( σ , t ) = ∂ ∂ t [ A * ( σ , t ) s i n ϑ ( σ , t ) ] = A t * ( σ , t ) s i n ϑ ( σ , t ) + A * ( σ , t ) ϑ t ( σ , t ) c o s ϑ ( σ , t ) = A * ( σ , t ) ϑ t ( σ , t ) s i n ϑ ( σ , t ) [ c o t ϑ ( σ , t ) − − A t * ( σ , t ) A * ( σ , t ) ϑ t ( σ , t ) ] = : A * ( σ , t ) ϑ t ( σ , t ) s i n ϑ ( σ , t ) [ c o t ϑ ( σ , t ) − g ( σ , t ) ] . (73)

We note that, as the extreme points concerned, the case s i n ϑ ( σ , t ) = 0 is excluded since at this time U t ( σ , t ) ≠ 0 . With this understanding, all the locations of extreme points are determined by the equation c o t ϑ ( σ , t ) = g ( σ , t ) . It is easily checked that, no matter g ( σ , t ) is positive or not, its curve intersects with that of c o t ϑ ( σ , t ) only once on every given period ( k π , ( k + 1 ) π ) . This indicates U ( σ , t ) has a unique extreme point on ( k π , ( k + 1 ) π ) .

Furthermore, notice that for a fixed σ , c o t ϑ ( σ , t ) → ∞ as t → k π from the right and c o t ϑ ( σ , t ) → − ∞ as t → ( k + 1 ) π from the left, and g ( σ , t ) = − A t * ( σ , t ) / A * ( σ , t ) ϑ t ( σ , t ) is bounded, for all k the function c o t ϑ ( σ , t ) − g ( σ , t ) always changes from positive to negative on ( k π , ( k + 1 ) π ) . In addition, in case k = 2 m it holds s i n ϑ ( σ , t ) > 0 , yet in case k = 2 m − 1 it holds s i n ϑ ( σ , t ) < 0 . It follows from Equation (73) that, on ( 2 m π , ( 2 m + 1 ) π ) the unique positive extreme point of U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) is the maximum one, and on ( ( 2 m − 1 ) π ,2 m π ) the unique negative extreme point is the minimum one.

Hence, on every half period, the uniqueness of the extreme point of U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) is verified. It only requires t > T = 2 π e 4 C 0 .

To summarize all the obtained results we have the final theorem:

Theorem 5. For every given σ ∈ [ 1 / 2 ,1 ] and every t ≥ t 0 with

θ ( t 0 ) = t 0 4 l o g ( t 0 2 π e ) + π 6 = π 2

( t 0 ≈ 20.8747 ), the real part of ξ ( σ + i t ) is in an explicit periodic form:

U ( σ , t ) = A * ( σ , t ) s i n [ t 2 l o g ( t 2 π e ) + 4 3 π − ε ( σ , t ) ] .

Here the amplitude decays within

4 π 5 t 2 < A * ( σ , t ) < 2 5 F 0 ( 1 ) t .

The perturbation term ε ( σ , t ) is bounded. ε t ( σ , t ) is also bounded, say | ε t ( σ , t ) | ≤ C 0 ,where C 0 is a positive constant. In case t > T = 2 π e 4 C 0 , there is an estimation | ε ( σ , t ) | < π / 2 ,and on every period U ( σ , t ) has a unique positive maximum extreme point and a unique negative minimum extreme point. In addition, to define t k * by θ ( t k * ) = k π + π / 2 ,then except the particular case with ε ( σ , t k * ) ≡ 0 for all σ ∈ [ 1 / 2 ,1 ] and all nonnegative integer k,there is a refined estimation:

| ε ( σ , t ) | ≤ C 1 t ,

here C 1 is a positive constant.

These theoretical results accord well with the numerical ones in Section 2 (see Figures 1-3), and the distribution law

t k 2 l o g ( t k 2 π e ) + 4 3 π = k π ,   k = 1,2, ⋯

of the zero points on the critical line σ = 1 / 2 is embodied perfectly. Here t k stands for the approximate value of the k -th zero point. This indicates Proposition 2 given by Riemann is true. The explicit illustration is as follows:

Theorem 6. For a big enough T which does not accord with a zero point, in the range 0 < t < T the number of zero points on the critical line σ = 1 / 2 obeys the exact formula bellow (understood as taking the integer part):

N ( T ) = T 2 π l o g ( T 2 π e ) + 4 3 − ε ( 1 / 2 , T ) π (74)

with | ε ( 1 / 2 , T ) / π | < 1 / 2 .

4. The Finial Proof of Riemann Hypothesis

To get an intuitive understanding on the variation of ξ ( σ + i t ) , we make numerical simulations on U ( σ , t ) and V ( σ , t ) according to the formulas in Equation (18) and Equation (23), that is,

U ( σ , t ) = ∫ 1 ∞     f ( σ , x ) cos ( t 2 l o g x ) d x , V ( σ , t ) = ∫ 1 ∞     g ( σ , x ) s i n ( t 2 l o g x ) d x (75)

with

f ( σ , x ) = 2 x 1 / 4 ( x δ + x − δ ) ∑ n = 1 ∞     a n ( a n x − 3 2 ) e − a n x ,

g ( σ , x ) = 2 x 1 / 4 ( x δ − x − δ ) ∑ n = 1 ∞     a n ( a n x − 3 2 ) e − a n x ,

where δ = ( σ − 1 / 2 ) / 2 and a n = π n 2 .

Essentially, they are infinite integrals with infinite series. How to make objective numerical simulations for them is still a question. Our try with finite approximations yield Figure 8 and Figure 9 which are just for reference. These figures roughly reflect the basic variation characteristics of U ( σ , t ) and V ( σ , t ) . We note that the simulated periodic signals are too strong, and they have depressed the random variations. There is an evidence for this, that is, the simulated zeros points of U ( 1 / 2 , t ) does not completely coincide with the known numerical ones (mentioned in Section 2).

Since ξ ( σ + i t ) is analytic, the well-known Cauchy-Riemann conditions hold

for its real and imaginary parts:

U σ = V t ,   U t = − V σ . (76)

One can check these with the two formulas in Equation (75) in a direct way.

Let Ω be any finite domain in two-dimensional real space ℝ 2 . The analytic property of ξ ( σ + i t ) implies good smoothness of U and V on Ω . So the se- cond-order partial derivatives of them exist and are continuous, that is, U , V ∈ C 2 ( Ω ) . It follows from Equation (76) that

U σ σ + U t t = 0 ,   V σ σ + V t t = 0 ,   ( σ , t ) ∈ Ω .

These mean bother U and V satisfy the two-dimensional Laplace equation, and the result below hold for them:

Lemma 2. (Extremum Principle)  : If u ( σ , t ) satisfies the Laplace equation u σ σ + u t t = 0 in Ω , then it has no extreme point in the interior of Ω , unless it is a constant on the entire region Ω ¯ .

It is well-known that there is no zero point outside the critical strip 0 ≤ σ ≤ 1 (see      ). Notice that the zero points of ξ ( σ + i t ) accord with the nontrivial zero points of ζ ( σ + i t ) , as the real-valued problem concerned, the equivalent statement of the Riemann hypothesis is as follows:

Except on the critical line σ = 1 / 2 , U ( σ , t ) and V ( σ , t ) have no other mutual zero point in the region [ 0,1 ] × ( 0, ∞ ) .

Just as mentioned in Section 1, for any given finite region bounded by 0 ≤ σ ≤ 1 and 0 < t < T , the check can be always done. What lacks is the knowledge for the case with t → ∞ . Till now the known record was set by X. Gourdon in 2004 . He had checked the first 1013 zero points and found that they all possess real part σ = 1 / 2 . To speculate with our distribution law, its upper bound is up to T = 2.4460 × 10 12 . As reviewed in , another calculation had ever been done by Odlyzko in 1989  in selected intervals show that the Riemann hypothesis holds for over 3 × 108 zeros at height up to T = 2 × 10 20 . Our aim is to settle down the infinite case, and the lower bound T may be much lower than this. The following is the final theoretical result:

Theorem 7. Except on the critical line σ = 1 / 2 , the real part U ( σ , t ) and imaginary part V ( σ , t ) of ξ ( σ + i t ) have no other mutual zero point in the region [ 0,1 ] × ( T , ∞ ) ∪ [ 0,1 ] × ( − ∞ , − T ) for big enough T, say T = m a x { 2 π e 4 C 0 , 4 C 1 / π } , where C 0 is the upper bound of the bounded function | ε t ( σ , t ) | , C 1 is the positive constant which ensures | ε ( σ , t ) | ≤ C 1 / t .

Proof. It follows from Lemma 1 that, U ( σ , t ) is symmetric about the horizontal line t = 0 and vertical line σ = 1 / 2 , and V ( σ , t ) is anti-symmetric about the horizontal line t = 0 and vertical line σ = 1 / 2 . Hence, as the mutual zero point of U and V concerned, it only needs to consider one quarter part of the concerned region. Here we choose the one defined by ( 1 / 2 ,1 ] × ( T , ∞ ) , which excludes the critical line σ = 1 / 2 .

Suppose U ( σ , t ) and V ( σ , t ) has a mutual zero point in ( 1 / 2 ,1 ] × ( T , ∞ ) , say ( σ 0 , t 0 ) , then it must not be an isolated zero point of U ( σ , t ) or V ( σ , t ) . On the contrary, if so, to draw a small circle around ( σ 0 , t 0 ) , then U ( σ , t ) or V ( σ , t ) maintains its sign on this circle, either positive or negative. For this case, at the interior point ( σ 0 , t 0 ) , U = 0 or V = 0 , and this point must be a minimum extreme one or a maximum extreme one of U ( σ , t ) or V ( σ , t ) on this circular region. This contradicts with the Extremum Principle. Hence, ( σ 0 , t 0 ) is not isolated, and the 2-dimensional surfaces defined by U ( σ , t ) and V ( σ , t ) must intersect the zero-valued plane with two curves, say l 1 and l 2 , which pass through the mutual zero point ( σ 0 , t 0 ) . For convenience, we call them “zero-valued curves” of U and V.

Firstly, we assert that all the zero-valued curves of U intersect with the critical line σ = 1 / 2 . It follows from the Extremum Principle that all the zero points of ξ ( σ + i t ) on σ = 1 / 2 can not be isolated ones. Explicitly, since V ( 1 / 2 , t ) ≡ 0 holds for all t ∈ ( − ∞ , ∞ ) , this line itself is a zero-valued curve of V. Meanwhile, for every given positive integer k , let ( 1 / 2 , t k ) be the k -th zero point of U ( 1 / 2 , t ) , then there is also a zero-valued curve across it. This curve is also unique. In fact, it follows from Theorem 5 that U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) , and this curve is actually one of the intersecting curves of this continuous periodic surface and the zero-valued plane, which crosses the point ( 1 / 2 , t k ) . Exactly, this zero-valued curve is determined by

ϑ ( σ , t ) = 2 θ ( t ) − ε ( σ , t ) = t 2 l o g ( t 2 π e ) + 4 3 π − ε ( σ , t ) = k π (77)

with | ε ( σ , t ) | < π / 2 . Under the condition that t > T ≥ 2 π e 4 C 0 , the inequality ϑ t ( σ , t ) > 0 holds, and for the case with σ ∈ [ 1 / 2 ,1 ] , the value of t is uniquely determined. The periodicity of U ( σ , t ) ensured that there is no other zero-valued curves. Hence, the zero-valued curve l 1 of U ( σ , t ) , which is across ( σ 0 , t 0 ) , should accord with Equation (77) with some k = k 0 .

Secondly, there exists a horizontal strip which includes l 1 , on which the sign of V σ maintains unchanged, and a contradiction is met. In the following we make case by case discussion.

For the particular case with ε ( σ , t 0 ) ≡ 0 for all σ ∈ [ 1 / 2 ,1 ] , the point t 0 accords with θ ( t 0 ) = m π + π / 2 for some integer m, and hence

ϑ ( σ , t 0 ) = 2 θ ( t 0 ) − ε ( σ , t 0 ) = 2 ( m π + π / 2 ) − 0 = ( 2 m + 1 ) π .

At this time, l 1 is a horizontal straight line which does not rely on σ . In view of U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) , along with the increasing of σ ∈ [ 1 / 2 ,1 ] , with respect this zero-valued line, the neighboring maximum extreme points and minimum extreme points are also linked into two curves, say h 1 and h 2 , which satisfy

ϑ ( σ , t _ ) = 2 θ ( t _ ) − ε ( σ , t _ ) = 2 m π + π / 2 , ϑ ( σ , t ¯ ) = 2 θ ( t ¯ ) − ε ( σ , t ¯ ) = 2 m π + 3 π / 2 , (78)

respectively. Since t > T ≥ 2 π e 4 C 0 , it follows from Theorem 5 that | ε ( σ , t _ ) | < π / 2 and | ε ( σ , t ¯ ) | < π / 2 . Hence, 2 θ ( t _ ) = 2 m π + π / 2 + ε ( σ , t _ ) < ( 2 m + 1 ) π and 2 θ ( t ¯ ) = 2 m π + 3 π / 2 + ε ( σ , t ¯ ) > ( 2 m + 1 ) π . This indicates that 2 θ ( t _ ) < 2 θ ( t 0 ) < 2 θ ( t ¯ ) , from which we get t _ < t 0 < t ¯ for all σ ∈ [ 1 / 2 ,1 ] . Actually, h 1 and h 2 possess the coordinate ( σ , t _ ) and ( σ , t ¯ ) , respectively.

Certainly, l 1 is included in the horizontal strip delimited by h 1 and h 2 . If there is a zero-valued curve l 2 of V ( σ , t ) which intersects with l 1 at ( σ 0 , t 0 ) , then V ( 1 / 2 , t 0 ) = V ( σ 0 , t 0 ) = 0 , and along l 1 there must be an extreme point σ = σ * , at which V σ = 0 . Yet notice that for every given σ ∈ [ 1 / 2 ,1 ] , with respect to ϑ ( σ , t ) on ( 2 m π + π / 2 ,2 m π + 3 π / 2 ) , U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) decreases, there should be U t < 0 on the whole strip. Furthermore, it follows from the Cauchy-Riemann condition in Equation (76) that V σ = − U t > 0 . This leads to a contradiction.

Except the particular case ε ( σ , t 0 ) ≡ 0 , it follows from Theorem 5 that, the refined estimation | ε ( σ , t ) | ≤ C 1 / t holds for some positive constant C 1 . On condition that t > T ≥ 4 C 1 / π , we have

| ε ( σ , t ) | < π 4 .

To take 2 θ ( t ) as a research object, from Equation (77) we see l 1 accords with

k 0 π − π 4 < 2 θ ( t ˜ ) = k 0 π + ε ( σ , t ˜ ) < k 0 π + π 4 (79)

for some positive integer k 0 . At this time, the neighboring two extreme-point curves, say h 1 and h 2 , are defined as in Equation (78), which yield

2 θ ( t _ ) = k 0 π − π / 2 + ε ( σ , t _ ) < k 0 π − π 4 , 2 θ ( t ¯ ) = k 0 π + π / 2 + ε ( σ , t ¯ ) > k 0 π + π 4 . (80)

As a default, here t _ > T is fulfilled. To combine the inequalities in Equation (79) we have

2 θ ( t _ ) < k 0 π − π 4 < 2 θ ( t ˜ ) < k 0 π + π 4 < 2 θ ( t ¯ ) . (81)

In the horizontal strip delimited by h 1 and h 2 , there are two horizontal lines, say j 1 and j 2 , defined by 2 θ ( t ) = k 0 π − π / 4 and 2 θ ( t ) = k 0 π + π / 4 , respectively. If there is a zero-valued curve l 2 of V ( σ , t ) which intersects with l 1 at ( σ 0 , t 0 ) , then the horizontal line t = t 0 is between the two horizontal lines j 1 and j 2 . Naturally, we have V ( 1 / 2 , t 0 ) = V ( σ 0 , t 0 ) = 0 , and along the line t = t 0 there should be an extreme point σ = σ * , at which V σ = 0 . Yet in the horizontal strip delimited by h 1 and h 2 , ϑ ( σ , t ) lies in the monotone interval ( k 0 π − π / 2 , k 0 π + π / 2 ) . In case k 0 is odd, say k 0 = 2 m 0 + 1 , for every given σ ∈ [ 1 / 2 ,1 ] , U ( σ , t ) = A * ( σ , t ) s i n ϑ ( σ , t ) decreases on this interval. Particularly, between the two horizontal lines j 1 and j 2 , there is U t < 0 . Furthermore, it follows from the Cauchy-Riemann condition in Equation (76) that V σ = − U t > 0 . This leads to a contradiction. In case k 0 is even, say k 0 = 2 m 0 , ( k 0 π − π / 2 , k 0 π + π / 2 ) is the increasing interval of U ( σ , t ) , and between the two horizontal lines j 1 and j 2 there must be V σ = − U t < 0 . This is also contradictory.

The contradiction above indicate that the real part U ( σ , t ) and imaginary part V ( σ , t ) of ξ ( σ + i t ) have no mutual zero point in the region ( 1 / 2 ,1 ] × ( T , ∞ ) for big enough T. Furthermore, it follows from the symmetric properties of U and V that, except on the critical line σ = 1 / 2 , they have no other mutual zero point in the whole region [ 0,1 ] × ( T , ∞ ) ∪ [ 0,1 ] × ( − ∞ , − T ) . The proof is finished.

5. Conclusions and Remarks

The Riemann hypothesis is a well-known unsolved problem. Its difficulty lies in the complexity of the Zeta function ζ ( s ) . Essentially, it is involved in an infinite integral which includes infinite series with complex variables.

From October 2018 we had tried to detour this infinite integral, and explore an abstract proof with the “symmetry” of ξ ( s ) [which is equivalent to ζ ( s ) ], the Extremum Principle and some new techniques. After about one year, we had found that all these tries are in vain and turned to the “monotonicity”. As depicted in Figure 9, V ( σ , t ) shows good horizontal monotonicity. Yet the proof is not available. In 2020 we had grasped the key point, that is, to prove the Riemann hypothesis it requires, U ( 1 / 2 , t ) has no positive minimum extreme point and negative maximum extreme point. Yet, to prove this is like climbing the steep precipices and cliffs, there is no viable path! Until March 2021, when we began to investigate the distribution law of zero points on the critical line σ = 1 / 2 , a bold idea appeared, that is the “periodicity”! With this understanding, all the scattered results in the literatures are linked together. With three-month effort, a brand-new artwork appeared.

For every given σ ∈ [ 1 / 2 ,1 ] , we have found that the real part of ξ ( σ + i t ) is actually in an explicit periodic form:

U ( σ , t ) = A * ( σ , t ) s i n [ t 2 l o g ( t 2 π e ) + 4 3 π − ε ( σ , t ) ] .

with 4 π / 5 t 2 < A * ( σ , t ) < 2 5 F 0 ( 1 ) / t and | ε ( σ , t ) | < π / 2 . Except a particular case, there is also a refined estimation | ε ( σ , t ) | ≤ C 1 / t for some positive constants C 1 .

These theoretical results have verified the observed distribution law

t k 2 log ( t k 2 π e ) + 4 3 π = k π ,   k = 1 , 2 , ⋯

for the zero points on the critical line σ = 1 / 2 . Here t k stands for the approximate value of the k -th zero point. This indicates Proposition 2 given by Riemann is true. Explicitly, for a big enough T which does not accord with a zero point, in the range 0 < t < T the number of zero points on the critical line σ = 1 / 2 obeys the exact formula bellow (understood as taking the integer part):

N ( T ) = T 2 π l o g ( T 2 π e ) + 4 3 − ε ( 1 / 2 , T ) π

with | ε ( 1 / 2 , T ) / π | < 1 / 2 .

It is well-known that there is no zero point outside the critical strip 0 ≤ σ ≤ 1 . For any given finite region bounded by 0 ≤ σ ≤ 1 and 0 < t < T , it can be always checked that, all the zero points of ξ ( σ , t ) are on the critical line σ = 1 / 2 . What lacks is the knowledge for the case with t → ∞ . In 2004, Gourdon  had checked the first 1013 zero points and found the Riemann hypothesis is true. To speculate with our distribution law, its upper bound is up to T = 2.4460 × 10 12 . As reviewed in , another calculation had ever been done by Odlyzko in 1989  in selected intervals show that the Riemann hypothesis holds for over 3 × 108 zeros at height up to T = 2 × 10 20 . Our request on the lower bound T may be much lower than this. The final result is as follows:

Except on the critical line σ = 1 / 2 , the real part U ( σ , t ) and imaginary part V ( σ , t ) of ξ ( σ + i t ) have no other mutual zero point in the region [ 0,1 ] × ( T , ∞ ) ∪ [ 0,1 ] × ( − ∞ , − T ) for big enough T.

This indicates the Riemann hypothesis is true, and it is the moment to draw full stop for this suspending problem.

In addition to the “Riemann hypothesis”, we had also explored another well- known problem named “P versus NP” in , where a brand-new fast algorithm for the “travelling salesman problem” was constructed, and a surprising result “P = NP” was obtained. Those who are interested can read it.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Wang, J.L. (2021) The Proof of Riemann Hypothesis, the Key to the Door Is the Periodicity. Applied Mathematics, 12, 723-765. https://doi.org/10.4236/am.2021.128050

ReferencesNew Scientist (2018) Has the Riemann Hypothesis Been Solved? https://www.newscientist.com/section/news/Titchmarsh, E.C. (1986) The Theory of the Riemann Zeta-Function. 2nd Edition, Edited and with a Preface by Heath-Brown, D.R., The Clarendon Press, Oxford University Press, Oxford.Edwards, H.M. (1974) Riemann’s Zeta Function. Academic Press, Waltham. (Reprinted by Dover Publications, Mineola, 2001).Borwein, P., Choi, S., Rooney, B. and Weirathmueller, A. (2007) The Riemann Hypothesis: A Resource for the Afficionado and Virtuoso Alike. Springer Press, New York. https://doi.org/10.1007/978-0-387-72126-2Lu, C.H. (2016) Ramble on the Riemann Hypothesis. Tsinghua University Press, Beijing, 20-38. (In Chinese)Devlin, K.J. (2003) The Millennium Problems: The Seven Greatest Unsolved Mathematical Puzzles of Our Time. Basic Books, New York, 19-60.Riemann, B. (1998) On the Number of Prime Numbers Less than a Given Quantity. Translated by Wilkins, D.R. http://www.claymath.org/library/historical/riemann/english/1.pdfRiemann, B. (2016) The Collected Works of Bernhard Riemann. Vol. 1, Chinese Edition, Translated by Li, P.L., Higher Education Press, Beijing, 127-135.Bombieri, E. (2000) Problems of the Millennium: The Riemann Hypothesis. http://www.claymath.org/millennium-problems/riemann-hypothesisLune, J., Riele, J.J. and Winter, D.T. (1986) On the Zeros of the Riemann Zeta Function in the Critical Strip, IV. Mathematics of Computation, 46, 667-681. https://doi.org/10.1090/S0025-5718-1986-0829637-3Gourdon, X. (2004) The 1013 First Zeros of the Riemann Zeta Function, and Zeros Computation at Very Large Height. http://numbers.computation.free.fr/Constants/Miscellaneous/zetazeros1e13-1e24.pdfStein, E.M. and Shakarchi, R. (2003) Princeton Lectures in Analysis I: Fourier Analysis, An Introduction. Princeton University Press, Princeton.Stein, E.M. and Shakarchi, R. (2003) Princeton Lectures in Analysis II: Complex Analysis. Princeton University Press, Princeton.Kudryavtseva, A.E., Saidak, F. and Zvengrowski, P. (2005) Riemann and His Zeta Function. Morfismos, 9, 1-48.Selberg A. ,et al. (1942)On the Zeros of the Zeta-Function of Riemann. Norske Vid. Selsk. Forh Trondhjem 15, 59-62.Levinson, N. (1974) More than One-Third of the Zeros of the Riemann Zeta-Function Are on σ =1/2. Advances in Mathematics, 13, 383-436. https://doi.org/10.1016/0001-8708(74)90074-7Conrey, J.B. (1989) More than Two Fifths of the Zeros of the Riemann Zeta Function Are On the Critical Line. Journal für die reine und angewandte Mathematik, No. 399, 1-26. https://doi.org/10.1515/crll.1989.399.1Odlyzko, A. (n.d.) Tables of Zeros of the Riemann Zeta Function. http://www.dtc.umn.edu/~odlyzko/zeta_tables/index.htmlGilbarg, D. and Trudinger, N.S. (2001) Laplace’s Equation. In: Elliptic Partial Differential Equations of Second Order, Vol. 224, Springer-Verlag, Berlin, Heidelberg, 13-30. https://doi.org/10.1007/978-3-642-61798-0_2Wang, J.L. (2018) Fast Algorithm for the Travelling Salesman Problem and the Proof of P = NP. Applied Mathematics, 9, 1351-1359. https://doi.org/10.4236/am.2018.912088Odlyzko, A.M. (1989) Supercomputers and the Riemann Zeta Function, Supercomputing 89: Supercomputing Structures & Computations. Proceedings of 4th International Conference on Supercomputing, International Supercomputing Institute, 348-352.