The zeta function and classical mechanics

In this paper, Guilherme França and André LeClair show that $$\gamma_{y}\sim 2 \pi \left(y-11/8\right)/W\left((y-11/8)e^{-1}\right)$$ where $W$ is the Lambert W function, and $\gamma_{y}$ is the imaginary part of the the $y$th non-trivial zero along the critical line. This formula can then be used in a similar way to Gram points for finding zeta zeros.

When this function replaces $y$ in $\zeta(1/2+iy)$, and the real and imaginary parts are plotted parametrically, the result shares superficial similarities with a limaçon:

where one complete orbit is made approximately in every integer interval.

Firstly, analysing the limaçon in terms of its orbit seen as a particle travelling at constant velocity around the origin, comparing curvature, distance from origin, and non-directional "speed", where

\begin{align} &x(t)\text{:=}\cos (t)+\cos (2 t)+1&\\ &y(t)\text{:=}\sin (t)+\sin (2 t)&\\ \end{align}

the solutions to

\begin{align} t\in\mathbb{R}:\frac{2 \left(x'(t)^2+y'(t)^2\right)^{3/2}}{\left| x'(t) y''(t)-x''(t) y'(t)\right| }=\sqrt{x'(t)^2+y'(t)^2} \end{align}

are clearly

\begin{align} t=2\pi n-\cos^{-1}(-5/4),n\in\mathbb{Z}\\ t=2\pi n+\cos^{-1}(-5/4),n\in\mathbb{Z}\\ t=\dfrac{2}{3}(3\pi n-\pi),n\in\mathbb{Z}\\ t=\dfrac{2}{3}(3\pi n+\pi),n\in\mathbb{Z}\\ \end{align}

Interestingly, the zeta function has similar properties when analysed in this way. Since

\begin{align} |Z(x)|=\sqrt{\Im\left(\zeta \left(\frac{1}{2}+i y\right)\right)^2+\Re\left(\zeta \left(\frac{1}{2}+i y\right)\right)^2} \end{align}

where $Z$ is the Riemann-Siegel Z function, $\sqrt{\Im(f(y))^2+\Re(f(y))^2}$ is zero at $\gamma=2 \pi \left(y-11/8\right)/W\left((y-11/8)e^{-1}\right),$ for

\begin{align} &f(\text{y})\text{:=}\zeta \left(2 i \pi \left(y-11/8\right)/W\left((y-11/8)e^{-1}\right)+1/2\right)&\\ &v(y)\text{:=}\sqrt{\Im\left(f'(y)\right)^2+\Re\left(f'(y)\right)^2}\\ &\kappa (y)\text{:=}\frac{\left| \Re\left(f'(y)\right) \Im\left(f''(y)\right)-\Im\left(f'(y)\right) \Re\left(f''(y)\right)\right| }{\left(\Im\left(f'(y)\right)^2+\Re\left(f'(y)\right)^2\right)^{3/2}}\\ &d(y)\text{:=}\sqrt{\Im(f(y))^2+\Re(f(y))^2} \end{align}

a plot of $\{v(y),2 \pi /\kappa (y),\pi d(y)\}$ with gridlines at $\gamma=2 \pi \left(y-11/8\right)/W\left((y-11/8)e^{-1}\right)$ looks like this:

where it can be seen that $\{v(y),2 \pi /\kappa (y)\}$ intersect at the zeta zeros.

f[y_] := InputForm[Zeta[Rational[1, 2] + Complex[0, 2] Pi (Rational[-11, 8]
+ y)/ProductLog[E^(-1) (Rational[-11, 8] + y)]]]
v[y_] := InputForm[(Im[Derivative[1][f][y]]^2 + Re[Derivative[1][f][y]]^2)^
Rational[1, 2]]
\[Kappa][y_] := InputForm[Abs[Im[Derivative[2][f][y]] Re[Derivative[1][f][y]] -
Im[Derivative[1][f][y]] Re[Derivative[2][f][y]]] (Im[Derivative[1][f][y]]^2 + Re[
Derivative[1][f][y]]^2)^Rational[-3, 2]]
d[y_] := InputForm[(Im[f[y]]^2 + Re[f[y]]^2)^Rational[1, 2]]
r = 2; rr = 15;
zeros = Rest@   Table[y /.      FindRoot[
2 Pi (-11/8 + y)/LambertW[(-11/8 + y)/E] == Im[ZetaZero[n]], {y,n}], {n, r, rr}];
Plot[{v[y], 2 \[Pi]/\[Kappa][y], \[Pi] d[y]}, {y, r, rr},
PlotRange -> All, PlotTheme -> "Garnet",
Filling -> {1 -> {2}, 2 -> {3}}, Axes -> False, Frame -> True,
ImageSize -> 400, GridLines -> {zeros, {}}, Epilog -> {Red, PointSize[Large],
Point[Transpose@{zeros, Table[v[y], {y, zeros}]}]


It can be seen then that there is clearly a close relationship between the curvature of $f(y)$ (ie $\kappa(y)$), the speed of a particle moving along the path of $f(y)$ (ie $v(y)$), and the distance from origin (ie $d(y)$).

My question then is how would I go about finding the solutions to

\begin{align} A=\left\{y\in\mathbb{R}:v(y)=\dfrac{2\pi}{\kappa(y)}\right\} \end{align}

analytically? And for

\begin{align} B=\left\{y\in\mathbb{R}:\gamma=2 \pi \left(y-11/8\right)/W\left((y-11/8)e^{-1}\right)\right\} \end{align}

is it true that $B\in A$ ?

Update

In response to მამუკა ჯიბლაძე's comment, it certainly seems from a quick numerical check that $v(y),2\pi/\kappa(y)$ don't actually intersect at the zeta zeros, but are related to their approximations:

yy = y /. Table[FindRoot[v[y] == 2 Pi/\[Kappa][y], {y, n}], {n, 1, 20}];
y1 = Table[w[yy[[y]]], {y, 1, Length@yy}] // N;
y2 = Table[Im@ZetaZero[y], {y, 1, Length@yy}] // N;
y3 = Table[w[y], {y, 1, Length@yy}] // N;
ListLinePlot[{Table[(Log[y2[[y]]] - Log[y3[[y]]])/(2 Pi y^(3/2)), {y, 1,
Length@yy}], y1 - y2}]


where for

\begin{align} &w(y):=2 \pi \left(y-11/8\right)/W\left((y-11/8)e^{-1}\right)\\ &r_{y}:=\text{solutions to }v(y)=2\pi/\kappa(y) \end{align}

the above is a plot of $(\log(\gamma_{y})-\log(w(y))/(2\pi y^{3/2})$ against $w(r_{y})-\gamma_{y},$ so it would seem that $B\notin A.$

However, I am still unsure how to show this analytically, and since $(\log(\gamma_{y})-\log(w(y))/(2\pi y^{3/2})$ is just a rough guess, I would be interested to know whether there is a more accurate definition of the relationship between $r_{y}$ and $\gamma_{y}$.

• The first equation would imply that asymptotically 100% of the nontrivial Riemann zeros are on the critical line. This is a famous open problem, so I doubt the authors have a valid unconditional proof. – GH from MO Oct 26 '14 at 10:24
• @martin: Yes, and that is how I understood it. The claimed equation for the zeros on the critical line would imply that only 0% of the nontrivial zeros are off the critical line. We certainly don't have a proof of that. – GH from MO Oct 26 '14 at 10:28
• Yes. If there are more than 0% nontrivial zeros off the critial line, then the equation fails for the zeros on the critical line. The reason is simple: the equation is true for all the zeros in the critical strip, because we know asymptotically the number of such zeros up to a given height. – GH from MO Oct 26 '14 at 10:35
• We have no asymptotic for the zeros on the critical line, because we have no idea how many zeros are on the critical line. We know more than 40%, but the ratio might fluctuate between say 67% and 93% as far as I know. Please delete the remark "as GH from MO pointed out below, this only applies to the zeros on the critical line", because I never said this. – GH from MO Oct 26 '14 at 10:41
• @joro: The alleged formula merely says, assuming standard notation, that the imaginary part of the $n$-th zero on the critical line is asymptotically $2\pi n/\log n$. This is true if and only if asymptotically 100% of all nontrivial zeros are on the critical line (when ordered by imaginary part). In particular, the alleged formula is true under RH. These are all simple facts. – GH from MO Oct 26 '14 at 20:41

protected by Todd Trimble♦Oct 26 '14 at 21:12

Thank you for your interest in this question. Because it has attracted low-quality or spam answers that had to be removed, posting an answer now requires 10 reputation on this site (the association bonus does not count).