Fonction thêta de Riemann-Siegel

En mathématiques, la fonction thêta de Riemann – Siegel est définie en termes de la fonction gamma :

$\theta (t)=\arg \left(\mathrm{\Gamma }(\frac{1}{4}+\frac{\mathrm{i}t}{2})\right)-\frac{\log \pi }{2}t$

pour t réel. Ici, l'argument est choisi de manière à obtenir une fonction continue et $\theta (0)=0$, c'est-à-dire de la même manière que la branche principale de la fonction log-gamma.

Son développement asymptotique est

$\theta (t)\sim \frac{t}{2}\log \frac{t}{2\pi }-\frac{t}{2}-\frac{\pi }{8}+\frac{1}{48t}+\frac{7}{5760t^3}+\cdots$

qui n'est pas convergent, mais dont les premiers termes donnent une bonne approximation pour $t\gg 1$. Sa série de Taylor en 0 qui converge pour $|t|<1/2$ est

$\theta (t)=-\frac{t}{2}\log \pi +\sum _{k=0}^{\mathrm{\infty }}\frac{(-1{)}^k{\psi }^{(2k)}\left(\frac{1}{4}\right)}{(2k+1)!}{\left(\frac{t}{2}\right)}^{2k+1}$

où ${\psi }^{(2k)}$ désigne la fonction polygamma d'ordre $2k$. La fonction thêta de Riemann–Siegel est intéressante pour étudier la fonction zêta de Riemann, car elle peut faire pivoter la fonction zêta de Riemann de sorte qu'elle devienne la fonction Z à valeur totalement réelle sur la droite critique $s=1/2+\mathrm{i}t$.

Analyse graphique

La fonction thêta de Riemann-Siegel est une fonction analytique réelle impaire sur les $t$ réels avec trois racines $0$ et $\pm 17,8455995405\dots$. C'est une fonction croissante pour $|t|>6.29$, et a des extrema locaux à $\pm 6,289835988\dots$, de valeur $\mp 3,530972829\dots$. Son seul point d'inflexion est en $t=0$ avec ${\theta }^{\prime }(0)=-\frac{\ln \pi +\gamma +\pi /2+3\ln 2}{2}=-2,6860917\dots$, qui est le minimum de sa dérivée.

Thêta comme fonction de la variable complexe

Nous avons une expression de série infinie pour la fonction log-gamma

$\log \mathrm{\Gamma }(z)=-\gamma z-\log z+\sum _{n=1}^{\mathrm{\infty }}(\frac{z}{n}-\log (1+\frac{z}{n})),$

où γ est la constante d'Euler. En remplaçant z par $(2\mathrm{i}t+1)/4$ et en prenant la partie imaginaire terme à terme donne la série suivante pour θ(t)

$\theta (t)=-\frac{\gamma +\log \pi }{2}t-\arctan 2t+\sum _{n=1}^{\mathrm{\infty }}(\frac{t}{2n}-\arctan \left(\frac{2t}{4n+1}\right)).$

Pour les valeurs avec une partie imaginaire entre −1 et 1, la fonction arctangente est holomorphe, et on voit facilement que la série converge uniformément sur les compacts dans la région de partie imaginaire comprise entre −1/2 et 1/2, donnant à une fonction holomorphe sur ce domaine. Il s'ensuit que la Modèle:Lien est également holomorphe dans cette région, qui est la bande critique.

Nous pouvons utiliser les identités

$\arg z=\frac{\log z-\log \overline{z}}{2\mathrm{i}}\text{et}\overline{\mathrm{\Gamma }(z)}=\mathrm{\Gamma }(\overline{z})$

pour obtenir l'expression de forme fermée

$\theta (t)=\frac{\log \mathrm{\Gamma }\left(\frac{2\mathrm{i}t+1}{4}\right)-\log \mathrm{\Gamma }\left(\frac{-2\mathrm{i}t+1}{4}\right)}{2i}-\frac{\log \pi }{2}t=-\frac{i}{2}(\ln \mathrm{\Gamma }(\frac{1}{4}+\frac{\mathrm{i}t}{2})-\ln \mathrm{\Gamma }(\frac{1}{4}-\frac{it}{2}))-\frac{\ln (\pi )t}{2}$

qui étend notre définition originale à une fonction holomorphe en t.

Fonction thêta de Riemann – Siegel dans le plan complexe

$-1<\mathrm{\Re }(t)<1$	$-5<\mathrm{\Re }(t)<5$	$-40<\mathrm{\Re }(t)<40$

Points de Gram

La fonction zêta de Riemann sur la ligne critique peut s'écrire $\zeta (\frac{1}{2}+\mathrm{i}t)={\mathrm{e}}^{-\mathrm{i}\theta (t)}Z(t),$ où $Z(t)={\mathrm{e}}^{\mathrm{i}\theta (t)}\zeta (\frac{1}{2}+\mathrm{i}t).$

Si $t$ est un nombre réel, alors la fonction $Z(t)$ renvoie des valeurs réelles. Par conséquent, la fonction zêta sur la ligne critique sera réelle lorsque $\sin (\theta (t))=0$. Les valeurs réelles positives de $t$ où cela se produit sont appelés points de Gram, d' après J. P. Gram, et peuvent aussi être décrits comme les points où $\frac{\theta (t)}{\pi }$ est un entier. Un point de Gram est une solution $g_n$ de $\theta (g_n)=n\pi .$ Ces solutions sont approchées par la suite :

$g{\text{'}}_n=\frac{2\pi (n+1-\frac{7}{8})}{W\left(\frac{1}{\mathrm{e}}(n+1-\frac{7}{8})\right)},$

où $W$ est la fonction W de Lambert. Voici les plus petits points de Gram positifs

$n$	$g_n$	$\theta (g_n)$
−3	0	0
−2	3,4362182261. . .	Modèle:MathPi
−1	9,6669080561. . .	Modèle:MathPi
0	17.8455995405. . .	0
1	23,1702827012. . .	Modèle:MathPi
2	27,6701822178. . .	2Modèle:MathPi
3	31,7179799547. . .	3Modèle:MathPi
4	35,4671842971. . .	4π
5	38,9992099640. . .	5Modèle:MathPi
6	42,3635503920. . .	6Modèle:MathPi
7	45,5930289815. . .	7Modèle:MathPi
8	48,7107766217. . .	8Modèle:MathPi
9	51,7338428133. . .	9Modèle:MathPi
10	54,6752374468. . .	10Modèle:MathPi
11	57,5451651795. . .	11Modèle:MathPi
12	60,3518119691. . .	12Modèle:MathPi
13	63,1018679824. . .	13Modèle:MathPi
14	65,8008876380. . .	14Modèle:MathPi
15	68,4535449175. . .	15Modèle:MathPi

Les points Gram sont utiles pour calculer les zéros de $Z(t)$. À un point de Gram $g_n,$

$\zeta (\frac{1}{2}+i{g}_n)=\cos (\theta (g_n))Z(g_n)=(-1{)}^{n}Z(g_n),$

et si celle-ci est positive en deux points de Gram successifs, $Z(t)$ doit avoir s'annuler dans l'intervalle.

Selon la loi de Gram, la partie réelle est généralement positive tandis que la partie imaginaire alterne avec les points de Gram, entre des valeurs positives et négatives à intervalles assez réguliers.

$(-1{)}^{n}Z(g_n)>0$

Le nombre de racines, $N(T)$, dans la bande 0 à T, est donné par

$N(T)=\frac{\theta (T)}{\pi }+1+S(T),$

où $S(T)$ est un terme d'erreur qui croît asymptotiquement comme $\log T$.

Si $g_n$ obéissait à la loi de Gram, alors le nombre de racines dans la bande deviendrait simplement

$N(g_n)=n+1.$

Or lorsque T croit, on sait que la loi de Gram échoue pour environ 1/4 de tous les intervalles de Gram, le premier échec est à l'indice 126, avant le 127Modèle:E zéro.

Références

• Modèle:Ouvrage
• Gabcke, W. (1979), Neue Herleitung und explizierte Restabschätzung der Riemann-Siegel-Formel . Thèse, Université de Göttingen. Version révisée (eDiss Göttingen 2015)
• Modèle:Ouvrage

Liens externes

• Modèle:MathWorld
• Wolfram Research – Riemann-Siegel Theta function