Quelltext der Seite Donald Newmans Beweis des Primzahlsatzes

'''Donald Newmans Beweis des Primzahlsatzes''' stammt aus dem Jahr 1972 und benutzt primär Mittel der [[Funktionentheorie]]. Das Manuskript von [[Donald Newman]] mit dem Titel ''Simple analytic proof of the prime number theorem'' erschien 1980 im [[American Mathematical Monthly]]. Der Beweis ist bekannt durch seine besondere Kürze und damit relativ hohe Zugänglichkeit.

Der [[Primzahlsatz]] trifft eine Aussage über die asymptotische Häufigkeitsverteilung von [[Primzahl]]en. Erstmals [[Beweis (Mathematik)|bewiesen]] wurde er im Jahr 1896, unabhängig von [[Jacques Hadamard]] und [[Charles-Jean de La Vallée Poussin]]. Deren Beweise verwendeten ebenfalls funktionentheoretische Mittel, gelten jedoch als aufwändig und überholt. Zwar wurden im Jahr 1948 sogenannte „elementare“ Beweise des Primzahlsatzes von [[Atle Selberg]] und [[Paul Erdös]] gegeben, wofür Selberg 1950 unter anderem die [[Fields-Medaille]] erhielt, jedoch bezieht sich diese Bezeichnung nur auf die ohne Funktionentheorie auskommende Methodik und nicht den Schwierigkeitsgrad.

Zum vollen Verständnis des Beweises von Newman sind lediglich Grundlagen aus der [[Analysis]] (wie [[Grenzwert (Folge)|Konvergenz]] von [[Reihe (Mathematik)|Reihen]] und [[Unendliches Produkt|Produkten]], [[Ungleichung]]en, [[Differentialrechnung|Differential-]] und [[Integralrechnung]]), der elementaren Zahlentheorie (wie [[Primzahl]]en und [[Teilbarkeit]]) und der Funktionentheorie (wie [[holomorphe Funktion]]en, die [[Cauchysche Integralformel]] und [[Weierstraßscher Konvergenzsatz]]) vonnöten. Ferner wurde er in verschiedene Standardwerke zur Funktionentheorie übernommen, etwa bei [[Theodore W. Gamelin]] und [[Serge Lang]].

== Der Primzahlsatz ==
{{Hauptartikel|Primzahlsatz}}

Bereits seit der [[Antike]] ist bekannt, dass es unendlich viele Primzahlen gibt, weshalb die Liste 2, 3, 5, 7, 11, 13, … niemals endet. Der erste Beweis dieser Aussage stammt von [[Euklid]], weshalb das Ergebnis auch [[Satz des Euklid]] genannt wird.

Die bloße Unendlichkeit einer Teilmenge der natürlichen Zahlen sagt noch nicht allzu viel über deren Natur aus. Zum Beispiel gibt es unendlich viele gerade Zahlen 2, 4, 6, 8 … und unendlich viele Quadratzahlen 1, 4, 9, 16 …, jedoch weisen beide Folgen bei genauem Hinsehen ein unterschiedliches Verhalten auf. Während zum Beispiel die Differenz zweier aufeinanderfolgender gerader Zahlen stets 2 ist, nehmen die Abstände der Quadratzahlen immer weiter zu, etwa 4 − 1 = 3, 9 − 4 = 5 und 16 − 9 = 7. Beide Folgen haben jedoch ein sehr reguläres Muster gemein, d.&nbsp;h., sie können über einfache Rechenoperationen bestimmt werden. Zum Beispiel ist die ''n''-te gerade Zahl einfach 2''n''. Im Gegensatz dazu ist bis heute kein einfaches Muster unter der Folge 2, 3, 5, 7, 11, …, 59, 61, 67, … der Primzahlen entdeckt worden. Zum Beispiel gibt es kein „schnelles“ Verfahren, die ''n''-te Primzahl zu berechnen. Die genaue Natur der Primzahlen ist bis heute ein bedeutendes offenes Problem.

Es zeigt sich jedoch, dass es ''auf lange Sicht'' Muster unter den Primzahlen zu erkennen gibt. Betrachtet man also haufenweise Primzahlen zur gleichen Zeit, so können durch „Mittelwertbildung“ reguläre Strukturen erkannt werden. Das Prinzip hinter dieser Tatsache ist von statistischer Natur. Statistik bedeutet hierbei, aus einer großen Menge von Daten Muster herauszufiltern, obwohl das „exakte Verhalten“ der einzelnen Datenobjekte (oder Subjekte) sehr kompliziert sein kann. So sind etwa alle Menschen sehr komplex, doch im Verhalten sehr vieler Menschen zur gleichen Zeit können Muster oftmals erkannt werden, die dann in Form von [[Wahrscheinlichkeit]]en auf Individuen zurück schließen lassen. Also geht es bei diesen Überlegungen zunächst um die Frage, wie die ''Verteilung'' der Primzahlen zu verstehen ist, mit anderen Worten, ''wie viele'' Primzahlen unterhalb einer vorgegebenen Schranke ''zu erwarten sind''. Zum Beispiel sind nur 4 Primzahlen, nämlich 2, 3, 5 und 7, kleiner als die Zahl 10. Im Falle der oberen Schranke 150 gibt es schon 35 kleinere Primzahlen, nämlich

:2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, <span style="color:green;">53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149</span>.
[[Datei:LogAreaDashed50150.svg|mini|hochkant=1.5|Der blau unterlegte Flächeninhalt zwischen dem Graphen der Funktion <math>\tfrac{1}{\log(t)}</math> und der ''t''-Achse im Intervall von 50 bis 150 schätzt die Anzahl der Primzahlen, die zwischen 50 und 150 liegen. Wegen der abgeflachten Kurve mit etwa 0,2 Längeneinheiten Höhe und einer Breite von 150 − 50 = 100 Längeneinheiten, schätzt „das bloße Auge“ einen Wert von 0,2 · 100 = 20 Primzahlen zwischen 50 und 150.]]
[[Datei:PiLi1000.svg|mini|hochkant=1.5|Schaubilder der Primzahl zählenden Funktion (blau) und des Integrallogarithmus (orange) im Bereich 3 bis 1000]]
Dabei sind die insgesamt 20 Primzahlen zwischen 50 und 150 in grün markiert. Eine Frage der Zahlentheorie ist, ob es ein universelles und ''einfaches'' Prinzip gibt, zumindest zu schätzen, wie viele Primzahlen es unter einer gegebenen Schranke gibt. Erkannt wurde ein solches erstmals in den Jahren 1792/93 vom damals 15-jährigen [[Carl Friedrich Gauß]],<ref>''Carl Friedrich Gauß. Werke.'' Zweiter Band. Herausgegeben von der königlichen Gesellschaft der Wissenschaften zu Göttingen, 1863, [https://books.google.de/books?id=fDMAAAAAQAAJ&pg=PA444#v=onepage&q&f=false (Brief)], S. 444–447.</ref> nachdem dieser [[Logarithmentafel]]n studiert hatte. Gauß vermutete, dass die Anzahl aller Primzahlen von 2 bis zu einer großen Zahl ''x'' ungefähr dem [[Flächeninhalt]] zwischen der ''t''-Achse und der Funktion <math>t \mapsto \tfrac{1}{\log(t)}</math> im Intervall von 2 bis <math>x</math> entspricht. Dabei ist <math>\log(t)</math> der [[Natürlicher Logarithmus|natürliche Logarithmus]]. Es gilt also die [[Integralrechnung|Integral]]-Näherung<ref name="JB41">Jörg Brüdern: ''Einführung in die analytische Zahlentheorie.'' Springer Verlag, S. 41.</ref>
:Anzahl der Primzahlen bis <math>x</math> <math>\approx \int_{2}^x \frac{1}{\log(t)} \mathrm{d}t</math>
und allgemeiner für <math>y > x \geq 2</math>:
:Anzahl der Primzahlen zwischen <math>x</math> und <math>y</math> <math>\approx \int_{x}^y \frac{1}{\log(t)} \mathrm{d}t.</math>
Zum Beispiel gilt
:<math> \int_{50}^{150} \frac{1}{\log(t)} \mathrm{d}t = 22{,}033\ldots,</math>
womit sich die Formel wegen des exakten Wertes von 20 Primzahlen zwischen 50 und 150 (siehe oben in grün) ca. um den Wert 2 verschätzt. Das Integral von <math>\tfrac{1}{\log(t)}</math> kann nicht [[Elementare Funktion|elementar]] geschlossen berechnet werden, da der Kehrwert des Logarithmus keine elementare [[Stammfunktion]] besitzt. Es definiert somit eine „eigenständige“ [[Funktion (Mathematik)|Funktion]], die auch als [[Integrallogarithmus]] <math>\operatorname{Li}</math> bekannt ist:
:<math>\operatorname{Li}(x) := \int_2^x \frac{1}{\log(t)} \mathrm{d}t.</math>
Bezeichnet <math>\pi(x)</math> die ''exakte'' Anzahl der Primzahlen unterhalb der Schranke <math>x</math>, so wird die obere Aussage <math>\pi(x) \approx \operatorname{Li}(x)</math> wie folgt präzisiert:
:<math> \lim_{x \to \infty} \frac{\pi(x)}{\operatorname{Li}(x)} = 1.</math>
Für wachsende Werte von <math>x</math> wird also der obige [[Quotient]] immer näher gegen 1 streben, also der [[Fehlerschranke#Relativer Fehler|relative Fehler]] der Schätzung <math>\pi(x) \approx \operatorname{Li}(x)</math> gegen 0 gehen. Auch bei der „Statistik der Primzahlen“ gilt demnach der Grundsatz, dass größer werdende Datenmengen prozentual eine zuverlässigere Prognose erlauben. Gauß legte keinen mathematischen Beweis für diese Vermutung über die Primzahlverteilung vor, und es dauerte noch über 100 Jahre, bis ein solcher –&nbsp;unabhängig voneinander von [[Jacques Hadamard]] und [[Charles-Jean de La Vallée Poussin]]&nbsp;– im Jahr 1896 erbracht wurde.<ref>Władysław Narkiewicz: ''The Development in Prime Number Theory.'' Springer Monographs in Mathematics, S. 183.</ref> Dabei bedeutet ''Beweis'' nicht, dass alle erdenklichen Werte durchprobiert wurden, was bei unendlich vielen Zahlen unmöglich ist, sondern dass ein auf den mathematischen Axiomen basierendes logisches Argument den Sachverhalt in voller Allgemeinheit belegt. Das damit gezeigte Theorem wird als ''Primzahlsatz'' bezeichnet.<ref name="JB41" />

== Der Beweis ==

Newmans Beweis kann in mehrere Schritte unterteilt werden. Anlässlich des 100. Beweisjahres des Primzahlsatzes wurden diese von [[Don Zagier]] in einem weiteren Artikel im [[American Mathematical Monthly]] zusammengestellt.

=== Schritt 1: Euler-Produkt der Riemannschen Zeta-Funktion ===
{{Hauptartikel|Riemannsche Zeta-Funktion|Euler-Produkt}}

Für [[komplexe Zahl]]en <math>s</math>, deren Realteil größer als 1 ist, ist die Zeta-Funktion definiert durch die Dirichlet-Reihe<ref name="FreitagBusamFunk432">E. Freitag, R. Busam: ''Funktionentheorie 1.'' Springer; 4. Aufl. 2006, ISBN 978-3-540-31764-7, S. 432.</ref>
:<math>\zeta(s) := \sum_{n=1}^\infty\frac1{n^s} = 1 + \frac1{2^s} + \frac1{3^s} + \frac1{4^s} + \frac1{5^s} + \frac1{6^s} + \frac1{7^s} + \dotsb, \, \, \, \, \, \, \, \, \, n^s := \exp(s \log(n)).</math>

Eine anschauliche Erklärung dieser Definition findet sich [[Riemannsche Zeta-Funktion#Wie „funktioniert“ die Zeta-Funktion?|hier]].

Wie man mittels des [[Integralkriterium|Integralkriteriums für unendliche Reihen]] zeigen kann, ist diese Reihe im angegebenen Bereich [[Absolute Konvergenz|absolut konvergent]]. Zudem ist die Konvergenz auf [[Kompakte Menge|kompakten]] Teilmengen gleichmäßig, weshalb nach dem [[Weierstraßscher Konvergenzsatz|Satz von Weierstraß]] die dargestellte Funktion [[Holomorphe Funktion|holomorph]] ist. Wegen der Divergenz der [[Harmonische Reihe|harmonischen Reihe]] ist diese Darstellung für alle komplexen Zahlen mit Realteil kleiner oder gleich 1 jedoch ''ungültig.'' Sie kann aber zu einer in ganz <math>\Complex \setminus \{1\}</math> holomorphen Funktion [[Riemannsche Zeta-Funktion#Methoden zur analytischen Fortsetzung|fortgesetzt werden]].

Das Euler-Produkt ist eine alternative Darstellung der Zeta-Funktion im Konvergenzbereich der Dirichlet-Reihe. Als Formel geschrieben lautet es:
:<math> \zeta(s) = \prod_{p \ \text{Primzahl}} \frac{1}{1 - \frac{1}{p^s}} = \frac{1}{1 - \frac{1}{2^s}} \cdot \frac{1}{1 - \frac{1}{3^s}} \cdot \frac{1}{1 - \frac{1}{5^s}} \cdot \frac{1}{1 - \frac{1}{7^s}} \cdot \frac{1}{1 - \frac{1}{11^s}} \cdot \frac{1}{1 - \frac{1}{13^s}} \cdot \cdots, \quad</math> wobei <math>\operatorname{Re}(s) > 1 </math>.
Dabei ist <math>\Pi</math> das [[Produkt (Mathematik)#Endliche und unendliche Produkte|Produktzeichen]], und das rechte [[Produkt (Mathematik)|Produkt]] erstreckt sich über genau alle Primzahlen. Für [[Produkt (Mathematik)#Unendliche Produkte|unendliche Produkte]] (nach Euklid gibt es unendlich viele Primzahlen) gelten ähnliche Anschauungen wie für Reihen, nur dass die Faktoren („Glieder des Produktes“) im Laufe der Zeit gegen 1 streben müssen, falls das betroffene Produkt konvergieren soll, da Faktoren nahe 1, genau wie Summanden nahe 0, nur noch wenig am Zwischenwert ändern. Das Euler-Produkt ergibt sich aus der [[Geometrische Reihe|geometrischen Reihe]] sowie dem [[Fundamentalsatz der Arithmetik]]. Es ist andersherum eine analytische Formulierung der Tatsache, dass jede natürliche Zahl <math>n</math> eine eindeutige Primfaktorzerlegung besitzt, wobei die Eindeutigkeit durch die <math>1</math> im Zähler des Terms <math>\tfrac{1}{n^s}</math> innerhalb der Zeta-Reihe ausgedrückt wird.

{| class="wikitable left mw-collapsible mw-collapsed font-size: 105.3%;"
|style="text-align:left; font-size: 95%;"| '''Für die detaillierte Herleitung'''&nbsp;&nbsp;
|-
|

Für die formale Herleitung des Euler-Produktes werden lediglich die geometrische Reihe (siehe oben), der Satz, dass jede natürliche Zahl <math>n</math> [[Fundamentalsatz der Arithmetik|''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 <math>\tfrac{1}{1-\frac{1}{p^s}}</math> als eine geometrische Reihe <math>1 + \tfrac{1}{p^s} + \left( \tfrac{1}{p^s}\right)^2 + \left( \tfrac{1}{p^s}\right)^3 + \cdots</math>, so ergibt sich im Falle nur einer Primzahl
:<math> \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</math>,
wobei das Potenzgesetz <math>\left( \tfrac{1}{p^s} \right)^n = \tfrac{1}{p^{ns}} = \tfrac{1}{(p^n)^s}</math> zu beachten ist. Zur Rechten stehen genau die Zahlen, die ausschließlich Zweien in ihrer Primfaktorzerlegung haben, also die [[Zweierpotenz]]en. Verfährt man weiter mit den ersten zwei Primzahlen, ergibt sich
:<math> \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).</math>
Multipliziert man beide Klammen aus, ergeben sich in der Summe ''alle'' Kombinationen von Termen der Form <math>\tfrac{1}{2^{ns}3^{ms}}</math> mit <math>m,n \geq 0</math>, es gilt also
:<math> \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^s3^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 </math>
und auf der rechten Seite stehen genau alle Terme <math>n^{-s}</math>, sodass <math>n</math> ''nur Zweien und Dreien'' in seiner Primfaktorzerlegung hat. Beim Ausmultiplizieren wird jeder Summand der einen Klammer mit einem Summanden der anderen Klammer verrechnet, und das in jeder Kombination, für <math>72 = 2^3 \cdot 3^2</math> sind die entsprechenden Terme in Rot markiert. Auf ähnliche Weise findet man, dass <math>\tfrac{1}{1-\frac{1}{2^s}} \cdot \tfrac{1}{1-\frac{1}{3^s}} \cdot \tfrac{1}{1-\frac{1}{5^s}}</math> zu der entsprechenden Dirichlet-Reihe korrespondiert, in der alle Zahlen mit Primfaktorzerlegung <math>2^a 3^b 5^c</math> auftauchen, und so weiter. Entsprechend gilt allgemein für die ersten <math>n</math> Primzahlen
:<math> \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}.</math>
Nun kann man in dieser Formel <math>n</math> gegen Unendlich laufen lassen, und erhält
:<math>\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),</math>
da ''jede'' Zahl <math>m</math> ''genau eine'' Zerlegung <math>m = 2^{a_1}3^{a_2}5^{a_3} \cdots </math> besitzt.
|}

=== Schritt 2: Fortsetzung der Zeta-Funktion ===
{{Siehe auch|Riemannsche Zeta-Funktion#Methoden zur analytischen Fortsetzung|titel1=Analytische Fortsetzung der Riemannschen Zeta-Funktion}}

Als Nächstes muss gezeigt werden, dass die Funktion <math>s \mapsto \zeta(s) - \tfrac{1}{s-1}</math> auf die Halbebene <math>\mathrm{Re}(s) > 0</math> fortgesetzt werden kann. Insbesondere hat <math>\zeta(s)</math> einen einfachen [[Polstelle|Pol]] in <math>s=1</math>. Dies geht zum Beispiel über die zunächst nur für <math>\mathrm{Re}(s) > 1</math> gültige Identität
:<math> \zeta(s) - \frac{1}{s-1} = \sum_{n=1}^\infty \frac{1}{n^s} - \int_1^\infty \frac{\mathrm{d}x}{x^s} = \sum_{n=1}^\infty \int_n^{n+1} \left( \frac{1}{n^s} - \frac{1}{x^s}\right) \mathrm{d}x.</math>
Über den [[Mittelwertsatz der Differentialrechnung|Mittelwertsatz]] kann man zeigen, dass die rechte Reihe sogar für <math>\mathrm{Re}(s) > 0</math> gleichmäßig und absolut konvergiert auf kompakten Teilmengen, denn
:<math> \sum_{n=1}^\infty \left| \int_n^{n+1} \left( \frac{1}{n^s} - \frac{1}{x^s}\right) \mathrm{d}x \right| = \sum_{n=1}^\infty \left| s \int_n^{n+1} \left( \int_n^x \frac{\mathrm{d}u}{u^{s+1}}\right) \mathrm{d}x \right| \leq \sum_{n=1}^\infty \max_{n \leq u \leq n+1} \left| \frac{s}{u^{s+1}}\right| = \sum_{n=1}^\infty \frac{|s|}{n^{\mathrm{Re}(s)+1}} < \infty.</math>
Damit ist die betrachtete Funktion wegen des [[Weierstraßscher Konvergenzsatz|Konvergenzsatzes von Weierstraß]] auf der Halbebene <math>\mathrm{Re}(s) > 0</math> holomorph.

=== Schritt 3: Beschränkung der Chebyshev-Funktion ===
{{Hauptartikel|Chebyshev-Funktion}}

Die erste Chebyshev-Funktion ist definiert durch
:<math> \vartheta(x) := \sum_{{p \leq x \atop p \, \text{Primzahl}}} \log(p).</math>
Es wird gezeigt, dass der [[Quotient]] <math>\tfrac{\vartheta(x)}{x}</math> für <math>x \to \infty</math> beschränkt ist. Um dies zu sehen, benutzt man zuerst den [[Binomischer Lehrsatz|binomischen Lehrsatz]]:
:<math> 2^{2n} = (1+1)^{2n} = \binom{2n}{0} + \cdots + \binom{2n}{2n} \geq \binom{2n}{n} = \frac{(2n)!}{n!^2} \geq \prod_{{n < p \leq 2n \atop p \, \text{Primzahl}}} p = e^{\vartheta(2n)-\vartheta(n)}.</math>
Ändert man in <math>\vartheta(x)</math> das Argument um einen beschränkten Term <math>O(1)</math> ab, so gilt
:<math> \vartheta(x + O(1)) - \vartheta(x) = O(\log(x)).</math>
Damit gilt mit obiger Ungleichung:
:<math> \vartheta(x) - \vartheta\left( \frac{x}{2}\right) \leq Cx</math>
mit einem <math>C > \log(2)</math> für alle <math> x \geq x_0 = x_0(C)</math>. Damit folgt durch eine [[Teleskopsumme]] mit <math>\tfrac{x}{2^R} \geq x_0 > \tfrac{x}{2^{R+1}}</math>:
:<math>\vartheta(x) = \sum_{r=0}^\infty \left( \vartheta\left( \frac{x}{2^r} \right) - \vartheta\left( \frac{x}{2^{r+1}}\right)\right) \leq \sum_{r=0}^R \left( \vartheta\left( \frac{x}{2^r} \right) - \vartheta\left( \frac{x}{2^{r+1}}\right)\right) + O(1) \leq Cx \sum_{r=0}^R \frac{1}{2^r} + O(1) = O(x).</math>

=== Schritt 4: Fortsetzbarkeit der reziproken Zeta-Funktion ===

Es ist zuerst zu zeigen, dass <math>\zeta(s) \not= 0</math> für alle <math>\mathrm{Re}(s) \geq 1</math> gilt. Damit folgt schon mit Schritt 2, dass sich die Funktion <math>s \mapsto \tfrac{1}{\zeta(s)}</math> holomorph nach <math>\mathrm{Re}(s) \geq 1</math> fortsetzen lässt. Es ist mit dem Euler-Produkt bereits <math>\zeta(s) \not= 0</math> im gesamten Bereich <math>\mathrm{Re}(s) > 1</math>. Der verbleibende Fall <math>\mathrm{Re}(s) = 1</math> lässt sich nun mit einem Satz von [[Franz Mertens (Mathematiker)|Franz Mertens]] behandeln, der zeigte, dass für <math>a_n \geq 0</math> im Bereich der Konvergenz (<math>s = \sigma + it</math> mit <math>\sigma, t \in \R</math>)
:<math> 3 \sum_{n=1}^\infty \frac{a_n}{n^\sigma} + 4 \mathrm{Re}\left(\sum_{n=1}^\infty \frac{a_n}{n^{\sigma + it}}\right) + \mathrm{Re}\left( \sum_{n=1}^\infty \frac{a_n}{n^{\sigma+2it}}\right) = \sum_{n=1}^\infty \frac{a_n (3 + 4 \cos(t\log(n)) + \cos(2t\log(n)))}{n^\sigma}\geq 0, </math>
denn es gilt <math>3 + 4 \cos(\theta) + \cos(2\theta) = 2(1 + \cos(\theta))^2 \geq 0</math>. Über das logarithmierte Euler-Produkt folgt für <math> \sigma > 1</math>
:<math> 3 \log \zeta(\sigma) + 4 \mathrm{Re}(\log \zeta(\sigma+it)) + \mathrm{Re}(\log \zeta(\sigma+2it)) \geq 0,</math>
denn
:<math> \log \zeta(s) = -\sum_{p \, \text{Primzahl}} \log\left( 1 - \frac{1}{p^s}\right) = \sum_{n=2}^\infty \frac{\Lambda(n)}{n^s \log(n)}</math>
mit der [[Mangoldt-Funktion]] <math>\Lambda</math>, und man hat <math> \tfrac{\Lambda(n)}{\log(n)} \geq 0</math>. Ergo ist
:<math> \left| e^{3 \log \zeta(\sigma) + 4 \mathrm{Re}(\log \zeta(\sigma+it)) + \mathrm{Re}(\log \zeta(\sigma+2it))} \right| = \zeta(\sigma)^3|\zeta(\sigma+it)|^4|\zeta(\sigma+2it)| \geq 1, \qquad \sigma > 1.</math>
Wäre <math>\zeta(1+it) = 0</math> für ein reelles <math>t \not= 0</math>, so wäre <math>\lim_{\sigma \to 1} \zeta(\sigma)^3|\zeta(\sigma+it)|^4|\zeta(\sigma+2it)| = 0</math> (denn der –&nbsp;siehe Schritt&nbsp;2&nbsp;– einfache Pol in <math>\sigma=1</math> wird nur dreifach gewertet), was ein Widerspruch ist.

Mit diesem Nichtverschwindungsresultat wird nun argumentiert, dass sich die Funktion <math>s \mapsto \Phi(s) - \tfrac{1}{s-1}</math> mit
:<math>\Phi(s) := \sum_{p \, \text{Primzahl}} \frac{\log(p)}{p^s} = \frac{\log(2)}{2^s} + \frac{\log(3)}{3^s} + \frac{\log(5)}{5^s} + \cdots </math>
holomorph auf die Halbebene <math>\mathrm{Re}(s) \geq 1</math> fortsetzen lässt. Dafür wird die aus der [[Logarithmische Ableitung|logarithmischen Ableitung]] des Euler-Produktes gewonnene Identität
:<math> -\frac{\zeta'(s)}{\zeta(s)} = \sum_{p \, \text{Primzahl}} \frac{\log(p)}{p^s-1} = \Phi(s) + \sum_{p \, \text{Primzahl}} \frac{\log(p)}{p^s (p^s-1)}</math>
genutzt. Die Funktion <math>s \mapsto -\tfrac{\zeta'(s)}{\zeta(s)} - \tfrac{1}{s-1}</math> dehnt sich nach dem Nichtverschwindungssatz und Schritt 2, und die Reihe <math>\textstyle s \mapsto \sum_{p \, \text{Primzahl}} \frac{\log(p)}{p^s (p^s-1)}</math> wegen gleichmäßiger Konvergenz auf Kompakta in <math>\mathrm{Re}(s) > \tfrac12</math>, holomorph auf <math>\mathrm{Re}(s) \geq 1</math> aus.

=== Schritt 5: Konvergenz des Integrals mit der Chebyshev-Funktion ===

Es wird gezeigt, dass das [[Integralrechnung|Integral]]
:<math> \int_1^\infty \frac{\vartheta(x)-x}{x^2}\mathrm{d}x = \int_0^\infty (e^{-t}\vartheta(e^t) - 1) \mathrm{d}t</math>
konvergiert. Dafür schreibt man die Funktion <math>\Phi</math> als [[Laplace-Transformation|Laplace-Transformierte]] der Chebyshev-Funktion mit exponentiellem Argument:
:<math> \Phi(s) = \sum_{p \, \text{Primzahl}} \frac{\log(p)}{p^s} = s \int_1^\infty \frac{\vartheta(x)}{x^{s+1}} \mathrm{d}x = s \int_0^\infty \vartheta(e^t)e^{-st} \mathrm{d}t.</math>
Ergo ist
:<math> \frac{\Phi(s+1)}{s+1} - \frac{1}{s} = \int_0^\infty (e^{-t}\vartheta(e^t)-1)e^{-st} \mathrm{d}t. </math>
Die Behauptung folgt nun mit Schritten 3, 4 und dem folgenden

''[[Taubersatz]] von Newman: Es sei <math>f(t)</math> mit <math>t \geq 0</math> eine beschränkte, lokal integrierbare Funktion. Es sei angenommen, die auf <math>\mathrm{Re}(s)>0</math> holomorphe Funktion''
:<math>g(s) := \int_0^\infty f(t)e^{-st}\mathrm{d}t</math>
''dehne sich holomorph auf <math>\mathrm{Re}(s) \geq 0</math> aus. Dann konvergiert <math>\textstyle \int_0^\infty f(t)\mathrm{d}t</math> und hat den Wert <math>g(0)</math>.''

=== Schritt 6: Asymptotisches Verhalten der Chebyshev-Funktion ===

Es wird <math>\vartheta(x) \sim x </math> für <math>x \to \infty </math> gezeigt. Das funktioniert zum Beispiel über einen [[Widerspruchsbeweis]]. Mit der Konvergenz des Integrals <math>\textstyle \int_0^\infty \tfrac{\vartheta(x)-x}{x^2}\mathrm{d}x</math> folgt
:<math> \lim_{T \to \infty} \int_T^{\lambda T} \frac{\vartheta(x)-x}{x^2}\mathrm{d}x = 0</math>
für alle <math>\lambda > 0</math>. Ist nun <math>\lambda > 1</math> und <math>\vartheta(T) \geq \lambda T</math> für beliebig groß werdende <math>T</math>, so folgt, da <math>\vartheta</math> nicht fallend ist,
:<math> \int_T^{\lambda T} \frac{\vartheta(x) - x}{x^2} \mathrm{d}x \geq \int_T^{\lambda T} \frac{\lambda T - x}{x^2} \mathrm{d}x = \int_{1}^\lambda \frac{\lambda-x}{x^2} \mathrm{d}x > 0</math>
und der letzte Wert hängt nicht von <math>T</math> ab, ein Widerspruch. Ähnlich wird mit <math>\lambda < 1</math> argumentiert: Ist <math>\vartheta(T) \leq \lambda T</math>, so gilt
:<math> \int_{\lambda T}^{T} \frac{\vartheta(x) - x}{x^2} \mathrm{d}x \leq \int_{\lambda T}^{T} \frac{\lambda T - x}{x^2} \mathrm{d}x = \int_\lambda^1 \frac{\lambda - x}{x^2} \mathrm{d}x < 0.</math>
Auch dies führt zu einem Widerspruch.

=== Schritt 7: Beweis des Primzahlsatzes ===

Mit Schritt 6 kann nun der Primzahlsatz bewiesen werden. Es gilt einerseits folgende Abschätzung nach oben:
:<math>\vartheta(x) = \sum_{p \leq x} \log(p) \leq \sum_{p \leq x} \log(x) = \pi(x) \log(x).</math>
Auf der anderen Seite hat man für jedes <math>\varepsilon > 0</math> folgende Abschätzung nach unten:
:<math> \vartheta(x) \geq \sum_{x^{1-\varepsilon} \leq p \leq x} \log(p) \geq \sum_{x^{1-\varepsilon} \leq p \leq x} (1-\varepsilon) \log(x) = (1-\varepsilon) \log(x)\left( \pi(x) + O(x^{1-\varepsilon}) \right).</math>
Daraus folgt <math>\vartheta(x) \sim \log(x) \pi(x) \sim x</math>, also <math>\pi(x) \sim \tfrac{x}{\log(x)}</math>.

=== Appendix: Beweis des Taubersatzes ===

Für <math>T > 0</math> wird
:<math>g_T(s) := \int_0^T f(t)e^{-st}\mathrm{d}t</math>
definiert. Dies ist holomorph für alle <math>s</math>. Es genügt, <math>\lim_{T \to \infty} g_T(0) = g(0)</math> zu zeigen. Für große <math>R > 0</math> bezeichne <math>C</math> den [[Rand (Topologie)|Rand]] des abgeschnittenen Kreises <math>\{ s \in \Complex : |s| \leq R, \mathrm{Re}(s) \geq -\delta\}</math>, wobei <math>\delta > 0</math> abhängig von <math>R</math> klein genug gewählt ist, sodass <math>g</math> holomorph in ganz <math>C</math> ist (siehe auch [[Lebesguezahl]]). Es ist dabei zu beachten, dass Holomorphie in einem (Rand-)Punkt stets in eine Umgebung dieses Punktes ausstrahlt, etwa wegen der vorhandenen [[Taylor-Entwicklung]]. Mit der [[Cauchysche Integralformel|Cauchyschen Integralformel]] folgt
:<math> g(0) - g_T(0) = \frac{1}{2\pi i} \int_C (g(s) - g_T(s))e^{sT}\left( 1 + \frac{s^2}{R^2}\right) \frac{\mathrm{d}s}{s},</math>
wobei der Term <math>\tfrac{s^2}{R^2}</math> die Funktion einer [[Komplizierte Null|komplizierten Null]] innehat und <math>e^0 = 1</math> zu beachten ist. Es wird nun argumentiert, dass dieses Integral für <math>T \to \infty</math> gegen 0 strebt. Auf dem Halbkreis <math>C_+ := C \cap \{ \mathrm{Re}(s) > 0\}</math> ist der Integrand durch <math>\tfrac{2B}{R^2}</math> beschränkt, wobei <math>B := \sup_{t \geq 0} |f(t)|</math>, denn einerseits gilt
:<math>|g(s) - g_T(s)| = \left| \int_T^\infty f(t) e^{-st} \mathrm{d}t\right| \leq B \int_{T}^\infty \left| e^{-st}\right| \mathrm{d}t = \frac{Be^{-\mathrm{Re}(s)T}}{\mathrm{Re}(s)},</math>
und andererseits
:<math> \left| e^{sT} \left( 1 + \frac{s^2}{R^2}\right)\frac{1}{s}\right| = e^{\mathrm{Re}(s)T} \cdot \frac{2\mathrm{Re}(s)}{R^2}.</math>
Nach der [[Standardabschätzung für Wegintegrale]] ist der Beitrag des Integrals über <math>C_+</math> an <math>g(0) - g_T(0)</math> durch <math>\tfrac{B}{R}</math> beschränkt. Für den Kurventeil <math>C_- := C \cap \{ \mathrm{Re}(s) < 0\}</math> werden <math>g</math> und <math>g_T</math> separat betrachtet. Da <math>g_T</math> eine [[ganze Funktion]] ist, kann der Weg <math>C_-</math> hier zum Halbkreis <math>\{ s \in \Complex : |s| = R, \mathrm{Re}(s) < 0\}</math> deformiert werden, ohne dessen Wert zu verändern, und auf dieser neuen Kurve gilt die Beschränkung durch <math>\tfrac{B}{R}</math>, denn nach gleicher Argumentation wie oben gilt
:<math>|g_T(s)| = \left| \int_0^T f(t) e^{-st}\mathrm{d}t\right| \leq B \int_{-\infty}^T |e^{-st}| \mathrm{d}t = \frac{Be^{-\mathrm{Re}(s)T}}{|\mathrm{Re}(s)|}.</math>
Zu guter Letzt geht das verbliebene Integral über <math>C_-</math> für <math>T \to \infty</math> gegen 0, denn der Integrand ist das Produkt von <math>g(s)(1 + \tfrac{s^2}{R^2})\tfrac{1}{s}</math>, was unabhängig von <math>T</math> ist, mit <math>e^{Ts}</math>, was auf kompakten Teilmengen der Halbebene <math>\{ \mathrm{Re}(s) \le 0\}</math> schnell und gleichmäßig gegen 0 strebt. Daraus folgt
:<math> \limsup_{T \to \infty} |g(0) - g_T(0)| \leq \frac{2B}{R}</math>
und mit der beliebigen Wahl <math>R \to \infty</math> die Behauptung.

== Literatur ==

* Donald Newman: ''A simple proof of the prime number theorem.'' American Mathematical Monthly, Band 87, 1980, S. 693–696.
* Don Zagier: ''Newman’s short proof of the prime number theorem.'' In: ''[[The American Mathematical Monthly]]'', Band 104, Nr. 8 (Oktober 1997), S. 705–708 [https://people.mpim-bonn.mpg.de/zagier/files/doi/10.2307/2975232/fulltext.pdf (PDF).]

== Einzelnachweise ==

<references />

[[Kategorie:Beweis (Mathematik)]]
[[Kategorie:Analytische Zahlentheorie]]