Riesz function

Riesz(x) for x from 0 to 50

In mathematics, the Riesz function is an entire function defined by Marcel Riesz in connection with the Riemann hypothesis, by means of the power series

${\displaystyle {\rm {Riesz}}(x)=-\sum _{k=1}^{\infty }{\frac {(-x)^{k}}{(k-1)!\zeta (2k)}}.}$

If we set ${\displaystyle F(x)={\frac {1}{2}}{\rm {Riesz}}(4\pi ^{2}x)}$ we may define it in terms of the coefficients of the Laurent series development of the hyperbolic (or equivalently, the ordinary) cotangent around zero. If

${\displaystyle {\frac {x}{2}}\coth {\frac {x}{2}}=\sum _{n=0}^{\infty }c_{n}x^{n}=1+{\frac {1}{12}}x^{2}-{\frac {1}{720}}x^{4}+\cdots }$

then F may be defined as

${\displaystyle F(x)=\sum _{k=1}^{\infty }{\frac {x^{k}}{c_{2k}(k-1)!}}=12x-720x^{2}+15120x^{3}-\cdots }$

The values of ζ(2k) approach one for increasing k, and comparing the series for the Riesz function with that for ${\displaystyle x\ \exp(-x)}$ shows that it defines an entire function. Alternatively, F may be defined as

${\displaystyle F(x)=\sum _{k=1}^{\infty }{\frac {k^{\overline {k+1}}x^{k}}{B_{2k}}}.\ }$

$n^{\overline{k}}$ denotes the rising factorial power in the notation of D. E. Knuth and the number B_n are the Bernoulli number. The series is one of alternating terms and the function quickly tends to minus infinity for increasingly negative values of x. Positive values of x are more interesting and delicate.

Riesz criterion

It can be shown that

${\displaystyle \operatorname {Riesz} (x)=O(x^{e})\qquad ({\text{as }}x\to \infty )}$

for any exponent e larger than 1/2, where this is big O notation; taking values both positive and negative. Riesz showed that the Riemann hypothesis is equivalent to the claim that the above is true for any e larger than 1/4.^[1] In the same paper, he added a slightly pessimistic note too: «Je ne sais pas encore decider si cette condition facilitera la vérification de l'hypothèse».

Mellin transform of the Riesz function

The Riesz function is related to the Riemann zeta function via its Mellin transform. If we take

${\displaystyle {\mathbf {M} }({\rm {Riesz}}(z))=\int _{0}^{\infty }{\rm {Riesz}}(z)z^{s}{\frac {dz}{z}}}$

we see that if ${\displaystyle \Re (s)>-1}$ then

${\displaystyle \int _{0}^{1}{\rm {Riesz}}(z)z^{s}{\frac {dz}{z}}}$

converges, whereas from the growth condition we have that if ${\displaystyle \Re (s)<-{\frac {1}{2}}}$ then

${\displaystyle \int _{1}^{\infty }{\rm {Riesz}}(z)z^{s}{\frac {dz}{z}}}$

converges. Putting this together, we see the Mellin transform of the Riesz function is defined on the strip ${\displaystyle -1<\Re (s)<-{\frac {1}{2}}}$. On this strip, we have (cf. Ramanujan's master theorem) ${\displaystyle {\frac {\Gamma (s+1)}{\zeta (-2s)}}={\mathbf {M} }({\rm {Riesz}}(z))}$

From the inverse Mellin transform, we now get an expression for the Riesz function, as

${\displaystyle {\rm {Riesz}}(z)=\int _{c-i\infty }^{c+i\infty }{\frac {\Gamma (s+1)}{\zeta (-2s)}}z^{-s}ds}$

where c is between minus one and minus one-half. If the Riemann hypothesis is true, we can move the line of integration to any value less than minus one-fourth, and hence we get the equivalence between the fourth-root rate of growth for the Riesz function and the Riemann hypothesis.

J. garcia (see references) gave the integral representation of $f(x)$ using Borel resummation as

${\displaystyle \exp(-x)-1=\int _{0}^{\infty }{\frac {f(t)}{t}}\rho ({\sqrt {x/t}})dt}$ and ${\displaystyle \rho (x)=x-\lfloor x\rfloor }$ is the fractional part of 'x'

Calculation of the Riesz function

The Maclaurin series coefficients of F increase in absolute value until they reach their maximum at the 40th term of -1.753×10¹⁷. By the 109th term they have dropped below one in absolute value. Taking the first 1000 terms suffices to give a very accurate value for ${\displaystyle F(z)}$ for ${\displaystyle |z|<9}$. However, this would require evaluating a polynomial of degree 1000 either using rational arithmetic with the coefficients of large numerator or denominator, or using floating point computations of over 100 digits. An alternative is to use the inverse Mellin transform defined above and numerically integrate. Neither approach is computationally easy.

Another approach is to use acceleration of convergence. We have

${\displaystyle {\rm {Riesz}}(x)=\sum _{k=1}^{\infty }{\frac {(-1)^{k+1}x^{k}}{(k-1)!\zeta (2k)}}.}$

Since ζ(2k) approaches one as k grows larger, the terms of this series approach

${\displaystyle \sum _{k=1}^{\infty }{\frac {(-1)^{k+1}x^{k}}{(k-1)!}}=x\exp(-x)}$. Indeed, Riesz noted that: ${\displaystyle \ {\sum _{n=1}^{\infty }{\rm {Riesz}}(x/n^{2})=x\exp(-x)}.}$

Using Kummer's method for accelerating convergence gives

${\displaystyle {\rm {Riesz}}(x)=x\exp(-x)-\sum _{k=1}^{\infty }\left(\zeta (2k)-1\right)\left({\frac {(-1)^{k+1}}{(k-1)!\zeta (2k)}}\right)x^{k}}$

with an improved rate of convergence.

Continuing this process leads to a new series for the Riesz function with much better convergence properties:

${\displaystyle {\rm {Riesz}}(x)=\sum _{k=1}^{\infty }{\frac {(-1)^{k+1}x^{k}}{(k-1)!\zeta (2k)}}=\sum _{k=1}^{\infty }{\frac {(-1)^{k+1}x^{k}}{(k-1)!}}\left(\sum _{n=1}^{\infty }\mu (n)n^{-2k}\right)}$

${\displaystyle \sum _{k=1}^{\infty }\sum _{n=1}^{\infty }{\frac {(-1)^{k+1}\left(x/n^{2}\right)^{k}}{(k-1)!}}=x\sum _{n=1}^{\infty }{\frac {\mu (n)}{n^{2}}}\exp \left(-{\frac {x}{n^{2}}}\right).}$

Here μ is the Möbius mu function, and the rearrangement of terms is justified by absolute convergence. We may now apply Kummer's method again, and write

${\displaystyle {\rm {Riesz}}(x)=x\left({\frac {6}{\pi ^{2}}}+\sum _{n=1}^{\infty }{\frac {\mu (n)}{n^{2}}}\left(\exp \left(-{\frac {x}{n^{2}}}\right)-1\right)\right)}$

the terms of which eventually decrease as the inverse fourth power of n.

The above series are absolutely convergent everywhere, and hence may be differentiated term by term, leading to the following expression for the derivative of the Riesz function:

${\displaystyle {\rm {Riesz}}'(x)={\frac {\rm {Riesz(x)}}{x}}-x\left(\sum _{n=1}^{\infty }{\frac {\mu (n)}{n^{4}}}\exp \left(-{\frac {x}{n^{2}}}\right)\right)}$

which may be rearranged as

${\displaystyle {\rm {Riesz}}'(x)={\frac {\rm {Riesz(x)}}{x}}+x\left(-{\frac {90}{\pi ^{4}}}+\sum _{n=1}^{\infty }{\frac {\mu (n)}{n^{4}}}\left(1-\exp \left(-{\frac {x}{n^{2}}}\right)\right)\right).}$

Marek Wolf in^[2] assuming the Riemann Hypothesis has shown that for large x:

${\displaystyle {\rm {Riesz}}(x)=\mathrm {const} \times x^{1/4}\sin \left(\phi -{\frac {1}{2}}\gamma _{1}\log(x)\right)}$

where ${\displaystyle \gamma _{1}=14.13472514...}$ is the imaginary part of the first nontrivial zero of the zeta function, ${\displaystyle \mathrm {const} =7.7750627...\times 10^{-5}}$ and ${\displaystyle \phi =-0.54916...=(-31,46447^{\circ })}$. It agrees with the general theorems about zeros of the Riesz function proved in 1964 by Herbert Wilf.^[3]

Appearance of the Riesz function

A plot for the range 0 to 50 is given above. So far as it goes, it does not indicate very rapid growth and perhaps bodes well for the truth of the Riemann hypothesis.

Notes

1. M. Riesz, «Sur l'hypothèse de Riemann», Acta Mathematica, 40 (1916), pp.185-90.». For English translation look here
2. M. Wolf, "Evidence in favor of the Baez-Duarte criterion for the Riemann Hypothesis", Computational Methods in Science and Technology, v.14 (2008) pp.47-54
3. H.Wilf, " On the zeros of Riesz' function in the analytic theory of numbers", Illinois J. Math., 8 (1964), pp. 639-641

References

• Titchmarsh, E. C., The Theory of the Riemann Zeta Function, second revised (Heath-Brown) edition, Oxford University Press, 1986, [Section 14.32]