_{1}

^{*}

The behavior of the zeros in finite Taylor series approximations of the Riemann Xi function (to the zeta function), of modified Bessel functions and of the Gaussian (bell) function is investigated and illustrated in the complex domain by pictures. It can be seen how the zeros in finite approximations approach to the genuine zeros in the transition to higher-order approximation and in case of the Gaussian (bell) function that they go with great uniformity to infinity in the complex plane. A limiting transition from the modified Bessel functions to a Gaussian function is discussed and represented in pictures. In an Appendix a new building stone to a full proof of the Riemann hypothesis using the Second mean-value theorem is presented.

The present paper tries to find out the common ground for the zeros of the Riemann zeta function ζ ( z ) = ζ ( x + i y ) and of the modified Bessel functions I ν ( z ) (or Bessel functions J ν ( z ) of imaginary argument z) for imaginary argument z and, furthermore, for the absence of zeros of the Gaussian Bell function exp ( z 2 ) . For the function now called Riemann zeta function ζ ( z ) which was known already to Euler but was extended by Riemann to the complex plane Riemann expressed the hypothesis that all nontrivial zeros of this function lie on the axis z = 1 2 + i y that means on the axis through x = 1 2 and parallel to the imaginary axis y (Riemann hypothesis) [

The content of this article was not intended as a proof of the Riemann hypothesis but during the work we found a further, as it seems, essential building stone for its proof by the second mean-value approach which is represented in Appendix. The article is merely intended as an illustration to the zeros of a function with a possible representation in an integral form given in Section 2 (Equation (2.8)) with monotonically decreasing functions Ω ( u ) satisfied by the Riemann Xi function and by the modified Bessel functions and explains why the Gauss bell function although it can be represented in such form does not possess zeros. Other kinds of interesting illustrations of the Riemann zeta function (and of other functions) by the Newton flow are given in [

The Xi function Ξ ( z ) to the Riemann zeta function ζ ( z ) is defined by

Ξ ( z ) ≡ ξ ( 1 2 + z ) , (2.1)

where ζ ( z ) is the Riemann Xi function [^{1}

ξ ( z ) ≡ ( z − 1 ) ( z 2 ) ! π − z 2 ζ ( z ) . (2.2)

The Riemann zeta function is basically defined by the following Euler product

ζ ( z ) ≡ ∏ n = 1 ∞ ( 1 − 1 ( p n ) z ) − 1 , ( p 1 = 2, p 2 = 3, p 3 = 5, p 4 = 7, ⋯ ) , (2.3)

where p n is the sequence of prime numbers beginning with p 1 = 2 . The definition of the Riemann Xi function (2.3) is equivalent to the definition by the following (Dirichlet) series for complex variable z = x + i y

ζ ( z ) = ∑ n = 1 ∞ 1 n z , ( x = Re ( z ) > 1 ) , (2.4)

which is convergent for x > 1 and arbitrary y and can be analytically continued to the whole complex z-plane. The function Ξ ( z ) is an entire function which excludes the only singularity of ζ ( z ) at z = 1 and its “trivial” zeros at z = − 2 , − 4 , − 6 , ⋯ .

Next we consider the whole class of modified Bessel functions I ν ( z ) of imaginary argument which is connected with the basic class of Bessel functions J ν ( z ) in the following slightly modified form by ( μ ! ≡ Γ ( μ + 1 ) )

( 2 z ) ν I ν ( z ) = ( 2 i z ) ν J ν ( i z ) ≡ ∑ m = 0 ∞ 1 m ! ( m + ν ) ! ( z 2 ) 2 m , ( ν ∈ ℝ ) . (2.5)

The functions ( 2 z ) ν I ν ( z ) are entire functions which satisfy the differential equation

{ ( z ∂ ∂ z ) 2 + 2 ν z ∂ ∂ z − z 2 } ( 2 z ) ν I ν ( z ) = 0 ⇔ { ( z ∂ ∂ z ) 2 − z 2 − ν 2 } I ν ( z ) = 0. (2.6)

In comparison to I ν ( z ) the functions ( 2 z ) ν I ν ( z ) exclude the zeros or infinities of the first ones at z = 0 but the other zeros remain the same for bothfunctions.

Finally, we consider the Gaussian functions exp ( a 2 z 2 4 ) with parameter a 2 > 0 which can be represented by the following integral representation (continued from the imaginary axis y to the whole complex z-plane)

exp ( a 2 z 2 4 ) = 2 π a 2 ∫ 0 + ∞ d u exp ( − u 2 a 2 ) ch ( u z ) , ( a 2 > 0 ) , (2.7)

which become Gaussian bell functions for imaginary argument z = i y . Clearly, the functions exp ( a 2 z 2 4 ) do not possess zeros on the imaginary axis and zeros at all.

The three mentioned types of functions written as Ξ ( z ) have in common that they are symmetrical functions in z and that they possess a representation by an integral of the type

Ξ ( z ) = ∫ 0 + ∞ d u Ω ( u ) ch ( u z ) ) = Ξ ( − z ) = ( Ξ ( z * ) ) * , Ω ( u ) = Ω ( − u ) , (2.8)

with monotonically decreasing functions Ω ( u ) for 0 ≤ u ≤ ∞ that means

0 ≤ u 1 ≤ u 2 ⇒ Ω ( u 1 ) ≥ Ω ( u 2 ) , Ω ( u → + ∞ ) = 0. (2.9)

The Taylor series of Ξ ( z ) can be written in the form

Ξ ( z ) = ∑ m = 0 ∞ Ω 2 m z 2 m , Ω n ≡ 1 n ! ∫ − ∞ + ∞ d u 1 2 ( Ω ( u ) + Ω ( − u ) ) u n

Ω 2 m = 1 ( 2 m ) ! ∫ 0 + ∞ d u Ω ( u ) u 2 m , Ω 2 m + 1 = 0, ( m = 0,1,2, ⋯ ) , (2.10)

where Ω n , ( n = 0 , 1 , 2 , ⋯ ) are defined as the moments of the symmetrical function Ω ( u ) with respect to the reference point u 0 = 0 . A consequence of the definitions in (2.8) and (2.10) is

Ξ ( 0 ) = ∫ 0 + ∞ d u Ω ( u ) ≡ Ω 0 . (2.11)

The odd moments Ω 2 m + 1 of the function Ω ( u ) in the definition (2.10) vanish. The function Ξ ( z ) at z = 0 is equal to the zeroth moment Ω 0 of the function Ω ( u ) and is independent of the chosen reference point. In the following the moments of the function Ω ( u ) play an important role.

That Ω ( u ) is a symmetrical function in u is, in principle, not necessary since the integration over u in (2.8) is restricted by u ≥ 0 but the symmetry permits to extend the integral over negative values of u using it in the form

Ξ ( z ) = 1 2 ∫ − ∞ + ∞ d u Ω ( u ) e u z = 1 2 ∫ − ∞ + ∞ d u Ω ( u ) ch ( u z ) + 1 2 ∫ − ∞ + ∞ d u Ω ( u ) sh ( u z ) ︸ = 0 , (2.12)

which for imaginary z = i y is a Fourier transformation of Ω ( u ) with the inversion

Ω ( u ) = 1 π ∫ − ∞ + ∞ d y Ξ ( i y ) e − i u y , (2.13)

if the integral exists in some sense (e.g. weak convergence). We consider this now more explicitly.

The explicit representation of the Xi function Ξ ( z ) to the Riemann Xi function ξ ( z ) = Ξ ( z − 1 2 ) in this form (2.8) together with (2.9) is

Ξ ( z ) = ∫ 0 + ∞ d u Ω ( u ) ch ( u z ) ,

Ω ( u ) ≡ 4 exp ( u 2 ) ∑ n = 1 ∞ π n 2 e 2 u ( 2 π n 2 e 2 u − 3 ) exp ( − π n 2 e 2 u ) = Ω ( − u ) , (2.14)

with the special values

Ξ ( 0 ) = ∫ 0 + ∞ d u Ω ( u ) ≈ 0.4971207782, Ξ ( ± 1 2 ) = 1 2 ,

Ω ( 0 ) ≈ 1.7867876019. (2.15)

This was derived in detail in [

The function Ω ( u ) is monotonically decreasing from u = 0 on and decreases rapidly for u → ± ∞ . The function Ξ ( z ) increases rapidly for real z = x → ± ∞ and decreases rapidly for imaginary z = i y → ± i ∞ . The symmetry

of Ω ( u ) is not easily to see from (2.14) and it was a genuine surprise for us to meet such a kind of a symmetrical function (see discussion in [

In case of the modified Bessel functions of imaginary argument the following basic integral representations are known which for the functions proportional to ( 2 z ) ν I ν ( z ) may be written as follows (taking into account ( 1 2 ) ! = 1 2 ( − 1 2 ) ! = π 2 )

Ξ ν ( z ) ≡ Ω ν ( 0 ) ∫ 0 1 d u ( 1 − u 2 ) ν − 1 2 ch ( u z ) = Ω ν ( 0 ) ( 1 2 ) ! ( ν − 1 2 ) ! ( 2 z ) ν I ν ( z ) , ( ν > − 1 2 ) , (2.16)

or expressed by the Confluent Hypergeometric function 1 F 1 ( a ; c ; z )

Ξ ν ( z ) = Ω ν ( 0 ) ( 1 2 ) ! ( ν − 1 2 ) ! ν ! e ∓ z 1 F 1 ( ν + 1 2 ; 2 ν + 1 ; ± 2 z ) . (2.17)

The Taylor series expansion is

Ξ ν ( z ) = Ω ν ( 0 ) ( 1 2 ) ! ( ν − 1 2 ) ! ν ! { 1 + ∑ m = 1 ∞ ν ! m ! ( ν + m ) ! ( z 2 ) 2 m } . (2.18)

The functions Ξ ν ( z ) possess the principal form

Ξ ν ( z ) = ∫ 0 + ∞ d u Ω ν ( u ) ch ( u z ) , ( ν > − 1 2 ) , (2.19)

with the following explicit expressions for Ω ν (u)

Ω ν ( u ) = Ω ν ( 0 ) ( 1 − u 2 ) ν − 1 2 θ ( 1 − u 2 ) = Ω ν ( − u ) , (2.20)

where θ ( x ) is Heaviside’s step function defined by

θ ( x ) ≡ { 1 , ( x > 0 ) , 0 , ( x < 0 ) . (2.21)

It restricts the upper limit of integration in (2.19) according to the choice in (2.20) to u = 1 . The functions Ω ν ( u ) are equivalent to functions of u 2 only and are in this sense symmetrical functions of u and for real u 2 < 1 we have to choose the real value of ( 1 − u 2 ) ν − 1 2 in case of non-integer ν − 1 2 .

The first six cases of the function Ξ ν ( z ) with integer or semi-integer index ν are

Ξ 0 ( z ) = Ω 0 ( 0 ) π I 0 ( z ) 2 , < /mtext> Ξ 1 2 ( z ) = Ω 1 2 ( 0 ) sh ( z ) z , Ξ 1 ( z ) = Ω 1 ( 0 ) π I 1 ( z ) 2 z , < /mtext> Ξ 3 2 ( z ) = Ω 3 2 ( 0 ) 2 ( z ch ( z ) − sh ( z ) ) z 3 , Ξ 2 ( z ) = Ω 2 ( 0 ) 3 π I 2 ( z ) 2 z 2 , Ξ 5 2 ( z ) = Ω 5 2 ( 0 ) 8 ( ( 3 + z 2 ) sh ( z ) − 3 z ch ( z ) ) z 5 . (2.22)

The function Ξ − 1 2 ( z ) does not exist since the corresponding integral in (2.16) is divergent or ( ν − 1 2 ) ! is not finite for ν = − 1 2 but lim ν → − 1 2 Ξ − 1 2 ( z ) ( ν − 1 2 ) ! = ch ( z ) 2 exists. Furthermore, we have introduced in (2.16)

amplitudes Ω ν ( 0 ) and we will soon see that it is favorable to choose them for our purpose as constant for the whole class of functions Ξ ν ( z ) . We illustrate the functions Ω ν ( u ) and the corresponding functions Ξ ν ( z ) with a certain important modification in next Section 3.

It is well known that the Bessel functions J ν ( z ) possess zeros only on the real axis x and, correspondingly, I ν ( z ) only on the imaginary axis y. For this there exists a direct proof using their differential equations and which is similar to the derivation of duality and orthogonality relations (e.g. [

The considered Omega functions Ω ( u ) for the Riemann Xi function and Ω ν ( u ) for the modified Bessel functions and also for the Gaussian bell function have in common that they are monotonically decreasing functions for u ≥ 0 up to Ω ( + ∞ ) = 0 .

We now calculate the even moments Ω ν ,2 m of the functions Ω ν ( u ) in (2.20) which lead to well-known integrals

Ω ν ,2 m ≡ Ω ν ( 0 ) ( 2 m ) ! ∫ 0 + ∞ d u Ω ν ( u ) u 2 m = 1 ( 2 m ) ! ∫ 0 1 d u ( 1 − u 2 ) ν − 1 2 u 2 m = Ω ν ( 0 ) ( ν − 1 2 ) ! ( m − 1 2 ) ! 2 ( 2 m ) ! ( ν + m ) ! = Ω ν ( 0 ) ( 1 2 ) ! ( ν − 1 2 ) ! m ! ( ν + m ) ! 2 2 m , ( ( 1 2 ) ! = π 2 ) , (3.1)

and in the special case m = 0 of zeroth moments

Ω ν , 0 = Ω ν ( 0 ) ( 1 2 ) ! ( ν − 1 2 ) ! ν ! . (3.2)

These are the areas under the curves Ω ν ( u ) on the positive u-axes. The following considerations show that it is favorable to make Ω ν ,0 equal independently on parameter ν and to keep Ω ν ( 0 ) constant and we choose

Ω ν ( 0 ) = 1. (3.3)

To keep in addition also Ω ν ,0 constant we have now the only possibility to introduce a stretch factor to the variable u and we make the transformation (now already with choice (3.3))

Ω ν ( u ) → Ω ¯ ν ( u ) = ( 1 − u 2 u ν ,0 2 ) ν − 1 2 θ ( 1 − u 2 u ν ,0 2 ) = Ω ν ( u u ν ,0 ) , (3.4)

with the definition of u ν ,0

u ν , 0 ≡ ν ! ( 1 2 ) ! ( ν − 1 2 ) ! u 1 2 , 0 = u 1 2 , 0 Ω ν , 0 , u 1 2 , 0 = 1 , (3.5)

from which follows for the zeroth moments of Ω ¯ ν (u)

Ω ¯ ν ,0 ≡ ∫ 0 + ∞ d u Ω ν ( u ) = ∫ 0 u ν ,0 d u ( 1 − u 2 u ν ,0 2 ) ν − 1 2 = 1. (3.6)

The functions Ω ¯ ν ( u ) are graphically represented in

values of ν up to ν = 9 2 .

The corresponding transformation of the functions Ξ ν ( z ) using (2.16) and (3.3) are

Ξ ν ( z ) → Ξ ¯ ν ( z ) = u ν ,0 ∫ 0 u ν ,0 d u u ν ,0 ( 1 − u 2 u ν ,0 2 ) ν − 1 2 ch ( u u ν ,0 u ν ,0 z ) = u ν ,0 Ξ ν ( u ν ,0 z ) = ν ! ( 2 u ν ,0 z ) ν I ν ( u ν ,0 z ) . (3.7)

Their Taylor series are

Ξ ¯ ν ( z ) = ∑ m = 0 ∞ ν ! m ! ( ν + m ) ! ( u ν ,0 z 2 ) 2 m = 1 + ∑ m = 1 ∞ ν ! m ! ( ν + m ) ! ( ν ! ( 1 2 ) ! ( ν − 1 2 ) ! z 2 ) 2 m = Ξ ¯ ν ν ( − z ) , (3.8)

or with substitution z = ∓ i y (see also (2.5))

Ξ ¯ ν ( ∓ i y ) = ν ! ( 2 u ν ,0 y ) ν J ν ( u ν ,0 y ) = 1 + ∑ m = 1 ∞ ( − 1 ) m ν ! m ! ( ν + m ) ! ( ν ! ( 1 2 ) ! ( ν − 1 2 ) ! y 2 ) 2 m . (3.9)

The first three functions Ξ ¯ ν ( z ) with integer and with semi-integer ν are explicitly

Ξ ¯ 0 ( z ) = I 0 ( 2 z π ) , Ξ ¯ 1 2 ( z ) = sh ( z ) z , Ξ ¯ 1 ( z ) = π 2 z I 1 ( 4 z π ) , Ξ ¯ 3 2 ( z ) = 4 ( 3 z ch ( 3 z 2 ) − 2 sh ( 3 z 2 ) ) 9 z 3 , Ξ ¯ 2 ( z ) = 9 π 2 32 z 2 I 2 ( 16 z 3 π ) , Ξ ¯ 5 2 ( z ) = 512 ( ( 64 + 75 z 2 ) sh ( 15 z 8 ) − 120 z ch ( 15 z 8 ) ) 16875 z 5 . (3.10)

The functions Ξ ¯ ν ( z ) on the imaginary axes z = i y for integer and semi-integer index ν up to ν = 9 2 are illustrated in

That the mentioned approach for ν → ∞ to a Gaussian bell function is really true we establish exactly and determine these limits.

We now show that the new functions Ξ ¯ ν ( z ) in (3.7) for ν → ∞ go to a Gaussian function that becomes a Gaussian bell function for imaginary z = i y . As auxiliary formulae for ν ! ( ν − 1 2 ) ! in case of ν ≫ 1 follows

ν ! ( ν − 1 2 ) ! ︸ ≡ φ ( ν ) ( ν − 1 2 ) ! ( ν − 1 ) ! ︸ ≡ φ ( ν − 1 2 ) = ν ≈ ν + 1 4 ν − 1 4 ⇒ ν ! ( ν − 1 2 ) ! ≈ ν + 1 4 → ν ,

ν ! ( 1 2 ) ! ( ν − 1 2 ) ! ≈ 2 π ν + 1 4 ≈ 2 π ν , ( ( 1 2 ) ! = π 2 ) , (3.11)

and analogously for ν ≫ 1

ν ! ( ν + m ) ! = 1 ( ν + 1 ) ( ν + 2 ) ⋯ ( ν + m ) = 1 ν m ( 1 + 1 ν ) ( 1 + 2 ν ) ⋯ ( 1 + m ν ) ≈ ν − m . (3.12)

Applying these approximations we find from (3.8)

lim ν → ∞ Ξ ¯ ν ( z ) = lim ν → ∞ ∑ m = 0 ∞ ν ! m ! ( ν + m ) ! ( ν ! ( 1 2 ) ! ( ν − 1 2 ) ! z 2 ) 2 m = lim ν → ∞ ∑ m = 0 ∞ ν − m m ! ( ν 2 π z 2 ) 2 m = ∑ m = 0 ∞ 1 m ! ( z 2 π ) m = exp ( z 2 π ) . (3.13)

This means that in the limiting transition ν → ∞ the functions Ξ ¯ ν ( z ) in (3.7) approach to a Gaussian function according to

lim ν → ∞ ν ! ( 2 ( 1 2 ) ! ( ν − 1 2 ) ! ν ! z ) ν I ν ( ν ! ( 1 2 ) ! ( ν − 1 2 ) ! z 2 ) = exp ( z 2 π ) , (3.14)

or by substitution z = i y with, in general, complex variable y

lim ν → ∞ ν ! ( 2 ( 1 2 ) ! ( ν − 1 2 ) ! ν ! y ) ν J ν ( ν ! ( 1 2 ) ! ( ν − 1 2 ) ! y 2 ) = exp ( − y 2 π ) , (3.15)

where we use the identity 2 ( 1 2 ) ! = ( − 1 2 ) ! = π . From the corresponding limiting transition lim ν → ∞ Ω ¯ ν ( u ) using the definition in (3.4) we find

lim ν → ∞ ( 1 − ( ( 1 2 ) ! ( ν − 1 2 ) ! ν ! ) 2 u 2 ) ν − 1 2 θ ( 1 − ( ( 1 2 ) ! ( ν − 1 2 ) ! ν ! ) 2 u 2 ) = exp ( − π u 2 4 ) . (3.16)

We may check the transition from Ω ( u ) = exp ( − π u 2 4 ) to Ξ ( z ) = exp ( z 2 π ) via (2.8) using the auxiliary formula (2.7) with a 2 = 4 π .

Furthermore, we find using the approximation (3.11) that the factor between u 1 2 ,0 and u ν ,0 in (3.5) approaches for ν ≫ 1 according to

u ν , 0 u 1 2 , 0 = ν ! ( 1 2 ) ! ( ν − 1 2 ) ! ≈ 2 π ν + 1 4 , ( ν ≫ 1 ) , (3.17)

with a high precision and monotonically increasing without a finite limit.

We now show that for high increasing ν ≫ 1 the first positive roots of Ξ ¯ ν ν ( i y ) and thus also the higher roots increase. For Ξ ¯ ν ( i y ) this is clear since it is proportional to J ν ( y ) and it is well known that their roots y ν , s , ( s = 1 , 2 , ⋯ ) are situated approximately at

y ν , s = ( s + ν 2 − 1 4 ) π − ⋯ . (3.18)

However in Ξ ¯ ν ( i y ) the argument in the Bessel functions is stretched by the factor u ν ,0 and therefore on the imaginary axis we have the functions J ν ( u ν ,0 y ) that diminishes the values for the roots (3.18) by the factor u ν ,0 − 1 and with the approximation (3.17) the zeros of Ξ ¯ ν ( i y ) are now situated approximately at

y ν , s u ν , 0 ≈ ( s + ν 2 − 1 4 ) π 3 2 4 ν + 1 . (3.19)

Thus the first roots ( s = 1 ) of Ξ ¯ ν ( i y ) go with ν → ∞ proportional to ν also to infinity. This explains more in detail why the first zeros ( s = 1 ) of the curves for Ξ ¯ ν ( i y ) in

We remind that the starting point for the derivation of the approximations was making equal the zeroth moments Ω ¯ ν ,0 by stretching the argument of the functions Ω ν ( u ) to new functions Ω ¯ ν ( u ) leaving constant the amplitudes Ω ν ( 0 ) . The zeros of Ω ¯ ν ( u ) for ν → ∞ move in this process also to infinity but “very slowly”. This is a good illustration why the Gaussian function exp ( z 2 π ) which on the imaginary axis y becomes a Gaussian bell function exp ( − y 2 π ) does not possess zeros at all although its Omega functio Ω ( u ) = exp ( − π u 2 4 ) is monotonically decreasing. The first zero and in this way all other zeros are moved by the limiting transition to infinity although very slowly. The same is the case with the discontinuities of derivatives in the functions Ω ¯ ν ( u ) .

In this Section we consider the Taylor series expansions of the Xi function Ξ ( z ) to the Riemann Xi function ξ ( s ) in powers of z defined in (2.14) according to

Ξ ( z ) = ∫ 0 + ∞ d u Ω ( u ) ch ( u z ) = ∑ m = 0 ∞ Ω 2 m z 2 m = ∑ m = 0 ∞ Ξ ( 2 m ) ( 0 ) ( 2 m ) ! z 2 m ,

Ω 2 m ≡ 1 ( 2 m ) ! ∫ 0 + ∞ d u Ω ( u ) u 2 m = Ξ ( 2 m ) ( 0 ) ( 2 m ) ! . (4.1)

The coefficients in this Taylor series are the even moments Ω 2 m of the function Ω ( u ) as defined.

We truncate the Taylor series of Ξ ( z ) at upper summation numbers m = M

Ξ 2 M ( z ) = ∑ m = 0 M Ω 2 m z 2 m , (4.2)

and calculate all zeros of these approximations up to a certain maximal M and make graphical representations of their zeros. As mentioned the coefficients of the series are determined by the even moments Ω 2 m of the function Ω ( u ) where the odd moments vanish automatically if we determine Ω ( u ) as symmetrical functions according to (2.8) and (2.10). In the following Sections we discuss the same also for modified Bessel functions I ν ( z ) (in particular for ν = 1 2 ) and for the limiting transition to a Gaussian function and compare this with the zeros for the Xi function (2.14) to the Riemann zeta function. In each approximation Ξ 2 M ( z ) of the Taylor series we find 2M complex zeros. It was very interesting to see which difference appears between the zeros on the imaginary axis and the bulk of the other zeros with increasing number M. This is best seen in graphical representations of the zeros.

For the function Ω ( u ) related to the Riemann hypothesis which is explicitly given in (2.14) we calculated this series numerically up to 2 M = 100 with sufficiently high precision and obtained the following numerical coefficients (we write here explicitly down only a few of the obtained coefficients)^{2}

Ξ ( z ) = 0.497120778188314 + 1.14859721575727 × 10 − 2 z 2 + 1.23452018070318 × 10 − 4 z 4 + 8.32355481385527 × 10 − 7 z 6 + 3.99222655134413 × 10 − 9 z 8 + ⋯ + 1.05272334981972 × 10 − 68 z 48 + 4.68273888631736 × 10 − 72 z 50 + 1.97190495510121 × 10 − 75 z 52

+ 7.87602897895070 × 10 − 79 z 54 + 2.98910175366638 × 10 − 82 z 56 + 1.07971892491186 × 10 − 85 z 58 + 3.71787541158066 × 10 − 89 z 60 + 1.22215971011433 × 10 − 92 z 62 + ⋯ + 3.09393689867254 × 10 − 158 z 98 + 5.11093592293365 × 10 − 162 z 100 + ⋯ (4.3)

It is interesting to mention that the coefficients ( 2 m ) ! Ω 2 m , ( m = 1 , 2 , ⋯ ) possess an absolute minimum for m = 21 with ( 42 ) ! Ω 42 ( 0 ) = 1.1837577 × 10 − 7 and, apparently, are (“slowly”) monotonically increasing for m > 21 = 42 .

The first term in (4.3) can also be written (see [

0.497120778188313661 ≈ Ξ ( 0 ) = 1 2 − ∫ 1 + ∞ d t 2 t ∑ n = 1 ∞ exp ( − π n 2 t 2 ) = 1 2 − ∫ 1 + ∞ d s ∑ n = 1 ∞ exp ( − π n 2 s 4 ) ≈ 0.5 − 0.002879221811686339 , (4.4)

corresponding to the splitting of terms in the representation (2.14) of Ξ ( z ) after the substitution e u = t = s 2 .

In the truncation of the series in powers of z with the highest term proportional to z 2 M we found by numerical solution of the corresponding algebraic equations of degree 2M the 2M complex solutions from which the following pairs of zeros lie on the imaginary axis z = i y ^{3}:

2M | zeros on imaginary axis | 2M | zeros on imaginary axis | |
---|---|---|---|---|

2 | ±i6.57881 | 42 | ±i14.1347 | (4.5) |

4 | - | 44 | ±i14.1347, ±i20.2528 | |

6 | ±i8.48024 | 46 | ±i14.1347 | |

8 | - | 48 | ±i14.1347, ±i20.8469 | |

10 | ±i10.0594 | 50 | ±i14.1347, ±i21.0934, ±i21.8540 | |

12 | - | 52 | ±i14.1347, ±i21.0123 | |

14 | ±i11.4267 | 54 | ±i14.1347, ±i21.0238, ±i23.0702 | |

16 | - | 56 | ±i14.1347, ±i21.0218 | |

18 | ±i12.6077 | 58 | ±i14.1347, ±i21.0221, ±i24.0109 | |

20 | - | 60 | ±i14.1347, ±i21.0220 | |

22 | ±i13.5513 | 62 | ±i14.1347, ±i21.0220, ±i24.7095 | |

24 | - | 64 | ±i14.1347, ±i21.0220 | |

26 | ±i14.0560 | 66 | ±i14.1347, ±i21.0220, ±i24.9857 | |

28 | ±i14.1530, ±i15.8936 | 68 | ±i14.1347, ±i21.0220, ±i25.0147, ±i26.7482 | |

30 | ±i14.1319 | 70 | ±i14.1347, ±i21.0220, ±i25.0093 | |

32 | ±i14.1352, ±i17.2309 | 72 | ±i14.1347, ±i21.0220, ±i25.0102, ±i27.8069 | |

34 | ±i14.1347 | 74 | ±i14.1347, ±i21.0220, ±i25.0101 | |

36 | ±i14.1347, ±i18.3765 | 76 | ±i14.1347, ±i21.0220, ±i25.0101, ±i28.8689 | |

38 | ±i14.1347 | 78 | ±i14.1347, ±i21.0220, ±i25.0101 | |

40 | ±i14.1347, ±i19.3895 | 80 | ±i14.1347, ±i21.0220, ±i25.0101, ±i30.5405 |

In the following we give graphical illustrations of all zeros of the Taylor approximations in the complex z-plane where all zeros up to a certain order 2M are taken into account and where one may see how the zeros change from an order to a higher order.

We explain first how the following figures are made. We take a certain order 2M of the Xi function Ξ ( z ) given by the truncated Taylor series (4.2) and determine numerically all of its zeros and represents them by points in the complex ( z = x + i y ) -plane where we choose the same scale on the x- and y-axis. Two variants are made, first the representation by isolated points and second the representation by connected neighbored points. The obtained partial pictures are a little different for odd and even M that is represented in

Then we calculate and represent all zeros of the Taylor approximations of Ξ 2 M ( z ) (i.e. of Ξ 2 M ( z ) = 0 ) up to a certain maximal M and represent the zeros in described way by isolated points and by connected neighbored points in each of the approximations up to the maximal M. This is made in

only for the first 40 Taylor approximations where this is clearer to see. These pictures show that the zeros on the imaginary axis stabilize from order to higher orders at the genuine zeros of the Riemann zeta function on the y-axis and separate themselves from the main bulk of zeros in a considered order which do not lie on the imaginary axis. That this remains in this way for M → ∞ is, clearly, only a conjecture but in Section 8 we try to understand this by some analytic approximations.

In

It is interesting to compare the functions Ω ( u ) and Ξ ( z ) to the Riemann zeta function with a corresponding Gaussian functions Ω G ( u ) and Ξ G ( z ) where the two parameters of the Gaussian function Ω G ( u ) which are the amplitude and the stretching of the parameter u are chosen in the way that the first two terms of the Taylor series approximation are equal. For the function Ω (u)

Ω ( u ) = Ω ( 0 ) + Ω ( 2 ) ( 0 ) 2 ! u 2 + Ω ( 4 ) ( 0 ) 4 ! u 4 + … ≈ 1.78678760 − 16.7305008 u 2 + 67.6802 u 4 − ⋯ , (4.6)

this means for the Gaussian function Ω G (u)

Ω G ( u ) = Ω ( 0 ) exp ( Ω ( 2 ) ( 0 ) 2 Ω ( 0 ) u 2 ) = Ω ( 0 ) + Ω ( 2 ) ( 0 ) 2 u 2 + ( Ω ( 2 ) ( 0 ) ) 2 2 ! 2 2 Ω ( 0 ) u 4 + ⋯ ≈ 1.78678760 ⋅ exp ( − 9.363452465 u 2 ) ≈ 1.78678760 − 16.7305008 u 2 + 78.3276 u 4 − ⋯ . (4.7)

One may check that generally

Ω G ( u ) ≥ Ω ( u ) > 0. (4.8)

For the moments of these functions result the inequalities ( m = 0 , 1 , 2 , ⋯ )

Ω G ,2 m ≡ 1 ( 2 m ) ! ∫ 0 ∞ d u Ω G ( u ) u 2 m > 1 ( 2 m ) ! ∫ 0 ∞ d u Ω ( u ) u 2 m = Ω 2 m > 0,

Ξ G ( 2 m ) ( 0 ) ≡ ( 2 m ) ! Ω G ,2 m > ( 2 m ) ! Ω 2 m ≡ Ξ ( 2 m ) ( 0 ) . (4.9)

The resulting function Ξ G ( z ) is a Gaussian function of the complex variable z = x + i y

Ξ G ( z ) = Ω ( 0 ) − π Ω ( 0 ) 2 Ω ( 2 ) ( 0 ) exp ( − Ω ( 0 ) 2 Ω ( 2 ) ( 0 ) z 2 ) ≈ 0.517487503 ⋅ exp ( 0.0266995535 z 2 ) . (4.10)

It possesses another amplitude in comparison to Ξ ( z ) (i.e. to 0.497120778). Clearly, as a Gaussian function it does not possess zeros on the y-axis and zeros at all. The Ω ( u ) function to the Riemann Xi function and the considered modified Bessel functions possess the common property that they vanish in infinity more rapidly (or are even finite) than the Omega function to a Gauss function. In principle, this does not exclude Omega functions which vanish less rapidly. For example, a function Ω ( u ) = e α u , ( u ≥ 0, Re ( α ) > 0 ) provides Ξ ( z ) = α α 2 − z 2 without zeros of Ξ ( z ) at all in finite regions of the complex plane but with poles that regarding the zeros is the same as for a Gauss function (see Section 8 where this is explained by motion to infinity from finite approximations). However, by comparison with Gauss functions we may get inequalities for the moments of the considered Omega functions.

In this Section we now come back to the modified Bessel functions Ξ ν ( z ) in (3.7) to the Omega functions Ω ν ( u ) in (3.4). From this class of functions which in the limit ν → ∞ go to a Gaussian function we choose the function sh ( z ) z

Ξ 1 2 ( z ) = Ξ ¯ 1 2 ( z ) = sh ( z ) z , (5.1)

and give graphical representations of the zeros for its Taylor series approximations

Ξ 1 2 ,2 M ( z ) = Ξ ¯ 1 2 ,2 M ( z ) = ∑ m = 0 M z 2 m ( 2 m + 1 ) ! = ∑ m = 0 M ( 1 2 ) ! m ! ( m + 1 2 ) ! ( z 2 ) 2 m . (5.2)

They are represented in

z 0 = i n π , ( n = ± 1 , ± 2 , ⋯ ) ⇒ sh ( z 0 ) z 0 = 0. (5.3)

We may expect that for other low values of index ν ≥ 0 in the functions Ξ ν ( z ) we get similar illustrations with small topological distortions of the pictures for Ξ 1 2 ( z ) = sh ( z ) z .

The function sh ( z ) z possesses a peculiar importance since for monotonically decreasing functions Ω ( u ) (for 0 ≤ u < + ∞ with Ω ( + ∞ ) = 0 ) we may apply the second mean-value theorem (Gauss-Bonnet theorem; see, e.g. Courant [

Ξ ( z ) = ∫ 0 + ∞ d u Ω ( u ) ch ( u z ) = Ω ( 0 ) sh ( w 0 ( z ) z ) z . (5.4)

Herein, w 0 ( z ) is a mean value function from which we may assume that it is an analytic function of z since the integral on the left-hand side depends analytically on z. In next Section we consider shortly the problem of zeros of the function Ξ ( z ) of the form (5.4).

The problem of zeros of functions Ξ ( z ) of the form (5.4) leads essentially to the problem of zeros of sh ( w 0 ( z ) z ) with exclusion of the zero at z = 0 . If we separate the real and imaginary part of w 0 ( z ) according to

w 0 ( z ) = u 0 ( x , y ) + i v 0 ( x , y ) , z = x + i y , (6.1)

and if we then separate the real and imaginary part of sh ( w 0 ( z ) z ) the function (5.4) can be written in the form (the separation of real and imaginary part in the first factor is uninteresting since a zero at z = 0 is excluded)

Ξ ( x + i y ) = Ω ( 0 ) x + i y sh ( u 0 ( x , y ) x − v 0 ( x , y ) y + i ( u 0 ( x , y ) y + v 0 ( x , y ) x ) ) = Ω ( 0 ) x + i y { sh ( u 0 ( x , y ) x − v 0 ( x , y ) y ) cos ( u 0 ( x , y ) y + v 0 ( x , y ) x ) + i ch ( u 0 ( x , y ) x − v 0 ( x , y ) y ) sin ( u 0 ( x , y ) y + v 0 ( x , y ) x ) } . (6.2)

From both forms of the right-hand side in (6.2) we find that for zeros of Ξ ( x + i y ) the following two conditions [

u 0 ( x , y ) y + v 0 ( x , y ) x = n π , n = ± 1, ± 2, ⋯ ,

u 0 ( x , y ) x − v 0 ( x , y ) y = 0 , (6.3)

have to be satisfied at the same time and this is necessary and sufficient.

We now consider the special case of Ξ ( z ) on the imaginary axis y and find from (6.2)

Ξ ( i y ) = Ω ( 0 ) y sin ( ( u 0 ( 0, y ) + i v 0 ( 0, y ) ) y ) = Ω ( 0 ) sin ( w 0 ( y ) y ) y . (6.4)

Since due to symmetries (2.8) the function Ξ ( i y ) has to be a real-valued function that for y ≠ 0 in (6.4) is only possible if v 0 ( 0, y ) vanishes we find on the imaginary axis

v 0 ( 0, y ) = 0 ⇒ Ξ ( i y ) = Ω ( 0 ) y sin ( u 0 ( 0, y ) y ) = ∫ 0 + ∞ d u Ω ( u ) cos ( u y ) . (6.5)

This follows also immediately by application of the second mean-value theorem to the integral for Ξ ( i y ) where the mean value can only take on real values here in dependence on y as parameter. Thus for zeros on the imaginary axis ( x = 0 ) the two conditions (6.3) reduce to only one condition

u 0 ( 0 , y ) y = n π , n = ± 1 , ± 2 , ⋯ . (6.6)

The condition in the second line of (6.3) is then identically satisfied. If one knows the function u 0 ( 0 , y ) one may determine the zeros on the imaginary axis analytically or, at least, numerically but for many problems including the considered one it is not necessary to know these zeros explicitly.

We return to the general case of general values x on the real axis with the conditions (6.3) for zeros. From (2.8) follows that Ξ ( z ) should be an analytic function for all z for which the integral exists that means to satisfy the Cauchy-Riemann equations. For a general analytic function w ( z ) = w ( x + i y ) = u ( x , y ) + i v ( x , y ) one can derive by integration of the Cauchy-Riemann equations the following relations in case of v ( 0, y ) = 0 [

u ( x , y ) = cos ( x ∂ ∂ y ) u ( 0, y ) , v ( x , y ) = − sin ( x ∂ ∂ y ) u ( 0, y ) . (6.7)

Now come into play the following operator identities (operator identities are such identities which can be applied to arbitrary functions to provide function identities) [

cos ( x ∂ ∂ y ) y = y cos ( x ∂ ∂ y ) − x sin ( x ∂ ∂ y ) ,

sin ( x ∂ ∂ y ) y = x cos ( x ∂ ∂ y ) + y sin ( x ∂ ∂ y ) . (6.8)

More general identities of such kind can be derived by representing the Cosine and Sine functions by Exponential functions and using that exp ( ± i x ∂ ∂ y ) applied to analytic functions f ( y ) displace the argument of these functions to f ( y ± i x ) that is discussed in [

n π = { y cos ( x ∂ ∂ y ) − x sin ( x ∂ ∂ y ) } u 0 ( 0, u ) = cos ( x ∂ ∂ y ) y u 0 ( 0, y ) ,

0 = { x cos ( x ∂ ∂ y ) + y sin ( x ∂ ∂ y ) } u 0 ( 0 , y ) = sin ( x ∂ ∂ y ) y u 0 ( 0 , y ) . (6.9)

If we now apply the operator cos ( x ∂ ∂ y ) to the first of these conditions and the operator sin ( x ∂ ∂ y ) to the second one we obtain

cos 2 ( x ∂ ∂ y ) y u 0 ( 0 , y ) = n π , sin 2 ( x ∂ ∂ y ) y u 0 ( 0 , y ) = 0. (6.10)

By addition of both equalities using the operator identity

cos 2 ( x ∂ ∂ y ) + sin 2 ( x ∂ ∂ y ) = 1 (6.11)

follows from (6.10)

y u 0 ( 0, y ) = n π , ( n = ± 1, ± 2, ⋯ ) . (6.12)

This consequence for zeros results from both conditions (6.3) and provides a necessary condition for the imaginary values y of all zeros for which their imaginary value has to agree with one of the solutions for zeros on the imaginary axis. Thus we have to take all solutions y for zeros on the imaginary axis which we denote now by y 0 and which satisfy the conditions

y 0 u 0 ( 0 , y 0 ) = n π ⇔ sin ( y 0 u 0 ( 0 , y 0 ) ) y 0 = 0 , (6.13)

and have to look for solutions x at the same time of both equations

u 0 ( x , y 0 ) y 0 + v ( x , y 0 ) x = n π ,

u 0 ( x , y 0 ) x − v ( x , y 0 ) y 0 = 0. (6.14)

If one finds such solutions x = x 0 of one of the conditions then it is almost in all cases not a solution of the other condition and, therefore, x 0 ± i y 0 is then not a zero. Only in very few cases of y 0 depending on the function Ω ( u ) in (2.8) one would expect to find values x = x 0 ≠ 0 as solution of both conditions (6.14). This is a very strong restriction to the satisfaction of both conditions (6.14). It seems that the set of functions Ω ( u ) with such a property can be only a set of measure zero within the set of all possible functions Ω ( u ) in (2.8). One such case is when Ω ( u ) is stepwise constant and monotonically decreasing with equal lengths u 0 of the steps and this is the only case (Appendix A). The corresponding functions Ξ ( z ) are then superpositions of functions sh ( m u 0 z ) z , ( m = 1,2, ⋯ ) with different amplitudes. This was not correctly discussed in [^{4} and we recognized it [

We now consider the Gauss function Ξ G ( z ) = exp ( z 2 ) and calculate all zeros of its low-order Taylor series approximations Ξ G ,2 M (z)

Ξ G ,2 M ( z ) = ∑ m = 0 M z 2 m m ! , ( M = 1 , 2 , ⋯ ) . (7.1)

In

essential differences in the pictures for the Gauss function to the other discussed functions with genuine zeros when we increase the order of approximation. It becomes obvious that in the limiting transition to the Gaussian function exp ( z 2 ) all zeros of the finite-order Taylor approximations go to infinity with great uniformity from order to higher order and the Gaussian function does not possess genuine zeros in finite regions of the complex z-plane. In Section 3 we established in (3.13) a limiting transition from modified Bessel functions with stretched arguments to a Gauss function and it was interesting to see how the zeros of finite-order Taylor series approximations move from approximation to the next higher approximation in the limiting case to infinity. The pictures for the Gaussian functions are very homogenous in contrast to the considered functions with genuine zeros. This provides an uncommon view onto the absence of zeros of Gaussian functions exp ( a 2 z 2 4 ) despite a representation in the form (2.7).

The zeros in each approximation for the considered functions are either pairs ( z k ) ≡ ( i y k , − i y k ) on the imaginary axis or in majority quadruples ( z k ) ≡ ( x k + i y k , x k − i y k , − x k + i y k , − x k − i y k ) and there was no doubt which are pairs and which are quadruples even in case that their real part x k is small compared with maximal modulus of the zeros and since the whole number of zeros has to be 2M. It is noticeable that new zeros on the imaginary axis when they first appear in the 2M-th approximation may disappear in the next higher ( 2 M + 2 ) -th approximation from the imaginary axis and reappear then in the ( 2 M + 4 ) -th approximation as can be also seen from

To understand the discussed behavior of the zeros from order to next higher order we try to discuss this now in some approximation. We suppose that we have the Taylor series approximation of a Xi function Ξ ( z ) in the order 2M in the form

Ξ 2 M ( z ) = ∑ m = 0 M Ω 2 m z 2 m . (8.1)

We assume that z 0 is an exact zero in 2M-th approximation that means a solution of the equation

∑ m = 0 M Ω 2 m z 0 2 m = 0. (8.2)

Then we try to calculate next higher solutions z = z 0 + Δ z which are near to z 0 from the next higher 2 M + 2 approximation which satisfies the equation

∑ m = 0 M Ω 2 m ( z 0 + Δ z ) 2 m + Ω 2 M + 2 ( z 0 + Δ z ) 2 M + 2 = 0. (8.3)

In full generality this equation would provide the 2 M + 2 solutions of zeros from only one arbitrary solution z 0 in the considered approximation but in such generality we cannot and do not want to solve it. Well soluble is the equation for additions to z 0 leading to a quadratic equation for Δ z . This provides two solutions in a neighborhood for the considered z 0 . Thus we now make an expansion of the left-hand side in (8.3) up to quadratic terms in Δ z

0 = ∑ m = 0 M Ω 2 m z 0 2 m ︸ = 0 + ∑ m = 0 M 2 m Ω 2 m z 0 2 m − 1 Δ z + ∑ m = 0 M m ( 2 m − 1 ) Ω 2 m z 0 2 m − 2 ( Δ z ) 2 + Ω 2 M + 2 ( z 0 2 M + 2 + ( 2 M + 2 ) z 0 2 M + 1 Δ z + ( M + 1 ) 2 M z 0 2 M ( Δ z ) 2 ) . (8.4)

If we neglect from the additional terms of the ( 2 M + 2 ) -th approximation the terms proportional to Δ z and ( Δ z ) 2 as small terms we obtain the following quadratic equation for Δ z z 0

( Δ z z 0 ) 2 + ∑ m = 0 M 2 m Ω 2 m z 0 2 m ∑ m = 0 M m ( 2 m − 1 ) Ω 2 m z 0 2 m Δ z z 0 + Ω 2 M + 2 z 0 2 M + 2 ∑ m = 0 M m ( 2 m − 1 ) Ω 2 m z 0 2 m = 0. (8.5)

The two solutions of this equation are

Δ z z 0 = 1 2 { − ∑ m = 0 M 2 m Ω 2 m z 0 2 m ∑ m = 0 M m ( 2 m − 1 ) Ω 2 m z 0 2 m ± ( ∑ m = 0 M 2 m Ω 2 m z 0 2 m ∑ m = 0 M m ( 2 m − 1 ) Ω 2 m z 0 2 m ) 2 − 4 Ω 2 M + 2 z 0 2 M + 2 ∑ m = 0 M m ( 2 m − 1 ) Ω 2 m z 0 2 m } , (8.6)

or more compactly written

Δ z z 0 = − Ξ 2 M ( 1 ) ( z 0 ) z 0 Ξ 2 M ( 2 ) ( z 0 ) ± ( Ξ 2 M ( 1 ) ( z 0 ) z 0 Ξ 2 M ( 2 ) ( z 0 ) ) 2 − 2 Ω 2 M + 2 z 0 2 M Ξ 2 M ( 2 ) ( z 0 ) . (8.7)

where Ξ 2 M ( n ) ( z 0 ) are the n-th derivatives of Ξ 2 M ( n ) ( z ) taken at z = z 0 . One may express this also by the logarithmic derivative of Ξ 2 M ( 1 ) ( z ) that, however, is inconvenient since it goes into the formula as a denominator.

For the most interesting case of points z 0 = i y 0 on the imaginary axis we find from (8.6)

Δ z y 0 = i 2 { − ∑ m = 0 M ( − 1 ) m 2 m Ω 2 m y 0 2 m ∑ m = 0 M ( − 1 ) m m ( 2 m − 1 ) Ω 2 m y 0 2 m ± ( ∑ m = 0 M ( − 1 ) m 2 m Ω 2 m y 0 2 m ∑ m = 0 M ( − 1 ) m m ( ( 2 m − 1 ) Ω 2 m y 0 2 m ) 2 − 4 ( − 1 ) M + 1 Ω 2 M + 2 y 0 2 M + 2 ∑ m = 0 M ( − 1 ) m m ( 2 m − 1 ) Ω 2 m y 0 2 m } . (8.8)

The sum terms in braces are real ones but likely change their signs at the zeros (we think that it can be proved). The first sum term in braces changes a little the imaginary value of the root but the second term with the root with the two possible signs must be real or imaginary. In case of negative values of the content of the root it gives two imaginary values and together with the whole expression it provides two corrections off the imaginary axis. This can be seen in many of the picture for the roots. However a full discussion of the behavior from approximation to next higher approximation by (8.8) is complicated and has to describe how the roots stabilize on the imaginary axis in dependence on the moments. Such a discussion we cannot give to this time.

We consider here shortly for comparison with the pictures for the up to now discussed functions with a representation of the principal form (2.8) the following unorthodox entire function of an essential other kind

f ( z ) = ∑ n = 0 ∞ z n n ! , (9.1)

with respect to the zeros in its finite Taylor series approximations. This function plays a role for the calculation of the properties of coherent phase states [

The function f ( z ) in (9.1) possesses even and odd powers of variable z and, therefore it is not symmetrical with mirror symmetry to the real and imaginary axis but only symmetrical to the real axis. Taking separately the even and odd powers of z and applying the duplication formula for the factorials one may represent (9.1) in the form

f ( z ) = π 1 4 ∑ m = 0 ∞ z 2 m 2 m m ! ( m ! ( m − 1 2 ) ! + z 2 m ! ( m + 1 2 ) ! ) , f ( 0 ) = 1 , (9.2)

with the possible approximations in the coefficients for m ≫ 1 (see (3.11))

m ! ( m − 1 2 ) ! ≈ ( m + 1 4 ) 1 4 , m ! ( m + 1 2 ) ! ≈ ( m + 3 4 ) − 1 4 . (9.3)

This shows that f ( z ) for real positive z = x > 0 grows moderately faster than the function exp ( x 2 2 ) . Therefore, the function exp ( − z 2 2 ) f ( z ) is moderately increasing for real z = x > 0 and moderately decreasing for z = x < 0 with the Taylor series

exp ( − z 2 2 ) f ( z ) = π 1 4 ∑ l = 0 ∞ z 2 l 2 l l ! ∑ m = 0 l ( − 1 ) l − m l ! m ! ( l − m ) ! ( m ! ( m − 1 2 ) ! + z 2 m ! ( m + 1 2 ) ! ) , (9.4)

and, clearly, possesses the same zeros as the function f ( z ) .

We now give a graphical representation of the zeros of the Taylor series approximations

f N ( z ) = ∑ n = 0 N z n n ! , (9.5)

of the function f ( z ) up to N = 120 in

In

Up to now we calculated only the first 4 pairs of zeros with sufficient accuracy [^{5}.

In this article we illustrated the behavior of the zeros for low-order Taylor series approximations of the Xi function Ξ ( z ) as equivalent to the Riemann zeta function ζ ( z ) (in view of their nontrivial zeros) and the same of the modified Bessel functions ( 2 z ) ν I ν ( z ) . By scaling of the variable z in the modified Bessel functions we derived by limiting transition ν → ∞ the Gaussian function proportional to exp ( z 2 ) (with additional scaling of variable z). The considerations do not pretend to be a proof of the Riemann hypothesis but support the belief to its validity.

There are some bothering problems connected with the zeros of the Riemann zeta function and of functions Ξ ( z ) of the kind (2.8) with symmetrical functions Ω ( u ) which are monotonically decreasing up to their vanishing on the u-axis and which are the differences between strictly monotonically and discontinuously vanishing of these functions. The application of the second mean-value theorem leads to the Bessel function of the kind π 2 z I 1 2 ( z ) = sh ( z ) z with only zeros on the imaginary axis. One may imagine how the mapping onto these function by the second mean-value theorem transforms continuously functions with zeros only on the imaginary axis to the mentioned Bessel function but it is difficult to imagine how this is implemented in this Bessel function for functions of the basic integral representation (2.8) with additional zeros outside the imaginary axis but in Appendix A it is shown how this apparent contradiction can be solved by a limiting transition. It seems that the presence of zeros off the imaginary axis is restricted to step-wise discontinuous monotonically decreasing functions with periodic steps where the Ω ( u ) function to sh ( z ) z itself is step-wise discontinuous and this vanishes in the limiting transition to step lengths zero. Is the second mean-value theore applicable to include such cases? We think that it is applicable.

The time is now mature for a final solution of all aspects of the Riemann hypothesis concerning the nontrivial zeros of the Riemann zeta function.

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

Wünsche, A. (2019) Common Properties of Riemann Zeta Function, Bessel Functions and Gauss Function Concerning Their Zeros. Advances in Pure Mathematics, 9, 281-316. https://doi.org/10.4236/apm.2019.93013

In this Appendix we deal with the general case of monotonically decreasing nonnegative real-valued ladder functions Ω ( u ) with equal step lengths which satisfy the condition

0 ≤ u 1 ≤ u 2 ⇒ Ω ( u 1 ) ≥ Ω ( u 2 ) , Ω ( ∞ ) = 0, (A.1)

with small improvements in comparison to [

Thus we consider now the following stepwise constant functions Ω ( u ) with equal step lengths u 0 and with monotonically decreasing discrete amplitude

platforms Ω ( ( n + 1 2 ) u 0 ) , ( n = 0 , 1 , 2 , ⋯ )

Ω ( u ) = ∑ n = 0 ∞ Ω ( ( n + 1 2 ) u 0 ) { θ ( u − n u 0 ) − θ ( u − ( n + 1 ) u 0 ) } , (A.2)

where θ ( x ) denotes the Heaviside step function. The function Ξ ( z ) possesses then the form

Ξ ( z ) = ∫ 0 + ∞ d u Ω ( u ) ch ( u z ) = 1 z ∑ n = 0 ∞ Ω ( ( n + 1 2 ) u 0 ) { sh ( ( n + 1 ) u 0 z ) − sh ( n u 0 z ) } = 1 z ∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } sh ( ( n + 1 ) u 0 z ) . (A.3)

Due to monotonicity all amplitudes Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) in front

of sh ( ( n + 1 ) u 0 z ) are nonnegative. We calculate now these functions more in detail. Using the relation

sh ( ( n + 1 ) w ) = sh ( w ) U n ( ch ( w ) ) , (A.4)

to the Chebyshev polynomial U n ( w ) we find from (A.3)

Ξ ( z ) = sh ( u 0 z ) z ∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } U n ( ch ( u 0 z ) ) = sh ( u 0 z ) z ∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } ∑ k = 0 [ n 2 ] ( − 1 ) k ( n − k ) ! k ! ( n − 2 k ) ! ( 2 ch ( u 0 z ) ) n − 2 k = sh ( u 0 z ) z ∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } ∑ l = 0 n ch ( ( n − 2 l ) u 0 z ) . (A.5)

The stable zeros on the imaginary axis under variation of u_{0} are determined by u 0 y = m π , ( m = ± 1, ± 2, ⋯ ) but due to symmetry y ↔ − y we discuss sometimes so as if we take into account only positive y. We now bring the factor which determines these zeros in (A.5) to the left-hand side and consider the zeros of the second factor by its separation in real and imaginary part

z Ξ ( z ) sh ( u 0 z ) = ∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } ⋅ ∑ l = 0 n { ch ( ( n − 2 l ) u 0 x ) cos ( ( n − 2 l ) u 0 y ) + i sh ( ( n − 2 l ) u 0 x ) sin ( ( n − 2 l ) u 0 y ) } . (A.6)

If we insert for y any zero y m on the imaginary axis determined by u 0 y m = m π with m ≠ 0 then the imaginary part in (A.6) vanishes and we find

( x + i m π ) Ξ ( x + i m π ) sh ( u 0 ( x + i m π ) ) = ∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } ⋅ ∑ l = 0 n ch ( ( n − 2 l ) u 0 x ) ( − 1 ) ( n − 2 l ) m . (A.7)

In dependence on the chosen m and the amplitudes

Ω ( ( n + 1 2 ) u 0 ) , ( n = 0 , 1 , 2 , ⋯ ) one may solve the equations

∑ n = 0 ∞ { Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) } ∑ l = 0 n ( − 1 ) n m ch ( ( n − 2 l ) u 0 x ) = 0 , (A.8)

and all solutions x = x 0 provide zeros z 0 = x 0 + i m π of Ξ ( z ) for the Omega function in (A.2) off the imaginary axis. If we here formally insert m = 0 on the

left-hand side of (A.8) we see using Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) > 0 that

this left-hand side becomes positive with no solution for x at all. If we make the steps u 0 smaller then the first zero z = y 1 on the imaginary axis grows

according to y 1 = π u 0 . For the limiting transition u 0 → 0 and y 1 → ∞ it is

not important that we choose the amplitude differences in the middle between

two steps as Ω ( ( n + 1 2 ) u 0 ) − Ω ( ( n + 3 2 ) u 0 ) and this could be, e.g. also

Ω ( n u 0 ) − Ω ( ( n + 1 ) u 0 ) .

In

y 1 of the ladder approximation becomes larger according to y m = m π u 0 and

goes for u 0 → 0 to infinity and since any zero z 0 = x 0 + i y m off the imaginary axis must possess an imaginary value which is the same as from zeros of the imaginary axis [

the complex plane. This concerns many functions among them also the Xi function to the Riemann hypothesis.

For the modified Bessel function where this is known from other more direct proofs, e.g. [

Ω ( u ) = { 1 − u , for 0 ≤ u ≤ 1 0, for u ≥ 0 ↔ Ξ ( z ) = ch ( z ) − 1 z 2 = 1 2 ( sh ( z 2 ) z 2 ) 2 . (A.9)

For this function Ω ( u ) it is directly seen that Ξ ( z ) possesses zeros only on the imaginary axis. Only for the monotonically decreasing ladder functions (A.2) as mentioned the described limiting transition u 0 → 0 does not lead to other functions than (A.2) and only these functions may possess zeros outside the imaginary axis which y values in addition have to agree with one of the zeros on the imaginary axis.

We could not find a reason why the Second mean-value theorem should not be applicable to the problem of zeros in the Riemann hypothesis and for other described functions.