Technische Samenvatting: Eenvoudige en Nauwkeurige Benaderingen van de Riemann-zètafunctie

Probleemstelling

De berekening van de Riemann-zètafunctie, $\zeta(s)$, voor grote waarden van het imaginaire deel, $t = \text{Im}(s)$, vormt aanzienlijke uitdagingen. Hoewel de Euler-Maclaurin-sommatiemethode strikte foutgrenzen biedt, schaalt de rekencomplexiteit als $O(|s|)$, waardoor deze onpraktisch is voor grote $t$. De standaardbenadering voor grote $t$ is de Riemann-Siegel-formule, die de complexiteit reduceert tot $O(t^{1/2})$ door gebruik te maken van een hoofdsom en een asymptotische expansie voor de restterm. Echter, hogere-orde correcties in de Riemann-Siegel-formule behelzen complexe afgeleiden van speciale functies (specifiek de functie $F(z)$), die moeilijk numeriek te evalueren zijn vanwege singulariteiten en ingewikkelde formules. Bovendien vereisen bestaande algoritmes met hoge precisie, zoals die van Galway, vaak het afstemmen van meerdere parameters om optimale nauwkeurigheid te bereiken.

Het artikel adresseert de behoefte aan benaderingen die zowel rekenkundig efficiënt zijn (met weinig termen) als uiterst nauwkeurig over willekeurige verticale stroken, inclusief de kritieke strook, zonder de implementatieproblemen die gepaard gaan met hogere-orde Riemann-Siegel-correcties.

Methodologie

De auteur ontwikkelt een familie van benaderingen, aangeduid als $\zeta_p(s)$, die de hoofdsom van de Riemann-Siegel-formule combineren met een nieuwe benadering van de restterm. De kernmethodologie rust op de volgende stappen:

1. Integraalvoorstelling: De restterm van de exacte Riemann-Siegel-formule wordt uitgedrukt als een integraal $I_M(s)$ die een complexwaardige maat $\eta(dx) = \frac{e^{-\pi x^2}}{\cosh(\pi \theta x)} dx$ bevat, waarbij $\theta = e^{-i\pi/4}$.
2. Gaussische Kwanturatie: In plaats van asymptotische reeksexpansies te gebruiken, benadert de auteur deze integraal met een Gaussische kwadratuurregel. Het doel is een set knopen $\{x_{p,j}\}$ en gewichten $\{\omega_{p,j}\}$ te vinden zodat de kwadratuur exact is voor een specifieke set van $4p+2$ testfuncties.
3. Afleiding van Coëfficiënten: De coëfficiënten $\{\omega_{p,j}\}$ en $\{\lambda_{p,j}\}$ (waarbij $\lambda_{p,j} = \theta x_{p,j}$) worden afgeleid door een stelsel lineaire vergelijkingen op te lossen. Dit stelsel wordt geconstrueerd door te eisen dat de kwadratuur exact is voor functies van de vorm $f_k(x) = e^{-\pi x^2 + 2\pi \theta y_{p,k} x}$, waarbij $y_{p,k}$ gelijkmatig verdeelde punten zijn in $[-1, 1]$.
   • Het probleem wordt getransformeerd tot het vinden van de knopen en gewichten van een Gaussische kwadratuur voor een lineaire functionaal $L$ gedefinieerd door de momenten van de Mordell-integraal $H(y)$.
   • Deze knopen en gewichten worden berekend met behulp van orthogonale polynomen gegenereerd via een drietermige recurrenterelatie.
4. Benaderingsformule: De uiteindelijke benadering wordt gedefinieerd als:
   $$\zeta_p(s) = \sum_{n=1}^{N_t} n^{-s} + \chi(s) \sum_{n=1}^{N_t} n^{s-1} - \frac{(-1)^{N_t}}{2} \left[ I_{N_t+1/2, p}(s) + \chi(s) I_{N_t+1/2, p}(1-s) \right]$$
   waarbij $I_{M,p}(s)$ de kwadratuurbenadering is van de integraalterm, die uitsluitend elementaire functies (exponentiële functies, logaritmen) en vooraf berekende complexe coëfficiënten bevat.

Belangrijkste Bijdragen

• Vereenvoudigde Implementatie: De voorgestelde benaderingen hangen af van één geheel getalparameter $p$ en een set vooraf berekende coëfficiënten $\{\omega_{p,j}\}$ en $\{\lambda_{p,j}\}$. In tegenstelling tot de Riemann-Siegel-formule, die voor hogere-orde termen de evaluatie vereist van complexe afgeleiden van speciale functies, vereist $\zeta_p(s)$ uitsluitend elementaire functies en de Gamma-functie.
• Vooraf Berekende Coëfficiënten: De auteur levert vooraf berekende coëfficiënten voor $1 \le p \le 30$ en geselecteerde waarden tot $p=150$. Deze coëfficiënten blijken een matige grootte te hebben en vertonen een snelle afname, wat een stabiele numerieke evaluatie faciliteert.
• Benadering van Afgeleiden: Het artikel breidt de methodologie uit tot het benaderen van de afgeleide $\zeta'(s)$ door de benaderingsformule te differentiëren naar $s$ (geëvalueerd bij vaste $N_t$), waardoor de discontinuïteiten worden vermeden die gepaard gaan met het differentiëren van de vloerfunctie in $N_t$.
• Hybride Algoritme: Er wordt een strategie voorgesteld voor het berekenen van $\zeta(s)$ over het gehele complexe vlak door te schakelen tussen directe sommatie, Euler-Maclaurin-sommatie en de $\zeta_p(s)$-benadering, afhankelijk van het gebied van $s$.

Resultaten

Het artikel presenteert uitgebreid numeriek bewijsmateriaal betreffende de nauwkeurigheid van $\zeta_p(s)$:

• Nauwkeurigheid in de Kritieke Strook: Numerieke experimenten in de kritieke strook ($0 \le \sigma \le 1$) tonen aan dat de fout $|\zeta_p(s) - \zeta(s)|$ aanzienlijk afneemt naarmate $p$ toeneemt.
  • Voor $p=10$ zijn de fouten $< 10^{-15}$ voor $t > 250$.
  • Voor $p=50$ zijn de fouten $< 10^{-100}$ voor $t > 4000$.
  • Voor $p=150$ zijn de fouten $< 10^{-300}$ voor $t > 6900$.
• Foutpatronen: De fout vertoont een duidelijk patroon waarbij deze verdwijnt aan de eindpunten van intervallen $[t_n, t_{n+1}]$ (waarbij $t_n = 2\pi n^2$) en op $4p$ gelijkmatig verdeelde punten binnen deze intervallen. Dit gedrag wordt analytisch verklaard door het feit dat de dominante foutterm evenredig is met $H_p(B(t)) - H(B(t))$, wat nul is op de kwadratuurknopen.
• Vergelijking met Riemann-Siegel: Voor matige $t$ (bijvoorbeeld $10^2 \le t \le 10^4$) blijkt de benadering $\zeta_7(s)$ nauwkeuriger te zijn dan de Riemann-Siegel-benadering met $K=5$ ($\zeta^{RS}_5$). Hoewel $\zeta^{RS}_5$ een snellere asymptotische afname heeft ($O(t^{-13/4})$) vergeleken met $\zeta_7$ ($O(t^{-1/4})$), vereist $\zeta_7$ aanzienlijk minder vooraf berekende coëfficiënten (15 versus ~100–170 voor hoge precisie) en is eenvoudiger te implementeren.
• Precisiegrenzen: In implementaties met dubbele precisie overheersen afrondingsfouten de benaderingsfout voor $t > 200$. Echter, in kwadruple precisie behouden de benaderingen hoge nauwkeurigheid tot veel grotere $t$ (bijvoorbeeld $t > 10^4$ met fouten $< 10^{-29}$).
• Prestaties: De MATLAB-implementatie van het hybride algoritme wordt gerapporteerd als honderden tot duizenden keren sneller dan de ingebouwde MATLAB-functie zeta.

Betekenis en Beweringen

Het artikel beweert dat deze benaderingen hoge precisie berekening mogelijk maken van $\zeta(s)$ en $\zeta'(s)$ in de kritieke strook en andere verticale stroken voor grote $t$. De betekenis ligt in de afweging tussen implementatiecomplexiteit en nauwkeurigheid: de methode biedt een "eenvoudig" alternatief voor de Riemann-Siegel-formule dat de rekenkundige last van hogere-orde correctietermen vermijdt terwijl het concurrerende of superieure nauwkeurigheid behoudt voor een breed scala aan $t$.

De auteur merkt op dat hoewel er in dit werk geen strikte foutgrenzen voor $\zeta_p(s)$ worden geboden, het numerieke bewijsmateriaal uitgebreid is. De methode wordt gepresenteerd als een praktisch instrument voor berekeningen met hoge precisie, die cruciaal zijn voor wiskundige ontdekkingen zoals het weerleggen van de Mertens-conjectuur. Het artikel sluit af met het beschikbaar stellen van open-source implementaties in MATLAB, Python en Fortran, wat direct gebruik door de onderzoeksgemeenschap faciliteert.