Riemannsche Zeta-Funktion

Die Riemannsche Zeta-Funktion, auch Riemannsche ζ-Funktion oder Riemannsche Zetafunktion (nach Bernhard Riemann), ist eine komplexwertige, spezielle mathematische Funktion, die in der analytischen Zahlentheorie, einem Teilgebiet der Mathematik, eine wichtige Rolle spielt. Erstmals betrachtet wurde sie im 18. Jahrhundert von Leonhard Euler, der sie im Rahmen des Basler Problems untersuchte. Bezeichnet wird sie üblicherweise mit dem griechischen Buchstaben ${\displaystyle \zeta }$ (Zeta).

Ihr Definitionsbereich umfasst alle komplexen Zahlen außer der Zahl ${\displaystyle 1}$. Für Werte ${\displaystyle s\in \mathbb {C} }$ mit Realteil größer als ${\displaystyle 1}$ wird die Riemannsche Zeta-Funktion über eine Dirichlet-Reihe definiert. Sie lautet:

$${\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}={\frac {1}{1^{s}}}+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+\cdots }$$ für ${\displaystyle \operatorname {Re} (s)>1}$.

Mittels analytischer Fortsetzung kann sie zu einer auf ${\displaystyle \mathbb {C} \setminus \{1\}}$ holomorphen Funktion ausgeweitet werden. Sie erfüllt eine wichtige Funktionalgleichung, mit deren Hilfe sie sogar charakterisiert werden kann.

Von großer Bedeutung für die Zahlentheorie ist der Zusammenhang der Zeta-Funktion mit der Primfaktorzerlegung natürlicher Zahlen. Auf dieser Basis konnte Riemann im Jahr 1859 die sehr enge und nicht offensichtliche Beziehung zwischen den Primzahlen und der Lage der Nullstellen der Zeta-Funktion nachweisen. So folgt aus der Tatsache ${\displaystyle \zeta (s)\neq 0}$ für alle komplexen Zahlen ${\displaystyle s}$ mit ${\displaystyle \operatorname {Re} (s)\geq 1}$ bereits, dass die ${\displaystyle n}$-te Primzahl ${\displaystyle p_{n}}$ „recht genau“ den Wert ${\displaystyle n\log(n)}$ hat – genauer gesagt folgt^[1][2]

${\displaystyle \lim _{n\to \infty }{\frac {p_{n}}{n\log(n)}}=1.}$

Hier bezeichnet ${\displaystyle \log(n)}$ den natürlichen Logarithmus von ${\displaystyle n}$. Genauere Informationen über nullstellenfreie Bereiche macht das Bild um die Primzahlverteilung deutlicher. Die bisher unbewiesene Riemannsche Vermutung sagt aus, dass alle nicht-trivialen Nullstellen der Riemannschen Zeta-Funktion den Realteil ${\displaystyle {\tfrac {1}{2}}}$ haben, also auf einer gemeinsamen Geraden liegen. Ob diese Vermutung zutrifft, ist eines der wichtigsten ungelösten Probleme der Mathematik. Aufgrund der Bedeutung der Primzahlen für moderne Kryptosysteme (wie in etwa der RSA-Verschlüsselung) genießt die Riemannsche Vermutung auch außerhalb der reinen Zahlentheorie Aufmerksamkeit.

Das Verhalten der Riemannschen Zeta-Funktion in den Bereichen ${\displaystyle \operatorname {Re} (s)\geq 1}$ und ${\displaystyle \operatorname {Re} (s)\leq 0}$ gilt als gut verstanden. Jedoch sind ihre Eigenschaften innerhalb des kritischen Streifens ${\displaystyle 0<\operatorname {Re} (s)<1}$ weitestgehend unbekannt und Gegenstand bedeutender Vermutungen. Dies betrifft unter anderem die Fragen nach asymptotischem Wachstum in imaginärer Richtung und der für die Zahlentheorie so wichtigen Nullstellenverteilung. Nach heutigem Wissensstand beschreibt die Zeta-Funktion im Streifen ${\displaystyle {\tfrac {1}{2}}<\operatorname {Re} (s)<1}$ im Wesentlichen Chaos. Die Werte der Nullstellen bauen nicht nur Brücken zur Theorie der Primzahlen, sondern höchstwahrscheinlich auch zur modernen Quantenphysik. Weitere Anwendungsgebiete sind die Wahrscheinlichkeitstheorie und die Theorie der automorphen Formen (insbesondere im Feld des Langlands-Programms).

Aus Sicht der algebraischen Zahlentheorie ist die Riemannsche Zeta-Funktion nur ein Spezialfall einer ganzen Klasse sogenannter L-Funktionen. So entspricht sie der zum Trivialen Charakter modulo 1 gehörigen Dirichletschen L-Funktion und der zum Zahlkörper ${\displaystyle \mathbb {Q} }$ (rationale Zahlen) korrespondierenden Dedekindschen Zeta-Funktion.

Wegen der überragenden Bedeutung der Riemannschen Vermutung für die Zahlentheorie und deren Anwendungen bleibt der Themenkreis der Riemannschen Zeta-Funktion ein Gebiet intensiver mathematischer Forschung. Entscheidende Fortschritte erzielten Mathematiker wie zum Beispiel Lindelöf, Hadamard, de La Vallée Poussin, Hardy, Littlewood, Selberg, Woronin und Conrey.

Im ganzen Artikel bezeichnet ${\displaystyle \mathrm {i} }$ die imaginäre Einheit und ${\displaystyle \mathrm {e} =2{,}71828\dots }$ die Eulersche Zahl. Zudem wird ${\displaystyle s}$ häufig als komplexe Variable verwendet, die oft in ${\displaystyle s=\sigma +\mathrm {i} t}$ zerlegt wird.

Des Weiteren wird oft die O-Notation von Landau für die Angabe von Fehlergrößen verwendet. Verhalten sich zwei (unbeschränkte) Funktionen ${\displaystyle f}$ und ${\displaystyle g}$ für wachsendes Argument gleich, gilt also ${\displaystyle \lim _{x\to \infty }{\tfrac {f(x)}{g(x)}}=1}$, so wird dies mit ${\displaystyle f(x)\sim g(x)}$ notiert.

Im Zentrum der Zahlentheorie, jenes Zweiges der Mathematik, der sich mit den Eigenschaften der natürlichen Zahlen 1, 2, 3, 4 … beschäftigt, stehen die Primzahlen 2, 3, 5, 7, 11 … Diese sind ausgezeichnet durch die Eigenschaft, genau zwei Teiler zu haben, nämlich die 1 und sich selbst. Die 1 ist keine Primzahl. Bereits Euklid konnte zeigen, dass es unendlich viele Primzahlen gibt, weshalb die Liste 2, 3, 5, 7, 11 … niemals enden wird.

Die Primzahlen sind gewissermaßen die Atome der ganzen Zahlen, da sich jede positive ganze Zahl eindeutig multiplikativ in solche zerlegen lässt. Zum Beispiel gilt 21 = 3 · 7 und 110 = 2 · 5 · 11. Trotz dieser elementaren Eigenschaft ist nach mehreren Jahrtausenden Mathematikgeschichte bis heute kein Muster bekannt, dem sich die Primzahlen in ihrer Folge unterwerfen. Ihre Natur ist eine der bedeutendsten offenen Fragen der Mathematik.

Auch wenn das detaillierte Verständnis der Sequenz 2, 3, 5, 7, 11 … unerreichbar fern ist, kann man nach Mustern suchen, wenn man den Blick ausweitet. Dabei hilft zum Beispiel die Vorstellung, dass mit Hilfe statistischer Methoden das Verhalten sehr vieler Menschen (zum Beispiel bezüglich des Konsum- und Wahlverhaltens) oft überraschend präzise beschrieben werden kann, obgleich ein einzelner Mensch äußerst komplex ist. Das hat grob gesagt damit zu tun, dass größer werdende relevante Datenmengen immer zuverlässigere Informationen liefern. Im Falle der Primzahlen führt eine solche Ausweitung unter anderem zu der Frage, wie viele Primzahlen es unterhalb einer fest gewählten Zahl gibt.

Zum Beispiel sind nur 4 Primzahlen, nämlich 2, 3, 5 und 7, kleiner als die Zahl 10. Im Falle von 50 gibt es schon 15 kleinere Primzahlen, nämlich

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47.

Ende des 19. Jahrhunderts konnte als Folgerung des Primzahlsatzes eine verblüffend genaue Abschätzung für die Verteilung der Primzahlen bewiesen werden. Der Primzahlsatz wurde bereits im 18. Jahrhundert vom 15-jährigen Gauß (in den Jahren 1792/93) vermutet.^[3] Die Abschätzung wurde bereits vor einem Beweis des Primzahlsatzes von Riemann gegeben und tritt als eine Formel in Erscheinung, die die schnelle Berechnung eines Vorhersagewertes zulässt. Mit dieser Formel kann zu einer gegebenen Zahl die Anzahl der Primzahlen, die kleiner als diese Zahl sind, in vernünftiger Zeit geschätzt werden. Die Formel zur Vorhersage wird prozentual immer genauer, je größer die Zahl gewählt wird (jedoch mit Schwankungen). Beispielsweise liefert sie für den Wert 50 die Prognose 14,97 (es sind tatsächlich 15 Primzahlen, siehe oben), womit der Fehler bei 0,16 Prozent liegt. Weiter sagt sie rund 78.527 Primzahlen unter der Zahl 1.000.000 voraus – tatsächlich sind es 78.498. Dies entspricht einer Abweichung von 0,037 Prozent.

Ein mögliches Werkzeug zum Beweis dieser Formel ist die Riemannsche Zeta-Funktion. Dabei wird ausgenutzt, dass sie das Gesetz der eindeutigen Primfaktorzerlegung in der Sprache der Analysis ausdrückt. Also werden die Eigenschaften der Primzahlen in dieser Funktion versteckt abgespeichert. Erhöht sich das Wissen um die Zeta-Funktion, so auch das Wissen um die Primzahlen, sogar in detaillierteren Fragestellungen. So können viele Primzahltests, wie der von Miller-Rabin unter Annahme der Riemannschen Vermutung bewiesen bzw. verbessert werden.^[4]

Die Nullstellen der Zeta-Funktion erzeugen einen Korrekturterm obiger Formel, der sie in einen exakten Ausdruck umwandelt. Diese dadurch entstehende exakte Formel kennt also die Verteilung der Primzahlen bis ins letzte Detail. Damit gelten die Fragen um die Primzahlen jedoch nicht als gelöst: der Rechenaufwand nimmt mit steigenden Werten sehr stark zu und somit sind praktische Berechnungen mit dieser Formel nicht effektiv. Für numerische Forschung eignen sich im Gegensatz dazu moderne Primzahltests besser. Die exakte Formel ist jedoch von theoretischem Interesse: sie birgt nämlich den Fehlerabstand zwischen der einfachen Vorhersage und der tatsächlichen Primzahlverteilung. Es wird vermutet, dass dieser Fehler (innerhalb des Spektrums aller Möglichkeiten) kleinstmöglich ist. Dabei wäre eine Entschlüsselung dieses Fehlers nicht so sehr für die Numerik von Relevanz. Vielmehr ist die reine Mathematik bestrebt, den bisher verborgenen Grund zu erfahren, weshalb der Fehler (falls zutreffend) so klein wie möglich ausfällt.

Die Primzahlen sind nicht nur Gegenstand der mathematischen Grundlagenforschung, sondern haben auch praktische Anwendungen. So kommen beispielsweise bei Kryptosystemen wie der RSA-Verschlüsselung sehr große Primzahlen zum Einsatz.

Eine mathematische Funktion ist im Grunde wie eine Rechenmaschine. Man gibt einen Wert in die Funktion ein, und diese liefert dann ein Ergebnis in Abhängigkeit vom Eingabewert, zumindest theoretisch. Damit ist gemeint, dass die Funktion an sich nicht rechnet, sondern meist nur eine Rechenvorschrift formelhaft festhält. Einfaches Beispiel für eine Funktion ist die quadratische Funktion, welche die Eingabe mit sich selbst multipliziert. Formelhaft schreibt man dies als ${\displaystyle f(x)=x^{2}}$. Somit ordnet die quadratische Funktion beispielsweise der Zahl ${\displaystyle -2}$ den Wert ${\displaystyle (-2)^{2}}$ zu. Rechnet man dies aus, ergibt sich ${\displaystyle 4}$, also ${\displaystyle f(-2)=4}$.

Im Prinzip funktioniert die Riemannsche Zeta-Funktion genau wie das obere Beispiel, nur dass die Rechenvorschrift etwas komplizierter ist. Um diese zu verstehen, muss das Konzept der unendlichen Reihen bekannt sein. Eine (konvergente) Reihe ist grob gesagt eine Summe von Zahlen, die niemals endet und die sich einer Zahl immer mehr annähert. Ein elementares nicht-triviales Beispiel einer Reihe fußt auf der Zahl ${\displaystyle {\tfrac {1}{9}}}$, die in Dezimalschreibweise nicht geschlossen, sondern nur durch die unendliche periodische Entwicklung

${\displaystyle {\frac {1}{9}}=0{,}11111111111\dots }$

geschrieben werden kann. Schaut man genauer hin, sieht man, dass dies gerade die Summe aller kehrwertiger 10er-Potenzen ist:

${\displaystyle {\frac {1}{10}}+{\frac {1}{100}}+{\frac {1}{1000}}+\cdots =0{,}1+0{,}01+0{,}001+\cdots =0{,}111111\dots }$

Damit sich die unendlich lange Summe einem Wert annähert, muss gewährleistet sein, dass die Summanden „schnell genug klein“ werden.

Die Riemannsche Zeta-Funktion ist nun vergleichbar mit einer Rechenmaschine, die zu einer gegebenen Zahl ${\displaystyle x}$ die unendliche Summe der Kehrwerte aller natürlichen Potenzen mit diesem Exponenten bildet. Als mathematische Formel lautet diese Vorschrift

${\displaystyle {\begin{aligned}\zeta (x)&={\frac {1}{1^{x}}}+{\frac {1}{2^{x}}}+{\frac {1}{3^{x}}}+{\frac {1}{4^{x}}}+{\frac {1}{5^{x}}}+{\frac {1}{6^{x}}}+{\frac {1}{7^{x}}}+{\frac {1}{8^{x}}}+{\frac {1}{9^{x}}}+\cdots \\&=\sum _{n=1}^{\infty }{\frac {1}{n^{x}}}.\end{aligned}}}$.

Um das besser zu verstehen, betrachtet man das Beispiel der Eingabe ${\displaystyle x=2}$. Die natürlichen Potenzen mit diesem Exponenten sind gerade die Quadratzahlen 1, 4, 9, 16, 25 … . Damit wäre das Ergebnis der Zeta-Rechenmaschine in Abhängigkeit vom Eingabewert 2 die Reihe

${\displaystyle \zeta (2)={\frac {1}{1}}+{\frac {1}{4}}+{\frac {1}{9}}+{\frac {1}{16}}+{\frac {1}{25}}+{\frac {1}{36}}+{\frac {1}{49}}+{\frac {1}{64}}+{\frac {1}{81}}+\cdots .}$

Es stellt sich heraus, dass die Ausdrücke ${\displaystyle {\tfrac {1}{N^{2}}}}$ schnell genug klein werden, so dass diese unendlich lange Summe sich einem gewissen Zahlenwert beliebig annähert, je weiter man summiert. Experimentell kann man festhalten:

${\displaystyle 1=1;\qquad 1+{\frac {1}{4}}=1{,}25;\qquad 1+{\frac {1}{4}}+{\frac {1}{9}}\approx 1{,}361;\qquad 1+{\frac {1}{4}}+{\frac {1}{9}}+{\frac {1}{16}}\approx 1{,}4236}$

und wenn man sogar bis zur millionsten Quadratzahl ${\displaystyle 1.000.000^{2}=1.000.000.000.000}$ geht

${\displaystyle 1+{\frac {1}{4}}+{\frac {1}{9}}+\cdots +{\frac {1}{1.000.000.000.000}}\approx 1{,}644933.}$

Da hier bereits sehr viele Terme addiert wurden, kann man vermuten, dass ${\displaystyle 1{,}644933}$ dem exakten Ergebnis schon recht nahe ist. Der exakte Grenzwert ist, das konnte Leonhard Euler begründen, die Zahl

${\displaystyle {\frac {\pi ^{2}}{6}}=1{,}644934066848226436472415\dotso }$.

Dabei ist ${\displaystyle \pi =3{,}1415926\dots }$ die Kreiszahl. Zwar kannte Euler den Beginn der Dezimalfolge ${\displaystyle 1{,}64493406684822\dotso }$, doch seine Begründung beruhte letztlich auf mathematischen Argumenten und nicht expliziten Rechnungen, da die Summe ja niemals endet. Somit liefert die Zeta-Funktion für die Eingabe 2 das Ergebnis ${\displaystyle {\tfrac {\pi ^{2}}{6}}}$. Analog müssen für die Eingabe 3, 4 … usw. entsprechend die Kehrwerte aller Kubikzahlen, Biquadratzahlen usw. addiert werden und es kommen neue Grenzwerte ${\displaystyle \zeta (3)}$, ${\displaystyle \zeta (4)}$ … usw. zustande.

Die Zeta-Funktion wird in der Literatur oft über ihre Darstellung als Dirichlet-Reihe definiert.

Für komplexe Zahlen ${\displaystyle s}$, deren Realteil größer als 1 ist, ist die Zeta-Funktion definiert durch die Dirichlet-Reihe^[5]

${\displaystyle \zeta (s):=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}=1+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+{\frac {1}{4^{s}}}+{\frac {1}{5^{s}}}+{\frac {1}{6^{s}}}+{\frac {1}{7^{s}}}+\dotsb ,\,\,\,\,\,\,\,\,\,n^{s}:=\exp(s\log(n)).}$

Wie man mittels des Integralkriteriums für unendliche Reihen zeigen kann, ist diese Reihe im angegebenen Bereich absolut konvergent. Zudem ist die Konvergenz auf kompakten Teilmengen gleichmäßig, weshalb nach dem Satz von Weierstraß die dargestellte Funktion holomorph ist. Wegen der Divergenz der harmonischen Reihe ist diese Darstellung für alle komplexen Zahlen mit Realteil kleiner oder gleich 1 jedoch ungültig. In besonderem Maße wird dies für negative Argumente ersichtlich, wenn man zum Beispiel versuchte, die Zeta-Funktion für ${\displaystyle s=-1}$ über die Dirichlet-Reihe auszuwerten. Man hätte dann

${\displaystyle {\begin{aligned}\zeta (-1)&={\frac {1}{1^{-1}}}+{\frac {1}{2^{-1}}}+{\frac {1}{3^{-1}}}+{\frac {1}{4^{-1}}}+{\frac {1}{5^{-1}}}+{\frac {1}{6^{-1}}}+{\frac {1}{7^{-1}}}+\dotsb \\&=1+2+3+4+5+6+7+\dotsb \end{aligned}}}$

und diese Reihe hat offensichtlich keinen endlichen Grenzwert.

Dennoch wird die Dirichlet-Reihe aufgrund ihrer Einfachheit und ihrer zahlentheoretischen Relevanz (siehe Euler-Produkt) als Basisdefinition verwendet. Mittels analytischer Fortsetzung (siehe unten) wird eine sinnvolle Berechnung für alle komplexen Zahlen ${\displaystyle s}$ mit ${\displaystyle s\neq 1}$ möglich. Damit kann schließlich auch Werten wie ${\displaystyle \zeta (-1)}$ ein Sinn gegeben werden, es gilt zum Beispiel ${\displaystyle \zeta (-1)=-{\tfrac {1}{12}}}$.

Eine wesentliche Eigenschaft der Zeta-Funktion ist ihre Verbindung zu den Primzahlen. Euler, der als Erster diesen Zusammenhang entdeckte, betrachtete dafür das später nach ihm benannte Euler-Produkt, das für alle ${\displaystyle s\in \mathbb {C} }$ mit ${\displaystyle \operatorname {Re} (s)>1}$ gültig ist:^[5]

${\displaystyle \zeta (s)={\frac {1}{\left(1-{\frac {1}{2^{s}}}\right)\left(1-{\frac {1}{3^{s}}}\right)\left(1-{\frac {1}{5^{s}}}\right)\cdots }}=\prod _{p\ {\text{Primzahl}}}{\frac {1}{1-{\frac {1}{p^{s}}}}}=\prod _{p\ {\text{Primzahl}}}\left(1+{\frac {1}{p^{s}}}+{\frac {1}{p^{2s}}}+{\frac {1}{p^{3s}}}+\dotsb \right).}$

Es ist genau äquivalent zur Dirichlet-Reihe und wird von manchen Autoren als Definition verwendet.^[6][7] Jeder einzelne Faktor des Produktes stellt eine geometrische Reihe ${\displaystyle \textstyle \sum q^{n}}$ gebildet über den Wert ${\displaystyle q=p^{-s}}$ dar, während sich das ganze Produkt über alle Primzahlen ${\displaystyle p}$ erstreckt. Das Euler-Produkt ist erstaunlich, weil Primzahlen aufgrund ihrer nicht genau vorhersehbaren Verteilung sehr schwer in analytischen Ausdrücken unterzubringen sind. Es stellt sich aber eine überraschend einfache Identität zwischen den „chaotischen“ Primzahlen und einer bekannten Reihe heraus.

Für die detaillierte Herleitung

Für die formale Herleitung des Euler-Produktes werden lediglich die geometrische Reihe, der Satz, dass jede natürliche Zahl ${\displaystyle n}$ genau eine Zerlegung als Produkt von Primzahlen besitzt, sowie Ausmultiplizieren von Klammern benötigt. Zu Beginn bewährt es sich, nur eine endliche Anzahl von Primzahlen im Produkt zu beachten. Entwickelt man jeden Term ${\displaystyle {\tfrac {1}{1-{\frac {1}{p^{s}}}}}}$ als eine geometrische Reihe ${\displaystyle 1+{\tfrac {1}{p^{s}}}+\left({\tfrac {1}{p^{s}}}\right)^{2}+\left({\tfrac {1}{p^{s}}}\right)^{3}+\cdots }$, so ergibt sich im Falle nur einer Primzahl ${\displaystyle {\frac {1}{1-{\frac {1}{2^{s}}}}}=1+{\frac {1}{2^{s}}}+{\frac {1}{2^{2s}}}+{\frac {1}{2^{3s}}}+\cdots =1+{\frac {1}{2^{s}}}+{\frac {1}{4^{s}}}+{\frac {1}{8^{s}}}+\cdots }$ wobei das Potenzgesetz ${\displaystyle \left({\tfrac {1}{p^{s}}}\right)^{n}={\tfrac {1}{p^{ns}}}={\tfrac {1}{(p^{n})^{s}}}}$ zu beachten ist. Zur Rechten stehen genau die Zahlen, die ausschließlich Zweien in ihrer Primfaktorzerlegung haben, also die Zweierpotenzen. Verfährt man weiter mit den ersten zwei Primzahlen, ergibt sich ${\displaystyle {\frac {1}{1-{\frac {1}{2^{s}}}}}\cdot {\frac {1}{1-{\frac {1}{3^{s}}}}}=\left(1+{\frac {1}{2^{s}}}+{\frac {1}{2^{2s}}}+{\color {red}{\frac {1}{2^{3s}}}}+\cdots \right)\cdot \left(1+{\frac {1}{3^{s}}}+{\color {red}{\frac {1}{3^{2s}}}}+{\frac {1}{3^{3s}}}+\cdots \right).}$ Multipliziert man beide Klammen aus, ergeben sich in der Summe alle Kombinationen von Termen der Form ${\displaystyle {\tfrac {1}{2^{ns}3^{ms}}}}$ mit ${\displaystyle m,n\geq 0}$, es gilt also ${\displaystyle {\frac {1}{1-{\frac {1}{2^{s}}}}}\cdot {\frac {1}{1-{\frac {1}{3^{s}}}}}=1+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+{\frac {1}{2^{s}3^{s}}}+{\frac {1}{2^{3s}}}+{\frac {1}{2^{2s}3^{s}}}+\cdots +{\color {red}{\frac {1}{2^{3s}3^{2s}}}}+\cdots =1+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+{\frac {1}{6^{s}}}+{\frac {1}{8^{s}}}+{\frac {1}{12^{s}}}+\cdots +{\color {red}{\frac {1}{72^{s}}}}+\cdots }$ und auf der rechten Seite stehen genau alle solchen Terme ${\displaystyle n^{-s}}$, sodass ${\displaystyle n}$ nur Zweien und Dreien in seiner Primfaktorzerlegung hat. Beim Ausmultiplizieren wird jeder Summand der einen Klammer mit einem Summand der anderen Klammer verrechnet, und das in jeder Kombination, für ${\displaystyle 72=2^{3}\cdot 3^{2}}$ sind die entsprechenden Terme in Rot markiert. Auf ähnliche Weise findet man, dass ${\displaystyle {\tfrac {1}{1-{\frac {1}{2^{s}}}}}\cdot {\tfrac {1}{1-{\frac {1}{3^{s}}}}}\cdot {\tfrac {1}{1-{\frac {1}{5^{s}}}}}}$ zu der entsprechenden Dirichlet-Reihe korrespondiert, in der alle Zahlen mit Primfaktorzerlegung ${\displaystyle 2^{a}3^{b}5^{c}}$ auftauchen, und so weiter. Entsprechend gilt für allgemein die ersten ${\displaystyle n}$ Primzahlen ${\displaystyle {\frac {1}{1-{\frac {1}{2^{s}}}}}\cdot {\frac {1}{1-{\frac {1}{3^{s}}}}}\cdot {\frac {1}{1-{\frac {1}{5^{s}}}}}\cdot \cdots \cdot {\frac {1}{1-{\frac {1}{p_{n}^{s}}}}}=\sum _{m=2^{a_{1}}3^{a_{2}}5^{a_{3}}\cdots p_{n}^{a_{n}} \atop a_{1},a_{2},a_{3},\dots ,a_{n}\geq 0}{\frac {1}{m^{s}}}.}$ Nun kann man in dieser Formel ${\displaystyle n}$ gegen Unendlich laufen lassen, und erhält ${\displaystyle \prod _{p\ {\text{Primzahl}}}{\frac {1}{1-{\frac {1}{p^{s}}}}}=\lim _{n\to \infty }\sum _{m=2^{a_{1}}3^{a_{2}}5^{a_{3}}\cdots p_{n}^{a_{n}}}{\frac {1}{m^{s}}}=\sum _{m=1}^{\infty }{\frac {1}{m^{s}}}=1+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+{\frac {1}{4^{s}}}+\cdots =\zeta (s),}$ da jede Zahl ${\displaystyle m}$ genau eine Zerlegung ${\displaystyle m=2^{a_{1}}3^{a_{2}}5^{a_{3}}\cdots }$ besitzt.

Das Euler-Produkt konvergiert im betrachteten Bereich ${\displaystyle \operatorname {Re} (s)>1}$ unbedingt.^[8] Mittels des Identitätssatzes für Dirichlet-Reihen lässt sich zeigen, dass das Euler-Produkt und der Fundamentalsatz der Arithmetik zueinander äquivalent sind. Daher wird es zuweilen auch als dessen analytische Version bezeichnet.^[9] Eine wichtige Folgerung des Euler-Produktes für die Analysis der Zeta-Funktion ist, dass ${\displaystyle \zeta (s)\not =0}$ für alle ${\displaystyle \operatorname {Re} (s)>1}$ gilt. Dies ist eine Konsequenz einer Erweiterung des Satzes vom Nullprodukt für unendliche Produkte: keiner der Faktoren des Euler-Produktes ist für irgendeinen Eingabewert aus diesem Bereich Null, also wird es auch nicht im Grenzwert Null sein.^[10] Weitaus nichttrivialer ist die Tatsache, dass das Euler-Produkt, im Gegensatz zur Dirichlet-Reihe, auch auf der Geraden ${\displaystyle \operatorname {Re} (s)=1}$, mit Ausnahme von ${\displaystyle s=1}$, an Gültigkeit behält. Es gilt^[11]

${\displaystyle \prod _{p\ {\text{Primzahl}}}{\frac {1}{1-{\frac {1}{p^{1+it}}}}}=\zeta (1+it),\qquad t\not =0,}$

was die Nullstellenfreiheit von ${\displaystyle \zeta (s)}$ im gesamten Bereich ${\displaystyle \operatorname {Re} (s)\geq 1}$ zur Folge hat. Als eine Folgerung der Funktionalgleichung ergibt sich, dass die einzigen Nullstellen von ${\displaystyle \zeta (s)}$ außerhalb des sog. kritischen Streifens ${\displaystyle 0<\operatorname {Re} (s)<1}$ die trivialen Nullstellen

${\displaystyle -2,\,-4,\,-6,\,-8,\,\dots }$

sind. Alle anderen Nullstellen bezeichnet man als nichttrivial, und diese liegen allesamt im kritischen Streifen.^[12]

Mithilfe des Euler-Produkts der Zeta-Funktion kann ein Beweis des Satzes von Euklid mit analytischen Methoden angegeben werden. Der Satz von Euklid besagt, dass es unendlich viele Primzahlen geben muss, und wurde etwa 300 Jahre vor Christus durch Euklid von Alexandria bewiesen. Unter der Annahme, es gäbe nur endlich viele Primzahlen, gilt

${\displaystyle \zeta (1)=\prod _{p\ {\text{Primzahl}}}{\frac {1}{1-{\frac {1}{p}}}}<\infty ,}$

was ein Widerspruch zur Divergenz der harmonischen Reihe ist. Ähnlich bemerkenswert ist die Argumentation über die Formel

${\displaystyle \prod _{p\ {\text{Primzahl}}}{\frac {1}{1-{\frac {1}{p^{2}}}}}={\frac {\pi ^{2}}{6}}.}$

Bei endlich vielen Primzahlen wäre die linke Seite eine rationale Zahl, die rechte Seite ist aber aufgrund der Transzendenz der Kreiszahl irrational.^[13]

Eine weitere direkte Folgerung des Euler-Produktes, durch Logarithmieren und anschließende Verwendung der Taylor-Reihe des Logarithmus, ist die für ${\displaystyle \operatorname {Re} (s)>1}$ gültige Formel

${\displaystyle \log \zeta (s)=\sum _{n=1}^{\infty }{\frac {P(ns)}{n}},}$

wobei ${\displaystyle P(s)}$ mit ${\displaystyle \textstyle P(s):=\sum _{p\ {\text{Primzahl}}}{\frac {1}{p^{s}}}}$ die Primzetafunktion bezeichnet.^[14] Mit Hilfe von Möbius-Inversion lässt sich daraus eine Möglichkeit ableiten, die Primzetafunktion schnell aus einer Reihe über logarithmierte Zeta-Funktionen zu gewinnen:^[15]

${\displaystyle P(s)=\sum _{n=1}^{\infty }{\frac {\mu (n)}{n}}\log \zeta (ns).}$

Unter anderem kann dieser Ausdruck für eine schnelle numerische Berechnung der Primzetafunktion herangezogen werden.^[16] Ferner folgt aus ${\displaystyle P(s)\sim \log \zeta (s)\to \infty }$ für ${\displaystyle s\to 1}$, dass die Reihe ${\displaystyle {\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{5}}+\dotsb }$ der reziproken Primzahlen divergiert.

Die nach der Definition als Dirichlet-Reihe und dem Euler-Produkt wohl elementarste und wichtigste Darstellung der Zeta-Funktion ist die mit Hilfe eines uneigentlichen Integral-Ausdrucks. Auch diese Darstellung geht direkt aus der Dirichlet-Reihe hervor.

Grundlage dieser Darstellung ist die eulersche Integral-Darstellung der Gamma-Funktion

${\displaystyle \Gamma (s)=\int \limits _{0}^{\infty }\mathrm {e} ^{-x}x^{s-1}\mathrm {d} x,}$

aus dem nach der Substitution ${\displaystyle x=tn}$ mit ${\displaystyle n=1,2,3,\dotsc }$ und Division durch ${\displaystyle n^{s}}$ nach beidseitigem Summieren der Ausdruck

${\displaystyle \zeta (s)\Gamma (s)=\sum _{n=1}^{\infty }\int \limits _{0}^{\infty }\mathrm {e} ^{-nt}t^{s-1}\mathrm {d} t=\int \limits _{0}^{\infty }{\frac {t^{s-1}}{\mathrm {e} ^{t}-1}}\mathrm {d} t}$

hervorgeht.^[17] Diese Darstellung von ${\displaystyle \zeta (s)\Gamma (s)}$ gilt naturgemäß nur auf der Halbebene ${\displaystyle \{s\in \mathbb {C} \mid \operatorname {Re} (s)>1\}}$. Die zweite Integraldarstellung von ${\displaystyle \zeta (s)\Gamma (s)}$ bezeichnet man auch als die Mellin-Transformation von ${\displaystyle {\tfrac {1}{\mathrm {e} ^{t}-1}}}$. Das mögliche Vertauschen von Summe und Integral kann mit absoluter Konvergenz und dem Satz von Lebesgue begründet werden. Eine dazu verwandte Form ist

${\displaystyle \pi ^{-{\frac {s}{2}}}\Gamma \left({\frac {s}{2}}\right)\zeta (s)={\frac {1}{2}}\int \limits _{0}^{\infty }(\vartheta (0,\mathrm {i} t)-1)t^{{\frac {s}{2}}-1}\mathrm {d} t}$

mit der Jacobischen Theta-Funktion (eine Modulform halbganzen Gewichts).

Die Darstellung der Zeta-Funktion mit Hilfe der Gamma-Funktion und der Mellin-Transformation von ${\displaystyle {\tfrac {1}{\mathrm {e} ^{t}-1}}}$ ist daher zentral, da sie ein Ausgangspunkt für die analytische Fortsetzung der Zeta-Funktion ist. Außerdem können mit ihr charakteristische Funktionalgleichungen und die Beziehung zur Theorie der Modulformen hergeleitet werden.

Die anfänglich nur für komplexe Zahlen ${\displaystyle \operatorname {Re} (s)>1}$ definierte Zeta-Funktion kann zu einer in ganz ${\displaystyle \mathbb {C} \setminus \{1\}}$ holomorphen Funktion ausgeweitet werden. Diese Tatsache mag zunächst ungewöhnlich wirken, da ihre Dirichlet-Reihe an vielen Stellen nicht mehr konvergiert. Tatsächlich aber steht die Dirichlet-Reihe (wie auch Euler-Produkt und Mellin-Transformation aus Gründen der Äquivalenz) nicht überall für die Definition der Zeta-Funktion zur Verfügung.

An der Stelle ${\displaystyle s=1}$ besitzt die Zeta-Funktion zunächst mit Sicherheit eine Definitionslücke, denn mit der Divergenz der harmonischen Reihe folgt

${\displaystyle \lim _{\sigma \to 1^{+}}\sum _{n=1}^{\infty }{\frac {1}{n^{\sigma }}}=\infty .}$

Also wird sie in jedem Intervall ${\displaystyle (1,1+\varepsilon )}$ beliebig anwachsen. Diese Lücke bildet gleichzeitig eine natürliche Barriere für die Konvergenz der Dirichlet-Reihe, was aus den Regeln für Abszissen von Dirichlet-Reihen folgt, denn die betrachtete Dirichlet-Reihe hat Konvergenzabszisse ${\displaystyle \sigma _{c}=1}$.

Eine analytische Fortsetzung der im Gebiet ${\displaystyle H:=\{s\in \mathbb {C} \mid \operatorname {Re} (s)>1\}}$ durch die Reihe ${\displaystyle \textstyle \sum n^{-s}}$ definierten holomorphen Funktion ist eine auf einem größeren Gebiet ${\displaystyle H\subsetneq D}$ holomorphe Funktion, die auf ganz ${\displaystyle H}$ mit dieser übereinstimmt. Nach dem Identitätssatz für holomorphe Funktionen ist eine solche Fortsetzung stets eindeutig bestimmt. Damit sind alle Werte der Zeta-Funktion im erweiterten Bereich ${\displaystyle D}$ bereits durch die Dirichlet-Reihe festgelegt, obwohl sie hier nicht mehr an allen Stellen konvergiert.

Obwohl es für den ganz allgemeinen Fall kein konstruktives Verfahren gibt, Berechnungsformeln für analytische Fortsetzungen anzugeben, ist es durch die Einfachheit der Dirichlet-Reihe ${\displaystyle \textstyle \sum n^{-s}}$ nicht schwierig, für die Zeta-Funktion eine zu finden. Besonders einfach erweist sich dies für die gelochte Halbebene

${\displaystyle D=\{s\in \mathbb {C} \mid \operatorname {Re} (s)>0,s\neq 1+{\tfrac {2\pi \mathrm {i} m}{\log 2}},m\in \mathbb {Z} \}}$

mittels folgender Beobachtung:^[18]

${\displaystyle (1-2^{1-s})\zeta (s)=\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}}{n^{s}}}=1-{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}-{\frac {1}{4^{s}}}+{\frac {1}{5^{s}}}-{\frac {1}{6^{s}}}+{\frac {1}{7^{s}}}-{\frac {1}{8^{s}}}+{\frac {1}{9^{s}}}-\dotsb .}$

Die Reihe zur Rechten konvergiert nachweislich in der Halbebene ${\displaystyle \{s\in \mathbb {C} \mid \operatorname {Re} (s)>0\}}$ gegen eine holomorphe Funktion und wird in der Literatur auch manchmal als Dirichletsche Etafunktion ${\displaystyle \eta (s)}$ bezeichnet. Damit lässt sich die Zeta-Funktion zu einer in ganz ${\displaystyle D}$ holomorphen Funktion fortsetzen. Die Lücke in ${\displaystyle s=1}$ wird mittels des Faktors ${\displaystyle 1-2^{1-s}}$ gehoben und muss daher ein Pol erster Ordnung sein. Das Residuum der Zeta-Funktion ist dort 1, das heißt, es gilt:^[19]

${\displaystyle \lim _{s\to 1}\,(s-1)\zeta (s)=1.}$

Alle Stellen ${\displaystyle s_{m}=1+{\tfrac {2\pi \mathrm {i} m}{\log 2}}}$ mit ${\displaystyle m\not =0}$ sind hingegen hebbare Singularitäten, denn es gilt dann ${\displaystyle \eta (s_{m})=0.}$ Dies zeigt man am besten mittels partieller Summation: Für alle ${\displaystyle x>3}$ gilt^[20]

${\displaystyle \sum _{n\leq x}{\frac {(-1)^{n}}{n^{s_{m}}}}=2^{1-s_{m}}\sum _{n\leq x/2}{\frac {1}{n^{s_{m}}}}-\sum _{n\leq x}{\frac {1}{n^{s_{m}}}}=-\sum _{x/2<n\leq x}{\frac {1}{n^{s_{m}}}}=O\left({\frac {1}{x}}\right).}$

Für eine weitere holomorphe Ausdehnung des Definitionsbereiches eignen sich nun viele Methoden, die jedoch nach dem Identitätssatz alle dieselbe Funktion darstellen. Eine davon bietet die Anwendung der Eulerschen Reihentransformation auf die obere alternierende Reihe. Man erhält damit eine von Konrad Knopp veröffentlichte und auf ganz ${\displaystyle \mathbb {C} \setminus \{1+{\tfrac {2\pi \mathrm {i} m}{\log 2}}\mid m\in \mathbb {Z} \}}$ definierte Reihenidentität

${\displaystyle \zeta (s)={\frac {1}{1-2^{1-s}}}\sum _{n=0}^{\infty }{\frac {1}{2^{n+1}}}\sum _{k=0}^{n}(-1)^{k}\ {n \choose k}\ {\frac {1}{(k+1)^{s}}}.}$

Diese wurde 1930 von Helmut Hasse bewiesen.^[21] Es treten daher während der weiteren Fortsetzung keine weiteren Lücken bzw. Pole mehr auf. Daraus folgt schließlich Holomorphie in ${\displaystyle \mathbb {C} \setminus \{1\}}$.

Die Idee der Theorie der Limitierungsverfahren ist es, einem divergenten Grenzprozess einen endlichen Wert zuzuordnen, indem man etwa zusätzliche Parameter einführt, die man anschließend gegen den Ursprungsausdruck „limitiert“. Dies geht bereits auf Leonhard Euler zurück, der berühmt für seinen sorglosen Umgang mit divergenten Reihen ist. Er berechnete einige Werte der Zeta-Funktion annähernd auch außerhalb des Bereichs, in welchem die Dirichlet-Reihe konvergiert. Auf diese Weise stieß er auch auf seine Vermutung hinsichtlich der Gültigkeit ihrer Funktionalgleichung, die er allerdings nicht beweisen konnte.

Der Gedanke ist, der für ${\displaystyle \mathrm {Re} (s)\leq 0}$ divergenten Reihe

${\displaystyle (1-2^{1-s})\zeta (s)=\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}}{n^{s}}}}$

für alle Werte ${\displaystyle s}$ einen „Grenzwert“ zuzuordnen. Dies gelingt durch Einführen eines weiteren Parameters ${\displaystyle 0<t\leq 1}$. Es ist die Reihe

${\displaystyle \phi _{s}(t):=\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}}{n^{s}}}t^{n}}$

für jedes ${\displaystyle s\in \mathbb {C} }$ für alle ${\displaystyle 0<t<1}$ konvergent. Auf dieser Basis kann der Konvergenzbegriff abgeschwächt werden: Eine Reihe ${\displaystyle \textstyle \sum _{n=0}^{\infty }a_{n}}$ heißt A-summierbar, falls die zugehörige Potenzreihe ${\displaystyle \textstyle t\mapsto \sum _{n=1}^{\infty }a_{n}t^{n}}$ für alle ${\displaystyle 0<t<1}$ konvergiert und der Grenzwert ${\displaystyle \textstyle \lim _{t\to 1^{-}}\sum _{n=0}^{\infty }a_{n}t^{n}}$ existiert.^[22] Ist eine Reihe bereits im klassischen Sinne konvergent, stimmen die jeweiligen Grenzwerte nach dem Abelschen Grenzwertsatz überein, doch es gibt A-summierbare Reihen, die nicht konvergieren. Es ist also A-Summierbarkeit eine wohldefinierte Verallgemeinerung der klassischen Reihenkonvergenz. Durch Bilden des Grenzwertes ${\displaystyle t\to 1^{-}}$, d. h., ${\displaystyle t}$ nähert sich von links, erhält man

${\displaystyle \lim _{t\to 1^{-}}\phi _{s}(t)=\lim _{t\to 1^{-}}\left(t-2^{-s}t^{2}+3^{-s}t^{3}-4^{-s}t^{4}+5^{-s}t^{5}-6^{-s}t^{6}+\cdots \right)=(1-2^{1-s})\zeta (s)}$

sogar für alle ${\displaystyle s}$, und die rechte Seite stellt eine ganze Funktion dar.^[23]

Durch die Hinzunahme einer zweiten zu limitierenden Variable entsteht zugleich eine Beziehung zum Polylogarithmus. Dieser verallgemeinert unter anderen den natürlichen Logarithmus und ist für ${\displaystyle |z|<1}$ gegeben durch die Potenzreihe

${\displaystyle \operatorname {Li} _{s}(z):=\sum _{n=1}^{\infty }{\frac {z^{n}}{n^{s}}}.}$

Ist zudem ${\displaystyle \operatorname {Re} (s)>0}$, so ist diese Reihe auch an den Randwerten ${\displaystyle |z|=1}$ (außer bei ${\displaystyle z=1}$) konvergent. Generell ist für feste Werte ${\displaystyle s}$ eine analytische Fortsetzung in ${\displaystyle z}$ auf das Gebiet ${\displaystyle \mathbb {C} \setminus \mathbb {R} _{\geq 1}}$ möglich. Es gilt für alle ${\displaystyle \operatorname {Re} (s)>1}$ die Beziehung ${\displaystyle \zeta (s)=\operatorname {Li} _{s}(1)}$, aber auch

${\displaystyle \zeta (s)={\frac {\operatorname {Li} _{s}(-1)}{2^{1-s}-1}}}$

für Werte ${\displaystyle s\in \mathbb {C} \setminus \{1\}}$. Etwa für ${\displaystyle s=-2}$ gilt ${\displaystyle \textstyle \mathrm {Li} _{-2}(z)=\sum _{n=1}^{\infty }n^{2}z^{n}={\tfrac {z^{2}+z}{(1-z)^{3}}}}$ zunächst für ${\displaystyle |z|<1}$ (siehe Bild rechts), also nach Grenzwertbildung

${\displaystyle \zeta (-2)={\frac {1}{2^{3}-1}}\lim _{z\to -1^{+}}{\frac {z^{2}+z}{(1-z)^{3}}}=0.}$

Eine weitere Möglichkeit, eine analytische Fortsetzung anzugeben, bietet die Euler-Maclaurin-Summenformel. Diese drückt diskrete Summen explizit in der Sprache der Integralrechnung aus und ist allgemein gegeben durch

${\displaystyle {\begin{aligned}\sum _{n=N+1}^{M}f(n)&=\int \limits _{N}^{M}f(x)\,\mathrm {d} x+{\frac {f(M)-f(N)}{2}}+\sum _{k=1}^{\lfloor p/2\rfloor }{\frac {B_{2k}}{(2k)!}}(f^{(2k-1)}(M)-f^{(2k-1)}(N))\\&\quad +(-1)^{p+1}\int \limits _{N}^{M}f^{(p)}(x){\frac {\mathrm {B} _{p}(x-\lfloor x\rfloor )}{p!}}\,\mathrm {d} x.\end{aligned}}}$

Hierbei ist ${\displaystyle f}$ eine auf dem Intervall ${\displaystyle [N,M]}$ mindestens ${\displaystyle p}$-mal differenzierbare Funktion und ${\displaystyle p\geq 1}$ eine natürliche Zahl. Es bezeichnen zudem ${\displaystyle \mathrm {B} _{\nu }(x)}$ die Bernoulli-Polynome und ${\displaystyle \lfloor x\rfloor }$ den ganzzahligen Anteil von ${\displaystyle x}$.^[24]

Mit ${\displaystyle f(x)=x^{-s}}$, ${\displaystyle N>0}$ und ${\displaystyle M\to \infty }$ folgt also^[25]

${\displaystyle \zeta (s)=\sum _{n=1}^{N-1}{\frac {1}{n^{s}}}+{\frac {N^{1-s}}{s-1}}+{\frac {1}{2}}N^{-s}+\sum _{r=1}^{\lfloor p/2\rfloor }{\frac {B_{2r}}{(2r)!}}s(s+1)\cdots (s+2r-2)N^{-s-2r+1}+R_{p}(s).}$

Dabei ist das Restglied gegeben durch^[25]

${\displaystyle R_{p}(s)=-{\frac {s(s+1)\cdots (s+p-1)}{p!}}\int \limits _{N}^{\infty }\mathrm {B} _{p}(x-\lfloor x\rfloor )x^{-s-p}\mathrm {d} x}$

und konvergiert in der gesamten Halbebene ${\displaystyle \operatorname {Re} (s)>1-p}$ (gleichmäßig auf kompakten Teilmengen). Daher stellt diese Formel eine holomorphe Fortsetzung der Zeta-Funktion in die Halbebene ${\displaystyle \operatorname {Re} (s)>1-p}$ dar. Lässt man ${\displaystyle p}$ gegen unendlich gehen, ergibt sich damit ein holomorpher Ausdruck für ganz ${\displaystyle \mathbb {C} \setminus \{1\}}$.^[26]

Setzt man zum Beispiel ${\displaystyle N=p=1}$, ergibt sich die in der Literatur häufig zitierte Darstellung

${\displaystyle \zeta (s)={\frac {s}{s-1}}+s\int \limits _{1}^{\infty }{\frac {\lfloor x\rfloor -x}{x^{s+1}}}\,\mathrm {d} x,}$

die für ${\displaystyle \operatorname {Re} (s)>0}$ gültig ist.^[27]

Eng verwandt mit der Darstellung der Zeta-Funktion mittels der Mellin-Transformation ist eine Darstellung der Funktion mittels eines Kurvenintegrals. Diese wurde von Riemann selbst verwendet, um die Zeta-Funktion in die komplexe Ebene fortzusetzen. Die Funktion ${\displaystyle \textstyle f(s;z)=z^{s-1}/(\mathrm {e} ^{z}-1)}$ ist je nach Wahl des Zweiges des Logarithmus in unterschiedlichen Bereichen holomorph. Für die Hankel-Kontur (einen speziellen Integrationsweg) ist es von Vorteil, die Gerade ${\displaystyle [0,\infty [}$ aus dem Gebiet auszuschließen via:

${\displaystyle {\frac {z^{s-1}}{\mathrm {e} ^{z}-1}}={\frac {\mathrm {e} ^{\pi \mathrm {i} (s-1)+\log(-z)(s-1)}}{\mathrm {e} ^{z}-1}}.}$

Nun definiert man für ${\displaystyle \pi >r>0}$ die Funktion ${\displaystyle I(s)}$ als ein Kurvenintegral über ${\displaystyle f(z)}$. Die gewählte Kurve ${\displaystyle C_{r}}$ kommt von ${\displaystyle +\infty }$, verläuft mit Abstand ${\displaystyle r}$ über der reellen Geraden, umläuft den Ursprung in einem Halbkreis und erstreckt sich dann wieder mit Abstand ${\displaystyle r}$ unterhalb der reellen Geraden gegen ${\displaystyle +\infty }$:

${\displaystyle I(s)=\int _{C_{r}}{\frac {z^{s-1}}{\mathrm {e} ^{z}-1}}\mathrm {d} z.}$

Wegen gleichmäßiger Konvergenz auf kompakten Mengen in ${\displaystyle \mathbb {C} }$ ist ${\displaystyle I(s)}$ eine ganze Funktion. Wählt man nun ${\displaystyle \operatorname {Re} (s)>1}$, so kann man wegen ${\displaystyle \textstyle |f(s;z)|\ll r^{\sigma -2}}$ die Schlaufe beliebig zusammenziehen und erhält mit der Mellin-Transformation

${\displaystyle I(s)=\lim _{r\to 0}(\mathrm {e} ^{2\pi \mathrm {i} s}-1)\int \limits _{r}^{\infty }{\frac {t^{s-1}}{\mathrm {e} ^{t}-1}}\mathrm {d} t=(\mathrm {e} ^{2\pi \mathrm {i} s}-1)\Gamma (s)\zeta (s).}$

Daraus ergibt sich mit dem Ergänzungssatz die Formel^[28]

${\displaystyle \zeta (s)={\frac {\mathrm {e} ^{-\pi \mathrm {i} s}}{2\pi \mathrm {i} }}\Gamma (1-s)I(s).}$

Ist ${\displaystyle s\in \mathbb {Z} }$, so ist ${\displaystyle f(s;z)}$ innerhalb des gelochten Streifens ${\displaystyle \{z\in \mathbb {C} \mid -2\pi <\operatorname {Im} (z)<2\pi \}\setminus \{0\}}$ holomorph. Damit lässt sich die Hankel-Kontur zu einer Kreiskurve zusammenziehen, ohne den Wert des Integrals zu verändern. Dies ermöglicht eine schnelle Berechnung der Werte ${\displaystyle I(n)}$ für ganze Zahlen ${\displaystyle n}$ mittels des Residuensatzes. Unter anderem folgt daraus

${\displaystyle I(n)=0}$

für alle ${\displaystyle n\geq 2}$ und die enge Beziehung der Werte der Zeta-Funktion an nicht-positiven ganzen Argumenten und den Bernoulli-Zahlen.

Diese Darstellungsform kann auch zu einer direkten Herleitung der Funktionalgleichung verwendet werden. Dabei wird die Kurve modifiziert und es kommt der Residuensatz zum Einsatz.^[29]

Erwähnenswert ist der Reihenausdruck

${\displaystyle \zeta (s)={\frac {1}{s-1}}+1-\sum \limits _{n=1}^{\infty }(\zeta (s+n)-1){\frac {s(s+1)\cdots (s+n-1)}{(n+1)!}}}$,

der für alle Werte ${\displaystyle s\neq 1,0,-1,\dotsc }$ definiert ist.^[30] Interessant daran ist, dass sich damit die Zeta-Funktion rekursiv auf die ganze Zahlenebene fortsetzen lässt, da für die Berechnung von ${\displaystyle \zeta (s)}$ lediglich die Werte ${\displaystyle \zeta (s+1),\zeta (s+2),\dotsc }$ benötigt werden.

Von Helmut Hasse stammt die global konvergente Reihe^[31]

${\displaystyle \zeta (s)={\frac {1}{s-1}}\sum _{n=0}^{\infty }{\frac {1}{n+1}}\sum _{k=0}^{n}\left({n \atop k}\right){\frac {(-1)^{k}}{(k+1)^{s-1}}}.}$

Ein exotischer und global konvergenter Ausdruck ergibt sich, wenn man direkt die elementare Reihendarstellung der Zeta-Funktion in die Abel-Plana-Summenformel einsetzt:^[32]

${\displaystyle \zeta (s)={\frac {1}{s-1}}+{\frac {1}{2}}+2\int \limits _{0}^{\infty }{\frac {\sin(s\arctan t)}{(1+t^{2})^{\frac {s}{2}}(\mathrm {e} ^{2\pi t}-1)}}\,\mathrm {d} t.}$

Im Gegensatz zu den Primzahlen oder der euklidischen Geometrie ist die mathematische Entdeckungsgeschichte der Riemannschen Zeta-Funktion sehr jung. So sind alle bis heute wesentlichen Entdeckungen zu dieser Funktion in den letzten 250 Jahren gemacht worden. Auf der einen Seite lässt sich die im Verhältnis zur Entstehung einer rigorosen (komplexen) Analysis frühe Entdeckung mit der Einfachheit der Reihe begründen. Auf der anderen Seite lassen sich die späten Resultate mit der Schwierigkeit von deren Eigenschaften erklären.

Einer der ersten Mathematiker, der sich mit einem Vorläufer der wie heute definierten Zeta-Funktion intensiv und ausführlich auseinandersetzte, war Leonhard Euler. Seit Mitte des 17. Jahrhunderts versuchten Mathematiker, den exakten Grenzwert der unendlichen Reihe

${\displaystyle 1+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+{\frac {1}{5^{2}}}+\dotsb }$

zu bestimmen. Persönlichkeiten wie Pietro Mengoli, der das Basler Problem (wie es später bezeichnet wurde) erstmals formulierte, aber auch Jakob I Bernoulli scheiterten mit ihren Lösungsversuchen. Erst um das Jahr 1734 fand Leonhard Euler die Lösung

${\displaystyle 1+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+\dotsb ={\frac {\pi ^{2}}{6}}=1{,}6449340668\dots }$

mit der Kreiszahl ${\displaystyle \pi =3{,}14159265\dots }$, indem er eine neuartige Technik zur Berechnung der Sinusfunktion entwickelte.^[33] Dieser Beweis wurde jedoch nach Veröffentlichung von seinen Zeitgenossen zunächst nicht akzeptiert. Daraufhin konterte er mit der Veröffentlichung eines alternativen Beweises im Jahr 1741.^[34] Natürlicherweise war Euler bald darauf an der Untersuchung von Reihen des Typs

${\displaystyle f(n)=1+{\frac {1}{2^{n}}}+{\frac {1}{3^{n}}}+{\frac {1}{4^{n}}}+{\frac {1}{5^{n}}}+\dotsb }$

interessiert. Er hatte die Hoffnung, weitere und außerdem weit bedeutendere Aussagen treffen zu können. Und tatsächlich sollte es nicht nur bei der Lösung des Basler Problems bleiben. Er fand unter anderem die Formeln

${\displaystyle f(4)=1+{\frac {1}{2^{4}}}+{\frac {1}{3^{4}}}+{\frac {1}{4^{4}}}+{\frac {1}{5^{4}}}+\cdots ={\frac {\pi ^{4}}{90}},}$

${\displaystyle f(6)=1+{\frac {1}{2^{6}}}+{\frac {1}{3^{6}}}+{\frac {1}{4^{6}}}+{\frac {1}{5^{6}}}+\cdots ={\frac {\pi ^{6}}{945}},}$

die 1735 in seiner Arbeit De Summis Serierum Reciprocarum erstmals veröffentlicht wurden. Obwohl mit steigenden Eingabezahlen die Funktionswerte komplizierter werden, berechnete Euler von Hand^[35] den Wert

${\displaystyle f(26)=1+{\frac {1}{2^{26}}}+{\frac {1}{3^{26}}}+{\frac {1}{4^{26}}}+{\frac {1}{5^{26}}}+\dotsb ={\frac {1.315.862}{11.094.481.976.030.578.125}}\pi ^{26}.}$

In seinem Buch Institutiones calculi differentialis cum eius usu in analysi finitorum ac doctrina serierum, das 1755 veröffentlicht wurde, bewies er schließlich eine allgemeine Formel für ${\displaystyle f(2n)}$.^[36] Diese zeigte auf, dass sich tatsächlich jedes ${\displaystyle f(2n)}$ stets als ein rationales Vielfaches der Potenz ${\displaystyle \pi ^{2n}}$ schreiben lässt. Nicht erfolgreich war er hingegen bei ungeraden Argumenten, also zum Beispiel bei der Reihe

${\displaystyle f(3)=1+{\frac {1}{2^{3}}}+{\frac {1}{3^{3}}}+{\frac {1}{4^{3}}}+{\frac {1}{5^{3}}}+\dotsb ,}$

da sich hier keine seiner Techniken anwenden ließ. Jedoch berechnete er die Werte ${\displaystyle f(2n+1)}$ für ${\displaystyle n=1,\dotsc ,5}$ bis auf mehrere Dezimalstellen. Außerdem schrieb er einheitlich ${\displaystyle f(n)=N\pi ^{n}}$, wobei ${\displaystyle N}$ im Falle, dass ${\displaystyle n}$ eine gerade Zahl ist, rational ist. Für den Fall, dass ${\displaystyle n}$ ungerade ist, vermutete Euler, ${\displaystyle N}$ sei „eine Funktion von ${\displaystyle \log(2)}$“.^[37] Dies konnte jedoch, ungeachtet der vagen Formulierung Eulers, bis heute nicht bestätigt werden. Die Werte der Reihen für ungerade Argumente größer als 1 sind bis heute (Stand 2020) weitestgehend unbekannt und Gegenstand zahlentheoretischer Vermutungen.^[38]

Euler gilt als Entdecker der Verbindung zwischen der Zeta-Funktion und den Primzahlen. Diese Verknüpfung wird bis heute als Euler-Produkt bezeichnet. So schrieb er in seiner Arbeit Variae observationes circa series infinitas:

Sofort war sich Euler der Beziehung zwischen den Primzahlen und der Geometrie bewusst, und er schrieb weiter:

Aus der schon damals gut bekannten Tatsache, dass die harmonische Reihe divergent ist, konnte Euler ebenfalls aus dem Euler-Produkt schließen, dass die Summe der Kehrwerte aller Primzahlen keinen endlichen Grenzwert hat.^[41] Bezeichnet wird dieses Resultat auch als Satz von Euler über die Summation der Kehrwerte der Primzahlen.

Auch war die von Riemann später bewiesene Funktionalgleichung bereits Euler bekannt. In seiner Arbeit Remarques sur un beau rapport entre les series des puissances tant directes que reciproques beschrieb er diese mathematisch nicht rigoros:^[42]

Dabei bezog sich Euler eigentlich auf die Dirichletsche Etafunktion, die jedoch bis auf einen Faktor der Riemannschen Zeta-Funktion entspricht. Euler gab keinen rigorosen Beweis einer Funktionalgleichung, sondern hatte diese nur für viele Werte überprüft und anschließend als universell gültig vermutet.

Im Jahre 1838 gelang dem Mathematiker Peter Dirichlet ein großer Beitrag zur Zahlentheorie. Er bewies eine Vermutung von Fermat, die nun Dirichletscher Primzahlsatz heißt. Diese besagt, dass jede arithmetische Progression ${\displaystyle ak+b}$ mit positiven, teilerfremden ${\displaystyle a,b}$ unendlich viele Primzahlen enthält.^[44] Ist hier zum Beispiel ${\displaystyle a=4}$ und ${\displaystyle b=1}$, folgt, dass die Liste 1, 5, 9, 13, 17 … unendlich viele Primzahlen beinhaltet.

Schlüssel zum Beweis waren neben der Riemannschen Zeta-Funktion eine ganze Klasse weiterer Funktionen, die ebenfalls in Primzahlprodukte zerfallen und somit eine „große Familie“ bilden. Erst ein Jahrhundert später konnten dank feinerer Methoden die Ergebnisse von Dirichlet durch Siegel und Walfisz deutlich präzisiert werden (Satz von Siegel-Walfisz).

Im Jahr 1859 arbeitete Bernhard Riemann in seiner Publikation Über die Anzahl der Primzahlen unter einer gegebenen Größe den bereits von Euler gegebenen Zusammenhang der Zeta-Funktion zu den Primzahlen entscheidend aus. Die große Leistung bestand darin, die Relevanz der Ausweitung des Definitionsbereichs auf komplexe Zahlen zu erkennen. Erst mit dieser Herangehensweise war es möglich geworden, konkrete Informationen über Primzahlen 2, 3, 5, 7 … selbst zu gewinnen. Das ist insofern bemerkenswert, als Primzahlen reelle Zahlen sind. Riemann, der ein Schüler von Carl Friedrich Gauß war, schrieb in seiner zehnseitigen Arbeit eine funktionentheoretische Interpretation und Auswertung des Euler-Produkts, die einen Zusammenhang zwischen Primzahlen und den nicht-trivialen Nullstellen der Zeta-Funktion schaffte. Das Hauptresultat war eine Formel, die ohne jeden Fehler die Anzahl der Primzahlen unter einer gegebenen (nicht ganzen) positiven Zahl abzählte. Damit war ihm ein völlig neuer Zugang zur Theorie der Primzahlen gelungen.

Er etablierte in seiner Arbeit das griechische ${\displaystyle \zeta }$ (Zeta) als Funktionssymbol und formulierte außerdem die bis heute unbewiesene nach ihm benannte Riemannsche Vermutung, die eine wichtige Aussage über die genaue Lage der Nullstellen der Zeta-Funktion behauptet.

Obwohl der Artikel heutzutage als Durch- und Aufbruch zur modernen analytischen Theorie um die Zeta-Funktion gesehen wird, stieß er damals in Mathematikerkreisen bei Weitem nicht nur auf Begeisterung. Schuld daran war in erster Linie, dass Riemann es an den meisten Stellen unterlassen hatte, Beweise für seine aufgestellten Formeln zu hinterlegen. So kam es, dass Godfrey Harold Hardy und John Edensor Littlewood Riemanns Arbeit lediglich als „beachtliche Ansammlung heuristischer Einsichten“ bezeichneten,^[45] die englischen Mathematiker waren allerdings in der analytischen Zahlentheorie zu Beginn des 20. Jahrhunderts anfangs noch so rückständig, dass Littlewood sich erinnerte, die Riemannsche Vermutung 1906 von seinem Professor als Übungsaufgabe gestellt bekommen zu haben. Auch Edmund Landau gehörte zu den lautesten Kritikern hinsichtlich der Bedeutung des Artikels. Zwar nannte er ihn zunächst „brillant und fruchtbar“, jedoch schlug sein Lob bald um:

Detlef Laugwitz bemerkt dazu in seiner Riemann-Biographie, dass Landau auch Eulers bahnbrechende Arbeiten in seinen Lehrbüchern wenig würdigte, da er tendenziell nur Arbeiten schätzte, in denen jedes Detail ausgearbeitet war.^[47] Auf der anderen Seite bewunderten Mathematiker wie Felix Klein, Riemann habe „mit großen allgemeinen Ideen“ gearbeitet und „oft auf seine Intuition“ vertraut.^[48] Das war noch bevor Carl Ludwig Siegel durch das Studium des Nachlasses zeigte, wie umfangreich Riemanns analytische Arbeiten zur Zeta-Funktion waren. Die Rechnungen im Nachlass waren allerdings schwer zu entschlüsseln und es bedurfte eines Mathematikers vom Kaliber Siegels, die Ideen Riemanns zu rekonstruieren.

Riemann arbeitete seit dieser Zeit bis zu seinem frühen Tode (er starb mit gerade mal 39 Jahren an den Folgen einer Tuberkulose) nicht mehr an der Zeta-Funktion, es blieb seine einzige Veröffentlichung zur Zahlentheorie.^[49] Der Aufsatz von 1859 war nur skizzenhaft ausgeführt, Riemann wollte sich damit für die Aufnahme in die Berliner Akademie der Wissenschaften bedanken.

Viele von Riemanns Aufzeichnungen wurden nach seinem Ableben von seiner Haushälterin verbrannt, bis sie von Mitarbeitern der Göttinger Fakultät gestoppt wurde. Die verbliebenen Schriften wurden seiner Witwe übergeben und verschwanden damit für viele Jahre. Über weitere Resultate zur Zeta-Funktion, die man ohne teilweise Vernichtung der Dokumente gefunden hätte, kann bis heute nur noch spekuliert werden.

Im Jahr 1893 publizierte der Mathematiker Jacques Hadamard eine Arbeit, in welcher der Grundstein eines detaillierteren Verständnisses von Riemanns Arbeit gelegt wurde. Hadamard war es gelungen, eine Formel für die Zeta-Funktion zu beweisen, die ihre Nullstellen beinhaltet. Genau genommen handelte es sich dabei um ein Verfahren, die Zeta-Funktion als Ganzes aus ihren Nullstellen zu konstruieren. Die Existenz einer solchen Formel wurde von Riemann bereits vorausgesetzt, war jedoch bis dato nicht rigoros bewiesen worden. Für die Verifikation der Riemannschen Ideen war sie aber ein substanzieller Teil: das Grundschema der Argumentation für Riemanns Hauptformel lautete nämlich „Primzahlprodukt (Euler) versus Nullstellenprodukt (Riemann/Hadamard)“. Unter anderem deshalb bezeichnete Hans von Mangoldt den Beitrag von Hadamard „als ersten wirklichen Fortschritt in diesem Gebiet seit 34 Jahren“.^[50][51]

Aufbauend auf der Arbeit von Hadamard gelang Hans von Mangoldt nur zwei Jahre später, im Jahre 1895, der Durchbruch zu Riemanns Hauptformel.^[52] Allerdings zeigte er diese in einer leicht veränderten Version, die heute als „natürlicher“^[53] angesehen wird. Seiner Leistung zu Ehren wird die Hauptformel heute als Riemann-von-Mangoldt-Formel bezeichnet.

Nachdem von Mangoldt im Jahr 1895 den Beweis der Riemannschen Hauptformel erbracht hatte, fehlte zum Beweis des Primzahlsatzes nicht mehr viel. Dieser Satz trifft eine Aussage darüber, wie häufig Primzahlen durchschnittlich auftauchen. Es blieb lediglich zu zeigen, dass die Zeta-Funktion keine Nullstellen in dem Bereich hat, in dem das Primzahlprodukt von Euler „gerade so nicht mehr gültig ist“. Unabhängig voneinander erbrachten Hadamard und der Belgier Charles-Jean de La Vallée Poussin den Beweis im Jahr 1896.^[54] Wichtige Punkte für den Beweis waren Ideen von Franz Mertens und die trigonometrische Identität ${\displaystyle 3+4\cos(\theta )+\cos(2\theta )=2(1+\cos(\theta ))^{2}}$.^[55]

Obwohl die Aufregung in der Mathematikerwelt groß war, gab es bezüglich der Beweismethode, die stark an die Eigenschaften der schwierigen Zeta-Funktion gebunden war, Bedenken über deren Natürlichkeit. Es wurde als seltsam erachtet, dass eine Aussage über Primzahlen sogar äquivalent zu einer gewissen Verteilung der Nullstellen einer komplexen Funktion war. So äußerte Albert Ingham 1932:

Im Jahr 1948 wurde schließlich ein elementarer (also gänzlich ohne funktionentheoretische Mittel auskommender) Beweis von Atle Selberg und Paul Erdös gegeben.^[58] Hierbei bedeutet „elementar“ jedoch keinesfalls „einfach“.^[59] Es wurden im Lauf der Zeit auch erheblich einfachere funktionentheoretische und elementare^[60] Beweise des Primzahlsatzes gefunden.

Im Rahmen des 2. Internationalen Mathematikerkongresses des Jahres 1900 in Paris hielt David Hilbert am 8. August einen Vortrag. In diesem formulierte er eine Liste von 23 mathematischen Problemen, die seiner Ansicht nach zu den wichtigsten des kommenden Jahrhunderts zählten. Hilbert zählte zu diesem Zeitpunkt bereits zu den führenden Mathematikern der Gegenwart. Problem Nr. 8 war die Riemannsche Vermutung:

Das hohe Ansehen, das Hilbert genoss, beflügelte die Mathematiker, sich mit seinen Problemen, darunter der Zeta-Funktion, auseinanderzusetzen. Bis heute gelten 15 der 23 Probleme als gelöst, jedoch nicht die Riemannsche Vermutung.^[62]

Im Jahre 1910 veröffentlichte der indische Mathematiker Srinivasa Ramanujan im Journal of the Indian Mathematical Society einen Artikel, in dem unter anderem die folgende Gleichung behauptet wurde:

${\displaystyle 1+2+3+4+5+\dotsb =-{\frac {1}{12}}.}$

Die meisten Mathematiker, die diese Gleichung zu Gesicht bekamen, hatten sie als offensichtlichen Schwachsinn gewertet. So kam es, dass Professor Hill vom University College in London schrieb:

Hill verhielt sich jedoch nicht völlig ablehnend und ermutigte Ramanujan, es weiter zu versuchen. Und so schickte dieser seine Ergebnisse direkt an einige Mathematiker in Cambridge. Zwei davon waren nicht in der Lage, die Aussagen hinter Ramanujans verschlüsselten Formeln zu erkennen, und lehnten die Bitte um Unterstützung ab. Als Ramanujan jedoch schließlich auch Godfrey Harold Hardy brieflich auf seine Ideen aufmerksam machte, wurde diesem in der Gleichung die korrekte Auswertung des Werts ${\displaystyle \zeta (-1)}$ bewusst, auch wenn sie bezüglich ihrer mathematischen Formalität natürlich inkorrekt war. In diesem Zusammenhang war die Formel schon Euler bekannt, da sie sich aus ${\displaystyle \zeta (2)}$ und der Funktionalgleichung ergibt (beides von Euler gefunden). Ferner erweckte Ramanujans Behauptung, eine Formel zu besitzen, die fast fehlerfrei voraussage, ob eine gegebene Zahl prim sei oder nicht (Ramanujans Formel verwendete jedoch nicht die Nullstellen der Zeta-Funktion^[63]), viel Neugier. Auch in einem zweiten Brief gab er jedoch dazu keinen Beweis. Littlewood äußerte:

Der anfänglich ausschließlich schriftliche Austausch gipfelte schließlich in einem Aufenthalt Ramanujans in England, wo sich das Duo aus Ramanujan und Hardy zu einer der produktivsten und außergewöhnlichsten mathematischen Korrespondenzen der Geschichte entwickelte.^[65]

Unter anderem nach Auswertung der Tagebücher Ramanujans durch George E. Andrews und Bruce Berndt offenbarten sich die zahlreichen Ideen Ramanujans zur Riemannschen Zeta-Funktion. So fand dieser unabhängig Eulers Formel für ${\displaystyle \zeta (2n)}$,^[66] das Euler-Produkt^[67] sowie zahlreiche unendliche Reihen und Integrale, die Zeta-Werte an ganzen und auch halbzahligen^[68] Stellen beinhalten.^[69][70][71]

Fünfzig Jahre nach Riemanns Tod tauchten einige nicht verbrannte Seiten wieder auf. Richard Dedekind, ein Kollege von Riemann hatte einige Seiten des Nachlasses von Riemanns Frau Elise erhalten und einige davon in der Bibliothek von Göttingen hinterlegt. Nachdem der Mathematikhistoriker Erich Bessel-Hagen 1926^[72] die Schriften fand und beim Versuch einer Entschlüsselung der wirren Aufzeichnungen keinen Erfolg hatte, gingen die Dokumente an Carl Ludwig Siegel. Dieser war erstaunt von der Tiefe der Gedanken Riemanns bezüglich der Zeta-Funktion.^[73] Dies entkräftete zugleich viel Kritik an Riemanns Originalarbeit, da die Notizen zeigten, dass Riemanns Behauptungen auf gründlichen Rechnungen beruhten.^[74] Jedoch beklagte auch Siegel das Chaos in den Notizen:

Siegel kam dahinter, dass Riemann ausschließlich über handschriftliche Rechnungen mindestens drei nicht-triviale Nullstellen der Zeta-Funktion relativ genau berechnet hatte. Die dafür verwendete Formel wurde durch Siegel ausgearbeitet, 1932 publiziert und wird seitdem auch Riemann-Siegel-Formel genannt.^[76]

In der Forschung rund um die Riemannsche Zeta-Funktion werden Computer vor allem dazu benutzt, die Korrektheit der Riemannschen Vermutung für möglichst viele Nullstellen zu überprüfen. Obwohl es sich bei allen Rechnungen um numerische Verfahren handelt, zeigen diese exakt und nicht nur annähernd, dass sich die untersuchten Nullstellen auf der kritischen Geraden befinden.^[77][78]

Bereits im Jahr 1936 hatte der in Oxford wirkende Mathematiker Edward Charles Titchmarsh mit einer Maschine, die ursprünglich für astronomische Berechnungen konstruiert worden war, die ersten 1.041 nicht-trivialen Nullstellen der Zeta-Funktion berechnet.^[79] Im Jahr 1953 wurden diese Berechnungen von Alan Turing fortgesetzt. Seine Methode wird bis heute benutzt. Erstmals kam dabei ein Computer zum Einsatz.^[80][81]

Ab Beginn der 80er Jahre wurden die Computer immer leistungsstärker. Bereits im Jahr 1979 hatte eine Gruppe aus Amsterdam um Herman te Riele und Richard P. Brent 200 Millionen Nullstellen überprüft (etwas später erhöhten sie ihre Rechnung auf 300 Millionen) – alle lagen auf der kritischen Geraden. Damit widersprachen sie einer Vorhersage von Don Zagier, der geäußert hatte, es sei „ein Wunder“, falls diese immer noch ausnahmslos auf der kritischen Geraden lägen.^[82] Zagier berief sich dabei auf theoretische Gründe, die zwar die Lage der ersten paar Tausend Nullstellen auf der Geraden bekräftigten, jedoch für steigende Zahlen schwächer – und letztlich sogar dagegen sprechend – auszulegen waren.

Bis 2005 wurden im Rahmen des sog. ZetaGrid Project durch verteilte Rechner die ersten 900 Milliarden Nullstellen überprüft. Um dieselbe Zeit berechnete Xavier Gourdon mit Unterstützung von Patrick Demichel die ersten 10 Billionen (${\displaystyle 10^{13}}$) Nullstellen.^[83][84] Alle lagen auf der kritischen Geraden.

Im Jahr 1972 offenbarte sich durch ein zufälliges Gespräch zwischen dem Physiker Freeman Dyson und dem Mathematiker Hugh Montgomery ein bis dahin unbemerkter Zusammenhang zwischen Quantenphysik und Zahlentheorie. Gegenstand der Diskussion waren die Nullstellen der Riemannschen Zeta-Funktion. In deren von Montgomery vermuteten Verteilung erkannte Dyson die Abstände zwischen Paaren von Eigenwerten von Hermitschen Zufallsmatrizen wieder. Diese werden von Quantenphysikern zur Vorhersage der Energieniveaus in einem schweren Atomkern verwendet, wenn dieser mit nieder-energetischen Neutronen bestrahlt wird. Als Montgomery sich die Abstände zwischen den Energieniveaus in dem Atomkern Erbium, dem 68. Element im Periodensystem, anschaute, erkannte er eine auffallende Ähnlichkeit. Die allgemeine Übereinstimmung zwischen einem gewissen Ausschnitt der Nullstellenorte auf der kritischen Geraden und den experimentell bestimmten Energieniveaus zeugte von großer Signifikanz.^[86]

Durch umfassenden Einsatz von Computern wurde die Vermutung von Montgomery über die Nullstellenabstände von Andrew Odlyzko überprüft. Die Zahlen sprachen dabei für Montgomerys Annahme. 1987 publizierte Odlyzko seine Ergebnisse.^[87]

Trotz starker Evidenzen wurden die Ergebnisse teilweise skeptisch gesehen. Es kam die Frage auf, ob damit irgendwelche Fortschritte in der reinen Mathematik erzielt worden waren. So äußerte der Zahlentheoretiker Peter Sarnak aus Princeton:

Jonathan Keating, ein Schüler des Physikers Michael Berry, lieferte schon bald eine zahlentheoretische Anwendung. Berry hatte sich zuvor schon mit Zusammenhängen zwischen Primzahlen und Quantenphysik beschäftigt (insbesondere Verbindungen zum Quantenchaos). Doch es waren schließlich Keating und seine Doktorandin Nina Snaith, die mit Hilfe statistischer Methoden (die häufig in der Quantenphysik zum Einsatz kommen) eine exakte Formel für das durchschnittliche Verhalten von Potenzen der Absolutwerte der Zeta-Funktion entlang der kritischen Geraden aufstellten. Diese Durchschnittswerte sind für die Zahlentheorie von Bedeutung und haben viele Anwendungen, zum Beispiel auf das Dirichletsche Teilerproblem. Wenige Minuten vor Keatings Vorstellung der Ergebnisse hatte er zusammen mit Snaith die Formel an einer Tafel „getestet“, indem sie schauten, ob ein bereits mühsam erarbeitetes Ergebnis korrekt vorausgesagt würde. Das Besondere an der Herangehensweise von Keating und Snaith, die unter anderem von Atle Selberg gelobt wurde^[89], war, dass sie Primzahlen dabei als Zufallsvariablen, also gewissermaßen Resultate eines