8
$\begingroup$

how could evaluate with accuracy the function $ \operatorname{arg} \zeta (1/2+is) $ here $ \zeta (s) $ is the 'Riemann Zeta function' on the critical line

I had thought that I could use the 'Riemann Siegel formula' for

$$ \zeta (1/2+ik)e^{i\theta (k)} \text{ with } \theta (k) = \left(\Im\log\Gamma (1/4+ik/2)\right)-(k/2)\log\pi $$

then I simply should take $\Im\log\zeta(1/2+ik) $ after the calculations of the Riemann Zeta function on the critical line.

  • 0
    Do you mean $\displaystyle\left(\Im \log \left(\Gamma(1/4+ik/2)\right)\right)-k\log(\pi)/2$?2012-03-28
  • 0
    yes , of course :) perhaps i did not give the expression well2012-03-28
  • 0
    There's really no better method than Riemann-Siegel on the critical line, I'm afraid. I presume you've looked at the usual references for implementing Riemann-Siegel numerically?2012-03-28

2 Answers 2

6

(This is a response to a comment of the OP dated july 31 from my other answer and not to the main question)

Perhaps that this picture will make things clearer concerning the argument of $\zeta$ :

arg

It represents $\ \frac 1{\pi}\rm{Arg}(\zeta(1/2+it))\ \ $ ($-\frac 1{\pi}\theta(t)$ would produce the same curve but continuous with the discontinues pieces larger than $14$ shifted down $1, 2, 3,\cdots$ units)

Note that for every zero of $\ \zeta\left(\frac 12+it\right)$ (near $14.135,\ 21.022,\ 25.011,\cdots$) we get a jump of $1$ but this jump is not at a fixed value of the argument $\bmod{2\pi}$. This implies that the $N(t)$ 'continuous zeros counting function' in your link (not exactly the argument of $\zeta$ by the way) won't have its fractional part exactly $0$ when crossing a zeta zero (i.e. $N(t)$ is used for the asymptotic proof using the argument principle so don't take the decimals of this 'number of zeros' too seriously!).

The values of $t$ such that $\ \theta(t)=n\pi\ $ (where $\zeta\left(\frac 12+it\right)$ will be real) are called the Gram points and the 'Gram law' is merely Gram's numerical observation that the zeros of the $Z$-function seem to alternate with the Gram points (on the picture $g_0\approx 17.85,\ g_1\approx 23.17,\ g_2\approx 27.67,\cdots$).
(for more see Edwards' book as well as some exceptions to the 'law')

To make a jump of $1$ exactly when crossing a zero we may use the method proposed in this link : subtract $-\theta(t)$ and $\rm{Arg}\left(\zeta\left(\frac 12+it\right)\right)$. This works because an $\rm{Arg}(z)$ function implementation will usually suppose that $z=r e^{i\phi}$ with $r\ge 0$ and $\phi$ in $(-\pi,\pi]$. So that when $z$ crosses $0$ the new small negative $r$ will be replaced by the positive $r=|Z(t)|$ and $\pi$ will be added/subtracted to the phase to take care of the change of sign. The phase will just change smoothly and furnish the 'crossing angle' at the origin where the fast radius modifications take as back as illustrated.

Orbit corresponding to $\zeta\left(\frac 12+it\right)$ in the complex plane : orbit

I must add that the constraint $\phi \in (-\pi,\pi]$ and the possibility for $\theta(t)$ to be near $k\pi$ will sometimes revert the $\pi$ displacement (near $t=300$ and $1000$ if I remember well) but the error doesn't seem to propagate...

I hope this will clarify more,

8

Let's define two real functions $\theta$ and $\operatorname{Z}$ by : $$\zeta\left(\frac 12+it\right)=e^{-i \theta(t)} \operatorname{Z}(t)$$

The $\theta$ function is named 'Riemann–Siegel theta function' (and shouldn't be confused with the 'Riemann–Siegel formula' used for numerical evaluation of the previous '$\operatorname{Z}$ function').

Then this $\theta$ function is directly minus the argument of $\zeta\left(\frac 12+it\right)$ as you hoped!

In the functional equation at $s=\dfrac 12 + i t\ $ we may replace $\zeta\left(\frac 12\pm it\right)$ by $e^{-i \theta(\pm t)} \operatorname{Z}(\pm t)$ to get (with some work and using the fact that $\operatorname{Z}(t)$ is real so that $\operatorname{Z}(t)$= $\overline{\operatorname{Z}(t)}$) :

$$\theta(t)= \Im\left(\log\left(\Gamma\left(\frac 14+\frac{it}2\right)\right)\right) - \frac {\log(\pi)}2 t$$

The searched argument is simply $-\theta(t)= -\arg\left(\Gamma\left(\frac 14+\frac{it}2\right)\right) + \frac {\log(\pi)}2 t$.

The difficult part is really the $\operatorname{Z}$ function!

  • 0
    so, i can use with no problem the Riemann-Siegel function to evaluate $ Imlog\zeta (1/2+is) $ on the critical line :) can't i ?2012-03-29
  • 0
    @JoseGarcia, looks like. +1 nice answer.2012-03-29
  • 1
    @JoseGarcia: Yes you can! :-) Let's add that it is important to understand that $\operatorname{Z}(t)$ will change sign while 'crossing' the zeros (with no visible effect on the very regular $\theta$ function). In fact we may use this behavior as opposed to having a $|\operatorname{Z}(t)|$ with the phase changed by $\pi$ at the zeros to count the zeros. See [explicit formula for Riemann zeros counting function](http://mathoverflow.net/questions/82635/explicit-formula-for-riemann-zeros-counting-function)2012-03-29
  • 0
    @draks: Thanks for that!2012-03-29
  • 0
    however since $ Z(t) $ is real then the argument of $ \zeta (1/2+it) $ is just the function $ \theta (t) $ but this theta function is just the same as the smooth part of the zeros ¡¡2012-07-31
  • 0
    @JoseGarcia: yes the argument is $-\theta(t)$. I don't know what you call the smooth part (since $0 e^{-i\theta}=0$) but the Arg (arg) at a zero will be the angle on this [picture](http://en.wikipedia.org/wiki/File:Zeta_polar.svg) when the $\zeta(s)$ curve crosses the origin.2012-07-31
  • 0
    but preciselyt this theta functio n $ \theta (t) $ is just the smooth part of the zeros :S as seen on http://en.wikipedia.org/wiki/Riemann_hypothesis#Number_of_zeros2012-07-31
  • 0
    @JoseGarcia: well it allows to count the zeros (as I said and linked to just a little higher) by dividing by $\pi$ or observing the 'jump of Arg'. So what do you want to know?2012-07-31
  • 0
    umm is there or are there graphics for $ arg\zeta (1/2+it) $ or graphical representaiton of the function $ S(T)= \frac{1}{\pi}arg\zeta (1/2+iT) $ for several values of T , can i see or evaluate the eigenvalue staircase for teh Riemann zeros ??2012-07-31
  • 0
    [A derivation of the formula for $\,\theta(t)$](http://math.stackexchange.com/questions/1483780/argument-of-the-riemann-zeta-function-on-res-1/1495533#1495533).2015-11-02