素数は 2 , 3 , 5 , 7 , 11 , … 2, 3, 5, 7, 11, \ldots 2 , 3 , 5 , 7 , 11 , … と不規則に並んでいるように見えます。ところが、x x x 以下の素数が何個あるかを正確に知ろうとすると、リーマン・ゼータ関数という一つの関数が「どこで値 0 0 0 になるか」という問いにつながります。この記事では、その関数の定義から出発して、素数の個数の誤差の大きさが零点の位置で決まる仕組みと、零点はすべて一本の直線の上にあるはずだという未解決の予想(リーマン予想)の中身を順に見ていきます。読み終えると、この予想がなぜ素数の問題として最重要とされるのかを、自分の言葉で説明できるようになります。
リーマン・ゼータ関数 ζ ( s ) = ∑ n ≥ 1 n − s \zeta(s) = \sum_{n \ge 1} n^{-s} ζ ( s ) = ∑ n ≥ 1 n − s は Re s > 1 \operatorname{Re} s > 1 Re s > 1 で絶対収束し、そこでオイラー積 ∏ p ( 1 − p − s ) − 1 \prod_p (1-p^{-s})^{-1} ∏ p ( 1 − p − s ) − 1 に等しくなります。この一つの等式が「和(すべての自然数)」と「積(素数)」を橋渡しし、ゼータ関数を素数の言語に変えます。
ζ \zeta ζ は s = 1 s = 1 s = 1 を唯一の極(単純極、留数 1 1 1 )として複素平面全体へ有理型に解析接続 され、s s s と 1 − s 1-s 1 − s を結ぶ関数等式 を満たします(有理型とは、孤立した極を除いて複素微分できる関数のこと。解析接続とは、定義域の外へ「滑らかなまま」延ばすことです)。
関数等式から s = − 2 , − 4 , − 6 , … s = -2, -4, -6, \ldots s = − 2 , − 4 , − 6 , … が零点になります(自明な零点 )。それ以外の零点は帯 0 ≤ Re s ≤ 1 0 \le \operatorname{Re} s \le 1 0 ≤ Re s ≤ 1 に閉じ込められます。
リーマン予想 とは「非自明零点の実部はすべて 1 / 2 1/2 1/2 である」という主張です。1859 年にリーマンが提示して以来、証明も反証もされていません。
重要性の核心は明示公式 にあります。素数の個数の誤差項は ∑ ρ x ρ / ρ \sum_\rho x^{\rho}/\rho ∑ ρ x ρ / ρ という零点にわたる和で書かれ、誤差の大きさは零点の実部そのものです。リーマン予想は π ( x ) = li ( x ) + O ( x log x ) \pi(x) = \operatorname{li}(x) + O(\sqrt{x}\log x) π ( x ) = li ( x ) + O ( x log x ) と同値 です。ここで O ( ⋯ ) O(\cdots) O ( ⋯ ) は「高々その程度の大きさ」を表す記号です。
この記事で使う事前知識
正則関数 : 複素数の変数の関数で、各点で複素微分できるもの。→ 正則関数
整関数 : 複素平面全体で正則な関数。→ 整関数
有理型関数 : 孤立した極を除いて正則な関数。→ 有理型関数
一様収束 : 各点ごとではなく、定義域全体で同じ速さで近づく収束。→ 一様収束
対数積分 : 素数の個数の近似に使う関数 li ( x ) \operatorname{li}(x) li ( x ) 。→ 対数積分
素数は 2 , 3 , 5 , 7 , 11 , … 2, 3, 5, 7, 11, \ldots 2 , 3 , 5 , 7 , 11 , … と、規則がありそうでない並び方をします。素数の魅力 - 素数定理 で見たように、x x x 以下の素数の個数 π ( x ) \pi(x) π ( x ) は大域的には x / log x x/\log x x / log x に近づきます(素数定理(定理 6.2)[素数の魅力と素数定理] )。しかし「近づく」だけでは足りません。数論の実際の問題では、π ( x ) \pi(x) π ( x ) と近似値の差 がどれくらい小さいかが決定的に効きます。差が x 0.9 x^{0.9} x 0.9 程度なのか x \sqrt{x} x 程度なのかで、そこから導ける結論はまるで変わります。
素数そのものを直接眺めても、この誤差の大きさは見えてきません。オイラーは 18 世紀に、まったく別の入口を見つけました。s > 1 s > 1 s > 1 の実数に対し
∑ n = 1 ∞ 1 n s = ∏ p : 素数 1 1 − p − s \sum_{n=1}^{\infty} \frac{1}{n^{s}} = \prod_{p:\ \text{素数}} \frac{1}{1 - p^{-s}} n = 1 ∑ ∞ n s 1 = p : 素数 ∏ 1 − p − s 1
が成り立つ、というのです。左辺には素数が一つも現れず、右辺は素数だけでできています。両者が等しいのは、素因数分解の一意性という算術の根本事実の言い換えにほかなりません(合同式とフェルマーの小定理 で使った一意分解、すなわち 算術の基本定理(定理 4.2)[素数の魅力と素数定理] が、ここで解析の道具に化けます)。
リーマンは 1859 年の論文で決定的な一歩を踏み出しました。s s s を複素数 に拡張し、ζ ( s ) \zeta(s) ζ ( s ) を複素平面全体で定義しなおしたのです。すると素数の個数を表す明示的な公式が得られ、その公式の中に ζ \zeta ζ の零点が現れました。素数の分布のゆらぎは、零点の位置によって完全に支配されていたのです。零点がどこにあるか——それがリーマン予想です。
この章では、ζ ( s ) \zeta(s) ζ ( s ) の定義から出発して、オイラー積・解析接続・関数等式を経て、リーマン予想の主張と、それが素数分布と同値である仕組みまでを追います。
複素数 s s s と正の実数 n n n に対し、n − s : = e − s log n n^{-s} := e^{-s \log n} n − s := e − s l o g n (log n \log n log n は実自然対数 )と定めます。s = σ + i t s = \sigma + it s = σ + i t (σ , t ∈ R \sigma, t \in \mathbb{R} σ , t ∈ R )と書くと
∣ n − s ∣ = ∣ e − ( σ + i t ) log n ∣ = e − σ log n = n − σ |n^{-s}| = |e^{-(\sigma + it)\log n}| = e^{-\sigma \log n} = n^{-\sigma} ∣ n − s ∣ = ∣ e − ( σ + i t ) l o g n ∣ = e − σ l o g n = n − σ
です。つまり絶対値は実部だけで決まります 。この一行が以後すべての収束判定の土台になります。
以下、s s s の実部を σ \sigma σ 、虚部を t t t と書く習慣に従います。また p p p は常に素数を走る変数とし、N = { 1 , 2 , 3 , … } \mathbb{N} = \{1, 2, 3, \ldots\} N = { 1 , 2 , 3 , … } (0 0 0 を含めない)とします。
定義 2.1 (リーマン・ゼータ関数(級数表示) ) ↗
複素数 s s s が Re s > 1 \operatorname{Re} s > 1 Re s > 1 を満たすとき、
ζ ( s ) : = ∑ n = 1 ∞ 1 n s \zeta(s) := \sum_{n=1}^{\infty} \frac{1}{n^{s}} ζ ( s ) := n = 1 ∑ ∞ n s 1 と定める。この級数は Re s > 1 \operatorname{Re} s > 1 Re s > 1 で絶対収束し、さらに任意の δ > 0 \delta > 0 δ > 0 に対し半平面 Re s ≥ 1 + δ \operatorname{Re} s \ge 1 + \delta Re s ≥ 1 + δ 上で一様収束する。したがって ζ \zeta ζ は Re s > 1 \operatorname{Re} s > 1 Re s > 1 で正則である。
収束を確かめておきます。σ ≥ 1 + δ \sigma \ge 1 + \delta σ ≥ 1 + δ のとき ∣ n − s ∣ = n − σ ≤ n − 1 − δ |n^{-s}| = n^{-\sigma} \le n^{-1-\delta} ∣ n − s ∣ = n − σ ≤ n − 1 − δ であり、∑ n ≥ 1 n − 1 − δ \sum_{n \ge 1} n^{-1-\delta} ∑ n ≥ 1 n − 1 − δ は積分判定法(∫ 1 ∞ x − 1 − δ d x = 1 / δ < ∞ \int_1^\infty x^{-1-\delta}dx = 1/\delta < \infty ∫ 1 ∞ x − 1 − δ d x = 1/ δ < ∞ )により収束します。よってワイヤシュトラスの M M M 判定法(各項を、収束する数列で一様に押さえられれば一様収束 する、という判定法)から一様収束が従い、正則関数列の一様収束極限は正則なので ζ \zeta ζ は正則です。一方 σ ≤ 1 \sigma \le 1 σ ≤ 1 では ∑ n − σ ≥ ∑ n − 1 = ∞ \sum n^{-\sigma} \ge \sum n^{-1} = \infty ∑ n − σ ≥ ∑ n − 1 = ∞ なので、この級数はもう使えません。σ = 1 \sigma = 1 σ = 1 が級数表示の限界です。
定理 3.1 (オイラー積表示 ) ✓
Re s > 1 \operatorname{Re} s > 1 Re s > 1 を満たすすべての複素数 s s s に対して
ζ ( s ) = ∏ p : 素数 ( 1 − p − s ) − 1 \zeta(s) = \prod_{p:\ \text{素数}} \left(1 - p^{-s}\right)^{-1} ζ ( s ) = p : 素数 ∏ ( 1 − p − s ) − 1 が成り立つ。ここで右辺は lim N → ∞ ∏ p ≤ N ( 1 − p − s ) − 1 \lim_{N \to \infty} \prod_{p \le N} (1 - p^{-s})^{-1} lim N → ∞ ∏ p ≤ N ( 1 − p − s ) − 1 の意味であり、この極限は存在する。
証明(定理 3.1) σ = Re s > 1 \sigma = \operatorname{Re} s > 1 σ = Re s > 1 を固定します。
第 1 段(各因子の展開). 任意の素数 p p p に対し ∣ p − s ∣ = p − σ ≤ 2 − σ < 1 |p^{-s}| = p^{-\sigma} \le 2^{-\sigma} < 1 ∣ p − s ∣ = p − σ ≤ 2 − σ < 1 ですから、等比級数の公式が使えて
( 1 − p − s ) − 1 = ∑ k = 0 ∞ p − k s \left(1 - p^{-s}\right)^{-1} = \sum_{k=0}^{\infty} p^{-ks} ( 1 − p − s ) − 1 = k = 0 ∑ ∞ p − k s が絶対収束します。
第 2 段(有限個の積の展開). N N N 以下の素数を p 1 < p 2 < ⋯ < p r p_1 < p_2 < \cdots < p_r p 1 < p 2 < ⋯ < p r とします。絶対収束する級数は有限個であれば自由に掛け合わせて項別に並べ替えられるので(2 個の積は ( k 1 , k 2 ) (k_1,k_2) ( k 1 , k 2 ) を添字とする絶対収束する二重級数で、定理 7.11[級数と収束判定] の 2 により和は並べ方によらない。r r r 個はこれを繰り返す)
∏ p ≤ N ( 1 − p − s ) − 1 = ∑ k 1 , … , k r ≥ 0 ( p 1 k 1 p 2 k 2 ⋯ p r k r ) − s = ∑ n ∈ A N n − s \prod_{p \le N} \left(1 - p^{-s}\right)^{-1} = \sum_{k_1, \ldots, k_r \ge 0} \left(p_1^{k_1} p_2^{k_2} \cdots p_r^{k_r}\right)^{-s} = \sum_{n \in A_N} n^{-s} p ≤ N ∏ ( 1 − p − s ) − 1 = k 1 , … , k r ≥ 0 ∑ ( p 1 k 1 p 2 k 2 ⋯ p r k r ) − s = n ∈ A N ∑ n − s となります。ここで A N A_N A N は「素因数がすべて N N N 以下であるような自然数」の集合です。二つ目の等号で素因数分解の一意性 (定理 4.2[素数の魅力と素数定理] )を使いました。A N A_N A N の各元 n n n は指数の組 ( k 1 , … , k r ) (k_1,\ldots,k_r) ( k 1 , … , k r ) とちょうど一対一に対応するので、右辺で各 n n n が過不足なく一度だけ現れます。
第 3 段(極限). n ≤ N n \le N n ≤ N ならば n n n の素因数はすべて n n n 以下、したがって N N N 以下ですから n ∈ A N n \in A_N n ∈ A N です。対偶を取れば n ∉ A N ⇒ n > N n \notin A_N \Rightarrow n > N n ∈ / A N ⇒ n > N 、すなわち A N A_N A N の補集合は { n : n > N } \{n : n > N\} { n : n > N } に含まれます。したがって
∣ ζ ( s ) − ∏ p ≤ N ( 1 − p − s ) − 1 ∣ = ∣ ∑ n ∉ A N n − s ∣ ≤ ∑ n > N n − σ . \left|\zeta(s) - \prod_{p \le N}\left(1 - p^{-s}\right)^{-1}\right| = \left|\sum_{n \notin A_N} n^{-s}\right| \le \sum_{n > N} n^{-\sigma}. ζ ( s ) − p ≤ N ∏ ( 1 − p − s ) − 1 = n ∈ / A N ∑ n − s ≤ n > N ∑ n − σ . 最後の量は収束級数 ∑ n n − σ \sum_n n^{-\sigma} ∑ n n − σ の剰余項なので N → ∞ N \to \infty N → ∞ で 0 0 0 に収束します。以上より主張が従います。
∎
この証明で使った算術の事実は素因数分解の一意性だけです。逆に言えば、オイラー積は一意分解の解析的な言い換え です。だからこそ、ζ \zeta ζ の解析的な性質を調べることが素数を調べることになるのです。
まず、この積表示から ζ \zeta ζ が右半平面で消えないことが出ます。
系 3.2 (右半平面での非零性 ) ✓
Re s > 1 \operatorname{Re} s > 1 Re s > 1 ならば ζ ( s ) ≠ 0 \zeta(s) \neq 0 ζ ( s ) = 0 である。
証明(系 3.2) σ > 1 \sigma > 1 σ > 1 とします。∑ n ∣ n − s ∣ \sum_n |n^{-s}| ∑ n ∣ n − s ∣ が収束するので並べ替えが自由にでき(定理 7.10[級数と収束判定] 、定理 7.11[級数と収束判定] )、任意の素数 p p p に対し
ζ ( s ) ( 1 − p − s ) = ∑ n ≥ 1 n − s − ∑ n ≥ 1 ( p n ) − s = ∑ p ∤ n n − s \zeta(s)\left(1 - p^{-s}\right) = \sum_{n \ge 1} n^{-s} - \sum_{n \ge 1} (pn)^{-s} = \sum_{p \,\nmid\, n} n^{-s} ζ ( s ) ( 1 − p − s ) = n ≥ 1 ∑ n − s − n ≥ 1 ∑ ( p n ) − s = p ∤ n ∑ n − s が成り立ちます(∑ n ( p n ) − s \sum_n (pn)^{-s} ∑ n ( p n ) − s はちょうど p p p の倍数にわたる和だからです)。これを p ≤ N p \le N p ≤ N なる素数について順に繰り返すと
ζ ( s ) ∏ p ≤ N ( 1 − p − s ) = ∑ n ∈ B N n − s , B N : = { n ∈ N : n の素因数はすべて N より大きい } \zeta(s) \prod_{p \le N} \left(1 - p^{-s}\right) = \sum_{n \in B_N} n^{-s}, \qquad B_N := \{n \in \mathbb{N} : n \text{ の素因数はすべて } N \text{ より大きい}\} ζ ( s ) p ≤ N ∏ ( 1 − p − s ) = n ∈ B N ∑ n − s , B N := { n ∈ N : n の素因数はすべて N より大きい } を得ます。B N B_N B N には 1 1 1 が属し、1 1 1 以外の元は N N N より大きい素因数を持つので N N N より大きい数です。よって
∣ ζ ( s ) ∏ p ≤ N ( 1 − p − s ) ∣ ≥ 1 − ∑ n > N n − σ . \left|\zeta(s)\prod_{p \le N}\left(1 - p^{-s}\right)\right| \ge 1 - \sum_{n > N} n^{-\sigma}. ζ ( s ) p ≤ N ∏ ( 1 − p − s ) ≥ 1 − n > N ∑ n − σ . ∑ n n − σ \sum_n n^{-\sigma} ∑ n n − σ が収束するので、N N N を十分大きく取れば ∑ n > N n − σ ≤ 1 / 2 \sum_{n>N} n^{-\sigma} \le 1/2 ∑ n > N n − σ ≤ 1/2 とできます。そのとき左辺は 1 / 2 1/2 1/2 以上、特に 0 0 0 ではありません。積の一因子である ζ ( s ) \zeta(s) ζ ( s ) が 0 0 0 なら左辺は 0 0 0 になってしまうので、ζ ( s ) ≠ 0 \zeta(s) \neq 0 ζ ( s ) = 0 です。
∎
証明(系 3.3) s s s を 1 1 1 より大きい実数に限って考えます。まず ζ ( s ) → ∞ \zeta(s) \to \infty ζ ( s ) → ∞ (s → 1 + 0 s \to 1+0 s → 1 + 0 )を示します。任意の N N N を固定すると、正項級数なので ζ ( s ) ≥ ∑ n = 1 N n − s \zeta(s) \ge \sum_{n=1}^{N} n^{-s} ζ ( s ) ≥ ∑ n = 1 N n − s です。右辺は s s s の連続関数で、s → 1 + 0 s \to 1+0 s → 1 + 0 のとき調和数 H N = ∑ n ≤ N 1 / n H_N = \sum_{n \le N} 1/n H N = ∑ n ≤ N 1/ n に収束します。よって lim inf s → 1 + 0 ζ ( s ) ≥ H N \liminf_{s \to 1+0} \zeta(s) \ge H_N lim inf s → 1 + 0 ζ ( s ) ≥ H N が任意の N N N で成り立ち、H N → ∞ H_N \to \infty H N → ∞ ですから ζ ( s ) → ∞ \zeta(s) \to \infty ζ ( s ) → ∞ です。
さて素数が有限個 p 1 , … , p r p_1, \ldots, p_r p 1 , … , p r しかないと仮定します。定理 3.1 の右辺は有限積になり、s → 1 + 0 s \to 1+0 s → 1 + 0 のとき
∏ i = 1 r ( 1 − p i − s ) − 1 ⟶ ∏ i = 1 r ( 1 − p i − 1 ) − 1 < ∞ \prod_{i=1}^{r} \left(1 - p_i^{-s}\right)^{-1} \longrightarrow \prod_{i=1}^{r}\left(1 - p_i^{-1}\right)^{-1} < \infty i = 1 ∏ r ( 1 − p i − s ) − 1 ⟶ i = 1 ∏ r ( 1 − p i − 1 ) − 1 < ∞ と有限値に収束します。これは ζ ( s ) → ∞ \zeta(s) \to \infty ζ ( s ) → ∞ に矛盾します。
∎
同じ道具をもう少し押すと、ユークリッドの証明(ユークリッドの定理(定理 3.2)[素数の魅力と素数定理] )では届かない「素数の密度」の情報が出ます。
命題 3.4 (素数の逆数和の発散 ) ✓
∑ p : 素数 1 p = ∞ \displaystyle\sum_{p:\ \text{素数}} \frac{1}{p} = \infty p : 素数 ∑ p 1 = ∞ である。
証明(命題 3.4) s > 1 s > 1 s > 1 を実数とします。定理 3.1 の両辺は正の実数なので対数が取れ、− log ( 1 − u ) = ∑ k ≥ 1 u k / k -\log(1-u) = \sum_{k \ge 1} u^{k}/k − log ( 1 − u ) = ∑ k ≥ 1 u k / k (0 < u < 1 0 < u < 1 0 < u < 1 )を各因子に適用して
log ζ ( s ) = ∑ p ∑ k = 1 ∞ p − k s k = ∑ p p − s + R ( s ) , R ( s ) : = ∑ p ∑ k ≥ 2 p − k s k \log \zeta(s) = \sum_{p} \sum_{k=1}^{\infty} \frac{p^{-ks}}{k} = \sum_{p} p^{-s} + R(s), \qquad R(s) := \sum_{p}\sum_{k \ge 2} \frac{p^{-ks}}{k} log ζ ( s ) = p ∑ k = 1 ∑ ∞ k p − k s = p ∑ p − s + R ( s ) , R ( s ) := p ∑ k ≥ 2 ∑ k p − k s を得ます(二重級数はすべて正項なので順序交換は自由です。定理 7.11[級数と収束判定] の 1)。剰余項 R ( s ) R(s) R ( s ) を評価します。1 / k ≤ 1 1/k \le 1 1/ k ≤ 1 と等比級数の和から
0 ≤ R ( s ) ≤ ∑ p ∑ k ≥ 2 p − k s = ∑ p p − 2 s 1 − p − s ≤ 2 ∑ p p − 2 s ≤ 2 ∑ n ≥ 2 1 n 2 = 2 ( π 2 6 − 1 ) < 1.3 0 \le R(s) \le \sum_{p} \sum_{k \ge 2} p^{-ks} = \sum_{p} \frac{p^{-2s}}{1 - p^{-s}} \le 2\sum_{p} p^{-2s} \le 2 \sum_{n \ge 2} \frac{1}{n^{2}} = 2\left(\frac{\pi^{2}}{6} - 1\right) < 1.3 0 ≤ R ( s ) ≤ p ∑ k ≥ 2 ∑ p − k s = p ∑ 1 − p − s p − 2 s ≤ 2 p ∑ p − 2 s ≤ 2 n ≥ 2 ∑ n 2 1 = 2 ( 6 π 2 − 1 ) < 1.3 です。ここで s > 1 s > 1 s > 1 , p ≥ 2 p \ge 2 p ≥ 2 より p − s ≤ 1 / 2 p^{-s} \le 1/2 p − s ≤ 1/2 、したがって ( 1 − p − s ) − 1 ≤ 2 (1-p^{-s})^{-1} \le 2 ( 1 − p − s ) − 1 ≤ 2 を使いました。よって
∑ p p − s ≥ log ζ ( s ) − 1.3. \sum_{p} p^{-s} \ge \log \zeta(s) - 1.3. p ∑ p − s ≥ log ζ ( s ) − 1.3. 系 3.3 の証明で示したとおり s → 1 + 0 s \to 1+0 s → 1 + 0 で ζ ( s ) → ∞ \zeta(s) \to \infty ζ ( s ) → ∞ 、すなわち log ζ ( s ) → ∞ \log \zeta(s) \to \infty log ζ ( s ) → ∞ ですから、∑ p p − s → ∞ \sum_p p^{-s} \to \infty ∑ p p − s → ∞ です。
もし ∑ p 1 / p \sum_p 1/p ∑ p 1/ p が収束していたとすると、s > 1 s > 1 s > 1 のとき p − s < p − 1 p^{-s} < p^{-1} p − s < p − 1 なので ∑ p p − s ≤ ∑ p p − 1 < ∞ \sum_p p^{-s} \le \sum_p p^{-1} < \infty ∑ p p − s ≤ ∑ p p − 1 < ∞ となり、s → 1 + 0 s \to 1+0 s → 1 + 0 で有界のままです。これは矛盾です。
∎
∑ n 1 / n 2 \sum_n 1/n^{2} ∑ n 1/ n 2 は収束するのに ∑ p 1 / p \sum_p 1/p ∑ p 1/ p は発散する。これは「素数は平方数よりずっと多い」という定量的な情報で、素数の無限性よりはるかに強い主張です。オイラー積が素数の分布 を掴んでいることの最初の証拠と言えます(この発散は 定理 5.2[素数の魅力と素数定理] としても扱いました)。
例 3.5 (バーゼル問題と「二数が互いに素である確率」 )
バーゼル問題の答え ζ ( 2 ) = π 2 / 6 = 1.644934 … \zeta(2) = \pi^{2}/6 = 1.644934\ldots ζ ( 2 ) = π 2 /6 = 1.644934 … を認めます。定理 3.1 を s = 2 s = 2 s = 2 で使うと
∏ p ( 1 − 1 p 2 ) = 1 ζ ( 2 ) = 6 π 2 = 0.607927 … \prod_{p}\left(1 - \frac{1}{p^{2}}\right) = \frac{1}{\zeta(2)} = \frac{6}{\pi^{2}} = 0.607927\ldots p ∏ ( 1 − p 2 1 ) = ζ ( 2 ) 1 = π 2 6 = 0.607927 … です。左辺の意味を考えます。1 ≤ a , b ≤ X 1 \le a, b \le X 1 ≤ a , b ≤ X からランダムに二つ選ぶとき、両方が p p p で割り切れる確率はおよそ 1 / p 2 1/p^{2} 1/ p 2 、割り切れない確率はおよそ 1 − 1 / p 2 1 - 1/p^{2} 1 − 1/ p 2 です。異なる素数についてこれらが独立に振る舞うと信じれば、a a a と b b b が互いに素である確率は ∏ p ( 1 − p − 2 ) = 6 / π 2 ≈ 60.79 % \prod_p (1-p^{-2}) = 6/\pi^{2} \approx 60.79\% ∏ p ( 1 − p − 2 ) = 6/ π 2 ≈ 60.79% になります。
実際、X = 100 X = 100 X = 100 で数えると、1 ≤ a , b ≤ 100 1 \le a, b \le 100 1 ≤ a , b ≤ 100 の 10000 10000 10000 組のうち gcd ( a , b ) = 1 \gcd(a,b)=1 g cd( a , b ) = 1 となるのは 6087 6087 6087 組で、60.87 % 60.87\% 60.87% です。予言との差は 0.1 0.1 0.1 ポイント未満でした。素数の「独立性」という素朴な直観が、オイラー積という等式によって正確な形を与えられています。
級数表示 定義 2.1 は Re s > 1 \operatorname{Re} s > 1 Re s > 1 でしか意味を持ちません。しかし素数の情報が濃く現れるのは、むしろその外側です。そこで解析接続 という操作が要ります。
定義 4.1 (解析接続 ) ↗
領域 U ⊂ C U \subset \mathbb{C} U ⊂ C 上の正則関数 f f f と、U U U を含む領域 V ⊃ U V \supset U V ⊃ U 上の正則関数 F F F が F ∣ U = f F|_{U} = f F ∣ U = f を満たすとき、F F F を f f f の V V V への解析接続 という。V V V が連結ならば、そのような F F F は存在すれば一意である(一致の定理)。有理型関数への拡張も同様に「有理型接続」として扱う。
一意性は重要です。「ζ ( − 1 ) = − 1 / 12 \zeta(-1) = -1/12 ζ ( − 1 ) = − 1/12 」といった等式が、勝手な定義の産物ではなく一意に決まった値 であることを保証するのがこの定理だからです。実際の接続には、Re s > 1 \operatorname{Re} s > 1 Re s > 1 で ζ \zeta ζ に一致し、より広い範囲で収束する式を一つ見つければ十分です。交代級数がその役を果たします。
補題 4.2 (ディリクレのエータ関数 ) ✓
η ( s ) : = ∑ n = 1 ∞ ( − 1 ) n − 1 n s \eta(s) := \sum_{n=1}^{\infty} \frac{(-1)^{n-1}}{n^{s}} η ( s ) := n = 1 ∑ ∞ n s ( − 1 ) n − 1 は半平面 Re s > 0 \operatorname{Re} s > 0 Re s > 0 の各コンパクト部分集合上で一様収束し、そこで正則関数を定める。
証明(補題 4.2) 第 2 n − 1 2n-1 2 n − 1 項と第 2 n 2n 2 n 項を組にして g n ( s ) : = ( 2 n − 1 ) − s − ( 2 n ) − s g_n(s) := (2n-1)^{-s} - (2n)^{-s} g n ( s ) := ( 2 n − 1 ) − s − ( 2 n ) − s とおきます。d d x ( − x − s ) = s x − s − 1 \frac{d}{dx}\left(-x^{-s}\right) = s x^{-s-1} d x d ( − x − s ) = s x − s − 1 ですから、微積分学の基本定理より
g n ( s ) = ∫ 2 n − 1 2 n s x − s − 1 d x . g_n(s) = \int_{2n-1}^{2n} s\, x^{-s-1}\,dx. g n ( s ) = ∫ 2 n − 1 2 n s x − s − 1 d x . ∣ x − s − 1 ∣ = x − σ − 1 |x^{-s-1}| = x^{-\sigma-1} ∣ x − s − 1 ∣ = x − σ − 1 なので、σ ≥ δ > 0 \sigma \ge \delta > 0 σ ≥ δ > 0 かつ ∣ s ∣ ≤ M |s| \le M ∣ s ∣ ≤ M を満たす s s s に対し
∣ g n ( s ) ∣ ≤ ∣ s ∣ ∫ 2 n − 1 2 n x − σ − 1 d x ≤ M ( 2 n − 1 ) − δ − 1 . |g_n(s)| \le |s| \int_{2n-1}^{2n} x^{-\sigma-1}dx \le M (2n-1)^{-\delta-1}. ∣ g n ( s ) ∣ ≤ ∣ s ∣ ∫ 2 n − 1 2 n x − σ − 1 d x ≤ M ( 2 n − 1 ) − δ − 1 . ∑ n ( 2 n − 1 ) − 1 − δ < ∞ \sum_n (2n-1)^{-1-\delta} < \infty ∑ n ( 2 n − 1 ) − 1 − δ < ∞ ですから、ワイヤシュトラスの M M M 判定法により ∑ n g n \sum_n g_n ∑ n g n はこの集合上で一様収束します。各 g n g_n g n は整関数なので、和は正則です。半平面 Re s > 0 \operatorname{Re} s > 0 Re s > 0 の任意のコンパクト集合はこの形の集合に含まれるので、広義一様収束が言えました。
∑ n ≤ 2 N ( − 1 ) n − 1 n − s = ∑ n ≤ N g n ( s ) \sum_{n \le 2N} (-1)^{n-1}n^{-s} = \sum_{n \le N} g_n(s) ∑ n ≤ 2 N ( − 1 ) n − 1 n − s = ∑ n ≤ N g n ( s ) ですから、η \eta η の偶数番目の部分和は収束します。奇数番目の部分和との差は ( 2 N + 1 ) − s (2N+1)^{-s} ( 2 N + 1 ) − s で、絶対値は ( 2 N + 1 ) − σ → 0 (2N+1)^{-\sigma} \to 0 ( 2 N + 1 ) − σ → 0 (σ > 0 \sigma > 0 σ > 0 )ですから、η ( s ) \eta(s) η ( s ) 自身も同じ値に収束します。
∎
定理 4.3 (ゼータ関数の有理型接続と単純極 ) ✓
ζ \zeta ζ は半平面 Re s > 0 \operatorname{Re} s > 0 Re s > 0 上の有理型関数に一意に接続され、その接続は s = 1 s = 1 s = 1 のみを極として持つ。この極は単純極で、留数は 1 1 1 である。すなわち
lim s → 1 ( s − 1 ) ζ ( s ) = 1. \lim_{s \to 1}\,(s-1)\zeta(s) = 1 . s → 1 lim ( s − 1 ) ζ ( s ) = 1.
証明(定理 4.3) 第 1 段(恒等式). Re s > 1 \operatorname{Re} s > 1 Re s > 1 では級数が絶対収束するので、偶数項だけを取り出して
∑ n 偶数 n − s = ∑ m ≥ 1 ( 2 m ) − s = 2 − s ζ ( s ) \sum_{n \text{ 偶数}} n^{-s} = \sum_{m \ge 1} (2m)^{-s} = 2^{-s}\zeta(s) n 偶数 ∑ n − s = m ≥ 1 ∑ ( 2 m ) − s = 2 − s ζ ( s ) と書けます。したがって
η ( s ) = ∑ n ≥ 1 n − s − 2 ∑ n 偶数 n − s = ζ ( s ) − 2 ⋅ 2 − s ζ ( s ) = ( 1 − 2 1 − s ) ζ ( s ) . \eta(s) = \sum_{n \ge 1} n^{-s} - 2\sum_{n \text{ 偶数}} n^{-s} = \zeta(s) - 2 \cdot 2^{-s}\zeta(s) = \left(1 - 2^{1-s}\right)\zeta(s). η ( s ) = n ≥ 1 ∑ n − s − 2 n 偶数 ∑ n − s = ζ ( s ) − 2 ⋅ 2 − s ζ ( s ) = ( 1 − 2 1 − s ) ζ ( s ) . 第 2 段(接続の定義). そこで Re s > 0 \operatorname{Re} s > 0 Re s > 0 、s ≠ 1 s \ne 1 s = 1 において
ζ ( s ) : = η ( s ) 1 − 2 1 − s \zeta(s) := \frac{\eta(s)}{1 - 2^{1-s}} ζ ( s ) := 1 − 2 1 − s η ( s ) と定めます。補題 4.2 より分子は Re s > 0 \operatorname{Re} s > 0 Re s > 0 で正則、分母は整関数ですから、右辺は分母の零点を除いて正則で、そこでは有理型です。第 1 段よりこの関数は Re s > 1 \operatorname{Re} s > 1 Re s > 1 で元の ζ \zeta ζ に一致します。一意性は 定義 4.1 の一致の定理によります。
第 3 段(s = 1 s=1 s = 1 での挙動). 分母 h ( s ) = 1 − 2 1 − s h(s) = 1 - 2^{1-s} h ( s ) = 1 − 2 1 − s は h ( 1 ) = 0 h(1) = 0 h ( 1 ) = 0 、h ′ ( s ) = 2 1 − s log 2 h'(s) = 2^{1-s}\log 2 h ′ ( s ) = 2 1 − s log 2 より h ′ ( 1 ) = log 2 ≠ 0 h'(1) = \log 2 \ne 0 h ′ ( 1 ) = log 2 = 0 なので、s = 1 s=1 s = 1 は h h h の単純零点です。一方、分子は交代調和級数の和 η ( 1 ) = log 2 ≠ 0 \eta(1) = \log 2 \ne 0 η ( 1 ) = log 2 = 0 です。よって s = 1 s=1 s = 1 は ζ \zeta ζ の単純極であり、
lim s → 1 ( s − 1 ) ζ ( s ) = lim s → 1 η ( s ) h ( s ) / ( s − 1 ) = η ( 1 ) h ′ ( 1 ) = log 2 log 2 = 1. \lim_{s \to 1}(s-1)\zeta(s) = \lim_{s\to 1}\frac{\eta(s)}{h(s)/(s-1)} = \frac{\eta(1)}{h'(1)} = \frac{\log 2}{\log 2} = 1 . s → 1 lim ( s − 1 ) ζ ( s ) = s → 1 lim h ( s ) / ( s − 1 ) η ( s ) = h ′ ( 1 ) η ( 1 ) = log 2 log 2 = 1. ∎
次の注意のねらいは、ζ ( s ) = η ( s ) / ( 1 − 2 1 − s ) \zeta(s)=\eta(s)/(1-2^{1-s}) ζ ( s ) = η ( s ) / ( 1 − 2 1 − s ) の分母が s = 1 s=1 s = 1 以外でも 0 0 0 になるため、そこにも極が立つのではないかという疑問を片づけることです。方針は、2 2 2 を 3 3 3 に替えた同じ形の式を作ると分母が消える点の集まりが別の場所になるので、二つの式を併用すれば s = 1 s=1 s = 1 以外に極がないと言える、というものです。
同じ手口をさらに進める(あるいはテータ関数のメリン変換を使う)と、ζ \zeta ζ は複素平面全体へ接続され、次の顕著な対称性が現れます。
次の関数等式にはガンマ関数 Γ \Gamma Γ (階乗 n ! n! n ! を複素数の変数へ延ばした関数で、Γ ( n ) = ( n − 1 ) ! \Gamma(n)=(n-1)! Γ ( n ) = ( n − 1 )! です)が係数として現れます。証明のねらいは、複素解析の記事で示された対称な形 Λ ( s ) = Λ ( 1 − s ) \Lambda(s)=\Lambda(1-s) Λ ( s ) = Λ ( 1 − s ) を、ガンマ関数の倍数公式と相反公式で ζ ( s ) \zeta(s) ζ ( s ) と ζ ( 1 − s ) \zeta(1-s) ζ ( 1 − s ) の比の形に書き直すことで、新しい零点や極を調べるわけではありません。
定理 4.5 (リーマンの関数等式 ) ✓
ζ \zeta ζ は C \mathbb{C} C 全体へ有理型に接続され、s = 1 s = 1 s = 1 を唯一の極(単純極、留数 1 1 1 )とする。さらに、すべての s ∈ C s \in \mathbb{C} s ∈ C に対して
ζ ( s ) = 2 s π s − 1 sin ( π s 2 ) Γ ( 1 − s ) ζ ( 1 − s ) \zeta(s) = 2^{s}\,\pi^{s-1}\,\sin\!\left(\frac{\pi s}{2}\right)\Gamma(1-s)\,\zeta(1-s) ζ ( s ) = 2 s π s − 1 sin ( 2 π s ) Γ ( 1 − s ) ζ ( 1 − s ) が有理型関数の等式として成り立つ。同値な形として、完備ゼータ関数
ξ ( s ) : = 1 2 s ( s − 1 ) π − s / 2 Γ ( s 2 ) ζ ( s ) \xi(s) := \frac{1}{2}\,s(s-1)\,\pi^{-s/2}\,\Gamma\!\left(\frac{s}{2}\right)\zeta(s) ξ ( s ) := 2 1 s ( s − 1 ) π − s /2 Γ ( 2 s ) ζ ( s ) は整関数であり、ξ ( s ) = ξ ( 1 − s ) \xi(s) = \xi(1-s) ξ ( s ) = ξ ( 1 − s ) を満たす。
証明(定理 4.5) 整関数とガンマ関数 の 定理 7.4[整関数とガンマ関数] で、ζ \zeta ζ が C \mathbb{C} C 全体へ有理型に接続され、極は s = 1 s=1 s = 1 の単純極(留数 1 1 1 )だけであること、Λ ( s ) : = π − s / 2 Γ ( s / 2 ) ζ ( s ) \Lambda(s):=\pi^{-s/2}\Gamma(s/2)\zeta(s) Λ ( s ) := π − s /2 Γ ( s /2 ) ζ ( s ) が Λ ( s ) = Λ ( 1 − s ) \Lambda(s)=\Lambda(1-s) Λ ( s ) = Λ ( 1 − s ) を満たすこと、ζ ( 0 ) = − 1 2 \zeta(0)=-\tfrac12 ζ ( 0 ) = − 2 1 、ζ ( − 2 k ) = 0 \zeta(-2k)=0 ζ ( − 2 k ) = 0 (k ≥ 1 k\ge1 k ≥ 1 )であることがテータ関数の変換則から証明されています。ここではそれをこの定理の 2 つの形に言い換えます。
ξ \xi ξ が整関数であること。 ξ ( s ) = 1 2 s ( s − 1 ) Λ ( s ) \xi(s)=\tfrac12s(s-1)\Lambda(s) ξ ( s ) = 2 1 s ( s − 1 ) Λ ( s ) です。Γ ( s / 2 ) \Gamma(s/2) Γ ( s /2 ) は s = 0 , − 2 , − 4 , … s=0,-2,-4,\dots s = 0 , − 2 , − 4 , … に単純極をもつほかは正則で(定理 6.2[整関数とガンマ関数] )、ζ \zeta ζ は s = 1 s=1 s = 1 以外で正則なので、ξ \xi ξ の特異点の候補は s = 1 s=1 s = 1 と s = 0 , − 2 , − 4 , … s=0,-2,-4,\dots s = 0 , − 2 , − 4 , … だけです。s = 1 s=1 s = 1 では ζ \zeta ζ の単純極を因子 s − 1 s-1 s − 1 が打ち消します。s = − 2 k s=-2k s = − 2 k (k ≥ 1 k\ge1 k ≥ 1 )では Γ ( s / 2 ) \Gamma(s/2) Γ ( s /2 ) の単純極を ζ \zeta ζ の零点 ζ ( − 2 k ) = 0 \zeta(-2k)=0 ζ ( − 2 k ) = 0 が打ち消します。s = 0 s=0 s = 0 では Γ ( s / 2 ) \Gamma(s/2) Γ ( s /2 ) の単純極を因子 s s s が打ち消します。いずれも除去可能特異点なので ξ \xi ξ は整関数です。s ( s − 1 ) s(s-1) s ( s − 1 ) は s ↦ 1 − s s\mapsto1-s s ↦ 1 − s で不変なので、Λ ( s ) = Λ ( 1 − s ) \Lambda(s)=\Lambda(1-s) Λ ( s ) = Λ ( 1 − s ) から ξ ( s ) = ξ ( 1 − s ) \xi(s)=\xi(1-s) ξ ( s ) = ξ ( 1 − s ) です。
非対称な形。 s ∉ Z s\notin\mathbb{Z} s ∈ / Z とします。Γ \Gamma Γ は零点をもたないので(定理 6.4[整関数とガンマ関数] )、Λ ( s ) = Λ ( 1 − s ) \Lambda(s)=\Lambda(1-s) Λ ( s ) = Λ ( 1 − s ) を書き直すと
ζ ( s ) = π s − 1 2 Γ ( 1 − s 2 ) Γ ( s 2 ) ζ ( 1 − s ) . \zeta(s)=\pi^{s-\frac12}\,\frac{\Gamma\bigl(\frac{1-s}2\bigr)}{\Gamma\bigl(\frac s2\bigr)}\,\zeta(1-s). ζ ( s ) = π s − 2 1 Γ ( 2 s ) Γ ( 2 1 − s ) ζ ( 1 − s ) . Legendre の倍数公式 Γ ( 2 ν + 1 ) = 4 ν π Γ ( ν + 1 2 ) Γ ( ν + 1 ) \Gamma(2\nu+1)=\frac{4^\nu}{\sqrt\pi}\Gamma(\nu+\frac12)\Gamma(\nu+1) Γ ( 2 ν + 1 ) = π 4 ν Γ ( ν + 2 1 ) Γ ( ν + 1 ) (補題 2.1[Bessel・Airy・Whittaker 関数] 。そこでは ν > − 1 2 \nu>-\frac12 ν > − 2 1 で示されていますが、両辺は ν \nu ν の有理型関数なので一致の定理により C \mathbb{C} C 全体で成り立ちます)に ν = − s / 2 \nu=-s/2 ν = − s /2 を代入すると Γ ( 1 − s ) = 2 − s π Γ ( 1 − s 2 ) Γ ( 1 − s 2 ) \Gamma(1-s)=\frac{2^{-s}}{\sqrt\pi}\Gamma(\frac{1-s}2)\Gamma(1-\frac s2) Γ ( 1 − s ) = π 2 − s Γ ( 2 1 − s ) Γ ( 1 − 2 s ) です。相反公式(定理 6.5[整関数とガンマ関数] )で Γ ( 1 − s 2 ) = π sin ( π s / 2 ) Γ ( s / 2 ) \Gamma(1-\frac s2)=\dfrac{\pi}{\sin(\pi s/2)\,\Gamma(s/2)} Γ ( 1 − 2 s ) = sin ( π s /2 ) Γ ( s /2 ) π と書けるので
Γ ( 1 − s 2 ) Γ ( s 2 ) = 2 s π − 1 2 sin ( π s 2 ) Γ ( 1 − s ) . \frac{\Gamma\bigl(\frac{1-s}2\bigr)}{\Gamma\bigl(\frac s2\bigr)}=2^{s}\pi^{-\frac12}\sin\Bigl(\frac{\pi s}2\Bigr)\Gamma(1-s). Γ ( 2 s ) Γ ( 2 1 − s ) = 2 s π − 2 1 sin ( 2 π s ) Γ ( 1 − s ) . 代入すると ζ ( s ) = 2 s π s − 1 sin ( π s 2 ) Γ ( 1 − s ) ζ ( 1 − s ) \zeta(s)=2^s\pi^{s-1}\sin(\frac{\pi s}2)\Gamma(1-s)\zeta(1-s) ζ ( s ) = 2 s π s − 1 sin ( 2 π s ) Γ ( 1 − s ) ζ ( 1 − s ) です。両辺は C \mathbb{C} C 上の有理型関数で、C ∖ Z \mathbb{C}\setminus\mathbb{Z} C ∖ Z で一致するので、一致の定理により有理型関数として等しくなります。
∎
ξ ( s ) = ξ ( 1 − s ) \xi(s) = \xi(1-s) ξ ( s ) = ξ ( 1 − s ) は、複素平面を直線 Re s = 1 / 2 \operatorname{Re} s = 1/2 Re s = 1/2 で折り返す鏡映対称性です。この直線が「臨界線」と呼ばれるのは、関数等式の対称軸だからです。さらに n − s ˉ = n − s ‾ n^{-\bar{s}} = \overline{n^{-s}} n − s ˉ = n − s より ζ ( s ˉ ) = ζ ( s ) ‾ \zeta(\bar{s}) = \overline{\zeta(s)} ζ ( s ˉ ) = ζ ( s ) が Re s > 1 \operatorname{Re} s > 1 Re s > 1 で成り立ち、一致の定理で全平面へ伝わります。したがって零点は実軸に関しても対称 です。零点は上下左右の二重の対称性を持って配置されているわけです。
命題 4.7 (ζ(0) の値 ) ✓
ζ ( 0 ) = − 1 2 \zeta(0) = -\frac12 ζ ( 0 ) = − 2 1 である。とくに s = 0 s=0 s = 0 は ζ \zeta ζ の零点ではない。
証明(命題 4.7) 定理 4.5 の等式で s → 0 s \to 0 s → 0 とします(ζ \zeta ζ は s = 0 s=0 s = 0 で正則なので、左辺は ζ ( 0 ) \zeta(0) ζ ( 0 ) に収束する)。sin ( π s / 2 ) ∼ π s / 2 \sin(\pi s/2) \sim \pi s/2 sin ( π s /2 ) ∼ π s /2 、また 定理 4.3 より ζ ( 1 − s ) ∼ − 1 s \zeta(1-s) \sim -\frac{1}{s} ζ ( 1 − s ) ∼ − s 1 ですから
ζ ( 0 ) = lim s → 0 2 s π s − 1 ⋅ π s 2 ⋅ Γ ( 1 − s ) ⋅ ( − 1 s ) = π − 1 ⋅ π 2 ⋅ ( − 1 ) = − 1 2 . \zeta(0) = \lim_{s\to 0} 2^{s}\pi^{s-1}\cdot\frac{\pi s}{2}\cdot \Gamma(1-s)\cdot\left(-\frac{1}{s}\right) = \pi^{-1}\cdot\frac{\pi}{2}\cdot(-1) = -\frac{1}{2}. ζ ( 0 ) = s → 0 lim 2 s π s − 1 ⋅ 2 π s ⋅ Γ ( 1 − s ) ⋅ ( − s 1 ) = π − 1 ⋅ 2 π ⋅ ( − 1 ) = − 2 1 . ∎
例 4.8 (関数等式から特殊値を計算する )
s = − 1 s = -1 s = − 1 を 定理 4.5 に代入します。sin ( − π / 2 ) = − 1 \sin(-\pi/2) = -1 sin ( − π /2 ) = − 1 、Γ ( 2 ) = 1 \Gamma(2) = 1 Γ ( 2 ) = 1 、ζ ( 2 ) = π 2 / 6 \zeta(2) = \pi^{2}/6 ζ ( 2 ) = π 2 /6 ですから
ζ ( − 1 ) = 2 − 1 π − 2 ⋅ ( − 1 ) ⋅ 1 ⋅ π 2 6 = − 1 12 . \zeta(-1) = 2^{-1}\pi^{-2}\cdot(-1)\cdot 1 \cdot \frac{\pi^{2}}{6} = -\frac{1}{12}. ζ ( − 1 ) = 2 − 1 π − 2 ⋅ ( − 1 ) ⋅ 1 ⋅ 6 π 2 = − 12 1 . s = 0 s = 0 s = 0 では sin 0 = 0 \sin 0 = 0 sin 0 = 0 と ζ ( 1 ) = ∞ \zeta(1) = \infty ζ ( 1 ) = ∞ が打ち消しあうので、極限で計算します。s → 0 s \to 0 s → 0 のとき sin ( π s / 2 ) ∼ π s / 2 \sin(\pi s/2) \sim \pi s/2 sin ( π s /2 ) ∼ π s /2 、また 定理 4.3 より ζ ( 1 − s ) ∼ 1 ( 1 − s ) − 1 = − 1 s \zeta(1-s) \sim \frac{1}{(1-s)-1} = -\frac{1}{s} ζ ( 1 − s ) ∼ ( 1 − s ) − 1 1 = − s 1 ですから
ζ ( 0 ) = lim s → 0 2 s π s − 1 ⋅ π s 2 ⋅ Γ ( 1 − s ) ⋅ ( − 1 s ) = 1 ⋅ π − 1 ⋅ π 2 ⋅ 1 ⋅ ( − 1 ) = − 1 2 . \zeta(0) = \lim_{s\to 0} 2^{s}\pi^{s-1}\cdot\frac{\pi s}{2}\cdot \Gamma(1-s)\cdot\left(-\frac{1}{s}\right) = 1 \cdot \pi^{-1}\cdot\frac{\pi}{2}\cdot 1 \cdot(-1) = -\frac{1}{2}. ζ ( 0 ) = s → 0 lim 2 s π s − 1 ⋅ 2 π s ⋅ Γ ( 1 − s ) ⋅ ( − s 1 ) = 1 ⋅ π − 1 ⋅ 2 π ⋅ 1 ⋅ ( − 1 ) = − 2 1 . 特に s = 0 s = 0 s = 0 は ζ \zeta ζ の零点ではありません 。この一点の確認が、次の命題で s = 0 s=0 s = 0 を除外する根拠になります。
命題 5.1 (自明な零点 ) ✓
k = 1 , 2 , 3 , … k = 1, 2, 3, \ldots k = 1 , 2 , 3 , … に対し ζ ( − 2 k ) = 0 \zeta(-2k) = 0 ζ ( − 2 k ) = 0 である。さらに、半平面 Re s < 0 \operatorname{Re} s < 0 Re s < 0 における ζ \zeta ζ の零点はこれらに限る。
証明(命題 5.1) σ = Re s < 0 \sigma = \operatorname{Re} s < 0 σ = Re s < 0 とし、定理 4.5 の右辺の各因子を調べます。
2 s π s − 1 2^{s}\pi^{s-1} 2 s π s − 1 は決して 0 0 0 になりません。
Re ( 1 − s ) = 1 − σ > 1 \operatorname{Re}(1-s) = 1 - \sigma > 1 Re ( 1 − s ) = 1 − σ > 1 なので、Γ ( 1 − s ) \Gamma(1-s) Γ ( 1 − s ) は極を持たず、しかもガンマ関数は零点を持たない(1 / Γ 1/\Gamma 1/Γ が整関数であることによります)ので Γ ( 1 − s ) ≠ 0 \Gamma(1-s) \neq 0 Γ ( 1 − s ) = 0 です。
同じく Re ( 1 − s ) > 1 \operatorname{Re}(1-s) > 1 Re ( 1 − s ) > 1 なので、系 3.2 より ζ ( 1 − s ) ≠ 0 \zeta(1-s) \neq 0 ζ ( 1 − s ) = 0 です。また級数表示があるので有限値です。
したがって右辺が 0 0 0 になるのは sin ( π s / 2 ) = 0 \sin(\pi s/2) = 0 sin ( π s /2 ) = 0 のとき、つまり s s s が偶数のときに限ります。σ < 0 \sigma < 0 σ < 0 の範囲でこれを満たすのは s = − 2 , − 4 , − 6 , … s = -2, -4, -6, \ldots s = − 2 , − 4 , − 6 , … です。逆にこれらの s s s では上の三つの因子が有限かつ非零で sin ( π s / 2 ) = 0 \sin(\pi s/2)=0 sin ( π s /2 ) = 0 なので、積は 0 0 0 になります。
∎
定義 5.2 (非自明零点と臨界帯 ) ↗
ζ \zeta ζ の零点のうち s = − 2 , − 4 , − 6 , … s = -2, -4, -6, \ldots s = − 2 , − 4 , − 6 , … 以外のものを非自明零点 と呼ぶ。系 3.2 (Re s > 1 \operatorname{Re} s > 1 Re s > 1 に零点なし)と 命題 5.1 (Re s < 0 \operatorname{Re} s < 0 Re s < 0 には自明零点のみ)により、非自明零点はすべて帯
0 ≤ Re s ≤ 1 0 \le \operatorname{Re} s \le 1 0 ≤ Re s ≤ 1 に属する。この帯を臨界帯 、直線 Re s = 1 / 2 \operatorname{Re} s = 1/2 Re s = 1/2 を臨界線 という。
Re Im −2 −4 −6 自明な零点 0 1 単純極 臨界帯 臨界線 Re s = 1/2 非自明零点 (実軸に関して対称)
臨界帯と零点の配置(虚部の目盛りは模式的)
公理 5.3 (リーマン予想(未解決の主張) ) ↗
ζ \zeta ζ のすべての非自明零点 ρ \rho ρ は
Re ρ = 1 2 \operatorname{Re}\rho = \frac{1}{2} Re ρ = 2 1 を満たす。同値な形で述べれば、ξ ( s ) = 0 ⇒ Re s = 1 / 2 \xi(s) = 0 \Rightarrow \operatorname{Re} s = 1/2 ξ ( s ) = 0 ⇒ Re s = 1/2 である。
例 5.5 (虚部 100 までの零点を数える )
虚部が 0 < γ < T 0 < \gamma < T 0 < γ < T の範囲にある非自明零点の個数 N ( T ) N(T) N ( T ) について、リーマン・フォン・マンゴルトの公式
N ( T ) = T 2 π log T 2 π − T 2 π + 7 8 + O ( log T ) N(T) = \frac{T}{2\pi}\log\frac{T}{2\pi} - \frac{T}{2\pi} + \frac{7}{8} + O(\log T) N ( T ) = 2 π T log 2 π T − 2 π T + 8 7 + O ( log T ) が知られています。T = 100 T = 100 T = 100 で主要項を計算すると、T / ( 2 π ) = 15.9155 T/(2\pi) = 15.9155 T / ( 2 π ) = 15.9155 、log ( 15.9155 ) = 2.7673 \log(15.9155) = 2.7673 log ( 15.9155 ) = 2.7673 より
15.9155 × 2.7673 − 15.9155 + 0.875 = 44.04 − 15.92 + 0.88 = 29.00. 15.9155 \times 2.7673 - 15.9155 + 0.875 = 44.04 - 15.92 + 0.88 = 29.00 . 15.9155 × 2.7673 − 15.9155 + 0.875 = 44.04 − 15.92 + 0.88 = 29.00. 実際に 0 < γ < 100 0 < \gamma < 100 0 < γ < 100 にある非自明零点は 29 個 で、主要項の予言とぴたりと一致します。虚部の小さいほうから並べると
γ 1 = 14.134725 … , γ 2 = 21.022040 … , γ 3 = 25.010858 … , γ 4 = 30.424876 … \gamma_1 = 14.134725\ldots,\quad \gamma_2 = 21.022040\ldots,\quad \gamma_3 = 25.010858\ldots,\quad \gamma_4 = 30.424876\ldots γ 1 = 14.134725 … , γ 2 = 21.022040 … , γ 3 = 25.010858 … , γ 4 = 30.424876 … で、対応する零点は ρ k = 1 2 + i γ k \rho_k = \tfrac{1}{2} + i\gamma_k ρ k = 2 1 + i γ k と、その複素共役 1 2 − i γ k \tfrac{1}{2} - i\gamma_k 2 1 − i γ k です。零点は虚部が大きくなるほど密に並んでおり、高さ T T T あたりの密度は 1 2 π log T 2 π \frac{1}{2\pi}\log\frac{T}{2\pi} 2 π 1 log 2 π T でゆっくり増えます。
リーマン予想が「ゼータ関数のマニアックな性質」ではなく素数の問題である理由は、明示公式 にあります。素数を数える代わりに、素数べきに log p \log p log p の重みを付けたチェビシェフ関数を使うと式が最も簡単になります。
定義 6.1 (チェビシェフ関数と対数積分 ) ↗
フォン・マンゴルト関数を
Λ ( n ) : = { log p ( n = p k , p は素数 , k ≥ 1 ) 0 ( それ以外 ) \Lambda(n) := \begin{cases} \log p & (n = p^{k},\ p \text{ は素数},\ k \ge 1) \\ 0 & (\text{それ以外}) \end{cases} Λ ( n ) := { log p 0 ( n = p k , p は素数 , k ≥ 1 ) ( それ以外 ) で定め、ψ ( x ) : = ∑ n ≤ x Λ ( n ) \psi(x) := \sum_{n \le x}\Lambda(n) ψ ( x ) := ∑ n ≤ x Λ ( n ) とおく。また対数積分(定義 6.1[素数の魅力と素数定理] )を主値積分
li ( x ) : = p . v . ∫ 0 x d t log t \operatorname{li}(x) := \mathrm{p.v.}\!\int_{0}^{x} \frac{dt}{\log t} li ( x ) := p.v. ∫ 0 x log t d t で定める。素数定理は ψ ( x ) ∼ x \psi(x) \sim x ψ ( x ) ∼ x 、あるいは π ( x ) ∼ li ( x ) \pi(x) \sim \operatorname{li}(x) π ( x ) ∼ li ( x ) と表される。
この式は驚くべきものです。左辺は素数を数える階段関数、右辺は滑らかな項 x x x と、零点ごとの波の重ね合わせです。ρ = β + i γ \rho = \beta + i\gamma ρ = β + iγ と書くと
∣ x ρ ρ ∣ = x β ∣ ρ ∣ \left|\frac{x^{\rho}}{\rho}\right| = \frac{x^{\beta}}{|\rho|} ρ x ρ = ∣ ρ ∣ x β
ですから、各零点の寄与の大きさは実部 β \beta β だけで決まります 。β \beta β が 1 1 1 に近い零点があれば誤差は x x x とほぼ同じ大きさになり、素数の分布は大きく乱れます。逆にすべての β \beta β が 1 / 2 1/2 1/2 なら、誤差はどれも x \sqrt{x} x の程度に揃います。リーマン予想とは「誤差項が可能な限り小さい」という主張なのです。
大きさの感覚をつかむために、各項の分母 ∣ ρ ∣ |\rho| ∣ ρ ∣ を無視して x β x^{\beta} x β だけを比べてみます。x = 10 6 x = 10^{6} x = 1 0 6 のとき、β = 1 / 2 \beta = 1/2 β = 1/2 なら x β = 10 3 x^{\beta} = 10^{3} x β = 1 0 3 ですが、β = 0.9 \beta = 0.9 β = 0.9 だと x β x^{\beta} x β は約 2.5 × 10 5 2.5 \times 10^{5} 2.5 × 1 0 5 になり、同じ零点一つでも誤差への効き方が二桁以上違います。
なお、素数を数える関数がなぜ複素積分と結びつくのかは次のように考えます。ψ ( x ) \psi(x) ψ ( x ) は級数 ∑ n Λ ( n ) n − s \sum_n \Lambda(n) n^{-s} ∑ n Λ ( n ) n − s の係数 Λ ( n ) \Lambda(n) Λ ( n ) を n ≤ x n \le x n ≤ x まで足したものです。この級数は − ζ ′ / ζ -\zeta'/\zeta − ζ ′ / ζ に等しく(演習 4)、級数から係数の和を取り出す道具が複素積分(ペロンの公式)です。その積分を留数定理で計算すると、− ζ ′ / ζ -\zeta'/\zeta − ζ ′ / ζ の極、つまり s = 1 s=1 s = 1 と ζ \zeta ζ の零点が拾われて、明示公式の各項になります。なお Λ ( n ) \Lambda(n) Λ ( n ) はフォン・マンゴルト関数で、n n n が素数べき p k p^k p k のときに log p \log p log p 、それ以外では 0 0 0 をとる関数です。
オイラー積: zeta(s) = prod (1 - p^-s)^-1
対数微分: -zeta'/zeta(s) = sum Lambda(n) n^-s
留数定理: s=1 の極 と 非自明零点 rho の寄与を拾う
明示公式: psi(x) = x - sum_rho x^rho / rho - ...
|x^rho| = x^(Re rho): 誤差の指数 = 零点の実部
リーマン予想 (Re rho = 1/2) は誤差 O(sqrt(x)) と同値
ゼータ関数から素数分布へ至る論理の流れ
この直観を厳密にしたのがフォン・コッホの定理です。証明は二段構えで、まず π ( x ) − li ( x ) \pi(x)-\operatorname{li}(x) π ( x ) − li ( x ) の大きさと ψ ( x ) − x \psi(x)-x ψ ( x ) − x の大きさが互いに言い換えられること(部分和の公式)を示し、そのうえで ψ \psi ψ の側で零点の位置と誤差の関係を調べます。
定理 6.3 (フォン・コッホの同値定理 ) ✓
次の二つは同値である。
リーマン予想(公理 5.3 )が成り立つ。
ある定数 C > 0 C > 0 C > 0 が存在して、すべての x ≥ 2 x \ge 2 x ≥ 2 に対し
∣ π ( x ) − li ( x ) ∣ ≤ C x log x . \left|\pi(x) - \operatorname{li}(x)\right| \le C\sqrt{x}\,\log x . ∣ π ( x ) − li ( x ) ∣ ≤ C x log x . 同値な別形として、リーマン予想は「任意の ε > 0 \varepsilon > 0 ε > 0 に対し ψ ( x ) = x + O ( x 1 / 2 + ε ) \psi(x) = x + O\!\left(x^{1/2+\varepsilon}\right) ψ ( x ) = x + O ( x 1/2 + ε ) 」とも同値である。
証明(定理 6.3) 補題は Appendix「フォン・コッホの同値定理の証明」に置きます。θ ( x ) = ∑ p ≤ x log p \theta(x)=\sum_{p\le x}\log p θ ( x ) = ∑ p ≤ x log p とし、x ≥ 2 x\ge2 x ≥ 2 で li ( x ) = li ( 2 ) + ∫ 2 x d t log t \operatorname{li}(x)=\operatorname{li}(2)+\int_2^x\frac{dt}{\log t} li ( x ) = li ( 2 ) + ∫ 2 x l o g t d t であることを使います(主値は t = 1 t=1 t = 1 でだけ必要)。
部分和の公式. x ≥ 2 x\ge2 x ≥ 2 について、log p = log x − ∫ p x d t t \log p=\log x-\int_p^x\frac{dt}t log p = log x − ∫ p x t d t を素数 p ≤ x p\le x p ≤ x について足し、Tonelli の定理で和と積分を交換すると θ ( x ) = π ( x ) log x − ∫ 2 x π ( t ) t d t \theta(x)=\pi(x)\log x-\int_2^x\frac{\pi(t)}t\,dt θ ( x ) = π ( x ) log x − ∫ 2 x t π ( t ) d t 。同様に 1 log p = 1 log x + ∫ p x d t t log 2 t \frac1{\log p}=\frac1{\log x}+\int_p^x\frac{dt}{t\log^2t} l o g p 1 = l o g x 1 + ∫ p x t l o g 2 t d t を log p \log p log p を掛けて足すと π ( x ) = θ ( x ) log x + ∫ 2 x θ ( t ) t log 2 t d t \pi(x)=\frac{\theta(x)}{\log x}+\int_2^x\frac{\theta(t)}{t\log^2t}dt π ( x ) = l o g x θ ( x ) + ∫ 2 x t l o g 2 t θ ( t ) d t 。一方 d d x ( li ( x ) log x ) = 1 + li ( x ) x \frac{d}{dx}(\operatorname{li}(x)\log x)=1+\frac{\operatorname{li}(x)}x d x d ( li ( x ) log x ) = 1 + x li ( x ) 、∫ 2 x d t log t = x log x − 2 log 2 + ∫ 2 x d t log 2 t \int_2^x\frac{dt}{\log t}=\frac x{\log x}-\frac2{\log2}+\int_2^x\frac{dt}{\log^2t} ∫ 2 x l o g t d t = l o g x x − l o g 2 2 + ∫ 2 x l o g 2 t d t (部分積分)なので
θ ( x ) − x = ( π ( x ) − li ( x ) ) log x − ∫ 2 x π ( t ) − li ( t ) t d t + c 1 , π ( x ) − li ( x ) = θ ( x ) − x log x + ∫ 2 x θ ( t ) − t t log 2 t d t + c 2 \theta(x)-x=\bigl(\pi(x)-\operatorname{li}(x)\bigr)\log x-\int_2^x\frac{\pi(t)-\operatorname{li}(t)}{t}\,dt+c_1,\qquad
\pi(x)-\operatorname{li}(x)=\frac{\theta(x)-x}{\log x}+\int_2^x\frac{\theta(t)-t}{t\log^2t}\,dt+c_2 θ ( x ) − x = ( π ( x ) − li ( x ) ) log x − ∫ 2 x t π ( t ) − li ( t ) d t + c 1 , π ( x ) − li ( x ) = log x θ ( x ) − x + ∫ 2 x t log 2 t θ ( t ) − t d t + c 2 です(c 1 = li ( 2 ) log 2 − 2 c_1=\operatorname{li}(2)\log2-2 c 1 = li ( 2 ) log 2 − 2 、c 2 = 2 log 2 − li ( 2 ) c_2=\frac2{\log2}-\operatorname{li}(2) c 2 = l o g 2 2 − li ( 2 ) )。
( 1 ) ⇒ ( 2 ) (1)\Rightarrow(2) ( 1 ) ⇒ ( 2 ) . 補題 7.11 と 補題 7.5 の 0 ≤ ψ ( x ) − θ ( x ) ≤ 2 x log x 0\le\psi(x)-\theta(x)\le2\sqrt x\log x 0 ≤ ψ ( x ) − θ ( x ) ≤ 2 x log x から、x ≥ 2 x\ge2 x ≥ 2 で ∣ θ ( x ) − x ∣ ≤ C x log 2 x |\theta(x)-x|\le C\sqrt x\log^2x ∣ θ ( x ) − x ∣ ≤ C x log 2 x です。第 2 の公式に代入すると
∣ π ( x ) − li ( x ) ∣ ≤ C x log x + C ∫ 2 x d t t + ∣ c 2 ∣ ≤ C ′ x log x ( x ≥ 2 ) . |\pi(x)-\operatorname{li}(x)|\le C\sqrt x\log x+C\int_2^x\frac{dt}{\sqrt t}+|c_2|\le C'\sqrt x\log x\qquad(x\ge2). ∣ π ( x ) − li ( x ) ∣ ≤ C x log x + C ∫ 2 x t d t + ∣ c 2 ∣ ≤ C ′ x log x ( x ≥ 2 ) . ( 2 ) ⇒ ( 1 ) (2)\Rightarrow(1) ( 2 ) ⇒ ( 1 ) . 第 1 の公式に ∣ π ( t ) − li ( t ) ∣ ≤ C t log t |\pi(t)-\operatorname{li}(t)|\le C\sqrt t\log t ∣ π ( t ) − li ( t ) ∣ ≤ C t log t を代入すると ∣ θ ( x ) − x ∣ ≤ C x log 2 x + C ∫ 2 x log t t d t + ∣ c 1 ∣ ≤ C ′ ′ x log 2 x |\theta(x)-x|\le C\sqrt x\log^2x+C\int_2^x\frac{\log t}{\sqrt t}dt+|c_1|\le C''\sqrt x\log^2x ∣ θ ( x ) − x ∣ ≤ C x log 2 x + C ∫ 2 x t l o g t d t + ∣ c 1 ∣ ≤ C ′′ x log 2 x 、したがって ∣ ψ ( x ) − x ∣ ≤ C ′ ′ ′ x log 2 x |\psi(x)-x|\le C'''\sqrt x\log^2x ∣ ψ ( x ) − x ∣ ≤ C ′′′ x log 2 x です。任意の ε > 0 \varepsilon\gt0 ε > 0 で log 2 x ≤ C ε x ε \log^2x\le C_\varepsilon x^\varepsilon log 2 x ≤ C ε x ε なので、補題 7.12 を α = 1 2 \alpha=\frac12 α = 2 1 で使うとリーマン予想(公理 5.3 )が従います。
別形. リーマン予想の下では、上のとおり ψ ( x ) = x + O ( x log 2 x ) = x + O ( x 1 / 2 + ε ) \psi(x)=x+O(\sqrt x\log^2x)=x+O(x^{1/2+\varepsilon}) ψ ( x ) = x + O ( x log 2 x ) = x + O ( x 1/2 + ε ) 。逆に、すべての ε > 0 \varepsilon\gt0 ε > 0 で ψ ( x ) = x + O ( x 1 / 2 + ε ) \psi(x)=x+O(x^{1/2+\varepsilon}) ψ ( x ) = x + O ( x 1/2 + ε ) なら、補題 7.12 からリーマン予想が成り立ちます。
∎
例 6.5 (素数計数関数と対数積分の差を実際に見る )
既知の値を並べます。
x x x π ( x ) \pi(x) π ( x ) li ( x ) \operatorname{li}(x) li ( x ) (概算)li ( x ) − π ( x ) \operatorname{li}(x)-\pi(x) li ( x ) − π ( x ) x \sqrt{x} x 10 4 10^{4} 1 0 4 1 229 1\,229 1 229 1 246.1 1\,246.1 1 246.1 17.1 17.1 17.1 100 100 100 10 6 10^{6} 1 0 6 78 498 78\,498 78 498 78 627.5 78\,627.5 78 627.5 129.5 129.5 129.5 1 000 1\,000 1 000 10 8 10^{8} 1 0 8 5 761 455 5\,761\,455 5 761 455 5 762 209.4 5\,762\,209.4 5 762 209.4 754.4 754.4 754.4 10 4 10^{4} 1 0 4 10 10 10^{10} 1 0 10 455 052 511 455\,052\,511 455 052 511 455 055 614.6 455\,055\,614.6 455 055 614.6 3 103.6 3\,103.6 3 103.6 10 5 10^{5} 1 0 5 10 12 10^{12} 1 0 12 37 607 912 018 37\,607\,912\,018 37 607 912 018 37 607 950 280.8 37\,607\,950\,280.8 37 607 950 280.8 38 262.8 38\,262.8 38 262.8 10 6 10^{6} 1 0 6
差はどの行でも x \sqrt{x} x より小さく、x x x が 10 8 10^{8} 1 0 8 倍になっても差は 2000 2000 2000 倍程度にしか増えていません。x \sqrt{x} x は 10 4 10^{4} 1 0 4 倍になっているので、経験的には x \sqrt{x} x よりむしろ緩やかです。定理 6.3 の主張する x log x \sqrt{x}\log x x log x の枠に、余裕をもって収まっています。
表を見ると li ( x ) > π ( x ) \operatorname{li}(x) > \pi(x) li ( x ) > π ( x ) が常に成り立ちそうに見えますが、これは誤りです。リトルウッド(1914)は li ( x ) − π ( x ) \operatorname{li}(x) - \pi(x) li ( x ) − π ( x ) が符号を無限回変える ことを証明しました。最初の符号変化が起きる x x x の位置は今も特定されておらず、上界(スキューズ数と呼ばれる系列の評価)が知られているだけです。有限の数値実験がどれほど説得的に見えても証明にはならない、という数論の教訓がここにあります。
10⁴ 10⁶ 10⁸ 10¹⁰ 10¹² x 10 10² 10³ 10⁴ 10⁵ 10⁶ 破線: √x 実線: li(x) − π(x) li(x) − π(x) と √x の比較(両軸とも対数目盛り。値は上の表のもの)
リーマン予想が素数以外にも波及することを示すため、同値な言明をいくつか挙げます。見た目の分野がまったく違うのに、すべて臨界線上の零点という一点に帰着します。
同値な言明 内容 出典 素数計数の誤差 π ( x ) = li ( x ) + O ( x log x ) \pi(x) = \operatorname{li}(x) + O(\sqrt{x}\log x) π ( x ) = li ( x ) + O ( x log x ) von Koch (1901) チェビシェフ関数 任意の ε > 0 \varepsilon > 0 ε > 0 で ψ ( x ) = x + O ( x 1 / 2 + ε ) \psi(x) = x + O(x^{1/2+\varepsilon}) ψ ( x ) = x + O ( x 1/2 + ε ) 明示公式から メルテンス関数 任意の ε > 0 \varepsilon > 0 ε > 0 で ∑ n ≤ x μ ( n ) = O ( x 1 / 2 + ε ) \sum_{n \le x}\mu(n) = O(x^{1/2+\varepsilon}) ∑ n ≤ x μ ( n ) = O ( x 1/2 + ε ) 古典的 約数和の上界 n > 5040 n > 5040 n > 5040 なるすべての n n n で σ ( n ) < e γ n log log n \sigma(n) < e^{\gamma} n \log\log n σ ( n ) < e γ n log log n Robin (1984) 調和数による判定 すべての n ≥ 1 n \ge 1 n ≥ 1 で σ ( n ) ≤ H n + e H n log H n \sigma(n) \le H_n + e^{H_n}\log H_n σ ( n ) ≤ H n + e H n log H n (等号は n = 1 n=1 n = 1 のみ) Lagarias (2002)
ここで μ \mu μ はメビウス関数(n n n が 1 1 1 より大きい平方数で割り切れるなら 0 0 0 、異なる k k k 個の素数の積なら ( − 1 ) k (-1)^k ( − 1 ) k をとる関数)、σ ( n ) \sigma(n) σ ( n ) は n n n の約数の総和 、H n = ∑ k ≤ n 1 / k H_n = \sum_{k \le n} 1/k H n = ∑ k ≤ n 1/ k 、γ \gamma γ はオイラー・マスケローニ定数です。最後の二つは複素解析をまったく含まない初等的な形をしていますが、リーマン予想と完全に同値です。とくに Lagarias の判定は、リーマン予想を算術の Π 1 \Pi_1 Π 1 文として書き下すための道具にもなります(例 8.2[ゲーデルの不完全性定理] )。
さらに、この構図は ζ \zeta ζ に限りません。数論に現れる多くの対象——ディリクレの L L L 関数、代数体のデデキント・ゼータ関数、そして 楕円曲線とモジュラー形式 で扱う楕円曲線の L L L 関数(定義 4.4[楕円曲線・モジュラー形式と現代数論] のトレース a p a_p a p から作られます)——にも同型の関数等式と「一般化リーマン予想 」があります。同じ形の予想が広い範囲で立つこと自体が、零点の位置が偶然ではないことを示唆しています。リーマン予想は、素数分布という一つの問題であると同時に、L L L 関数という一族全体の設計図についての問いでもあるのです。
演習 7.1 易
Re s > 1 \operatorname{Re} s > 1 Re s > 1 のとき ( 1 − 2 1 − s ) ζ ( s ) = ∑ n ≥ 1 ( − 1 ) n − 1 n − s \left(1 - 2^{1-s}\right)\zeta(s) = \sum_{n\ge1}(-1)^{n-1}n^{-s} ( 1 − 2 1 − s ) ζ ( s ) = ∑ n ≥ 1 ( − 1 ) n − 1 n − s が成り立つことを、絶対収束を根拠として示してください。また s = 2 s = 2 s = 2 でこの等式を確かめ、∑ n ≥ 1 ( − 1 ) n − 1 / n 2 \sum_{n \ge 1}(-1)^{n-1}/n^{2} ∑ n ≥ 1 ( − 1 ) n − 1 / n 2 の値を求めてください。
解答 Re s = σ > 1 \operatorname{Re} s = \sigma > 1 Re s = σ > 1 で ∑ n ∣ n − s ∣ = ∑ n n − σ < ∞ \sum_n |n^{-s}| = \sum_n n^{-\sigma} < \infty ∑ n ∣ n − s ∣ = ∑ n n − σ < ∞ なので、級数の項を自由に並べ替え・分割できます。偶数 n = 2 m n = 2m n = 2 m にわたる和は
∑ m ≥ 1 ( 2 m ) − s = 2 − s ∑ m ≥ 1 m − s = 2 − s ζ ( s ) \sum_{m \ge 1}(2m)^{-s} = 2^{-s}\sum_{m\ge1} m^{-s} = 2^{-s}\zeta(s) m ≥ 1 ∑ ( 2 m ) − s = 2 − s m ≥ 1 ∑ m − s = 2 − s ζ ( s ) です。したがって
∑ n ≥ 1 ( − 1 ) n − 1 n s = ∑ n ≥ 1 n − s − 2 ∑ n 偶数 n − s = ζ ( s ) − 2 ⋅ 2 − s ζ ( s ) = ( 1 − 2 1 − s ) ζ ( s ) . \sum_{n\ge1}\frac{(-1)^{n-1}}{n^{s}} = \sum_{n \ge 1}n^{-s} - 2\sum_{n \text{ 偶数}} n^{-s} = \zeta(s) - 2\cdot 2^{-s}\zeta(s) = \left(1-2^{1-s}\right)\zeta(s). n ≥ 1 ∑ n s ( − 1 ) n − 1 = n ≥ 1 ∑ n − s − 2 n 偶数 ∑ n − s = ζ ( s ) − 2 ⋅ 2 − s ζ ( s ) = ( 1 − 2 1 − s ) ζ ( s ) . s = 2 s = 2 s = 2 とすると 1 − 2 − 1 = 1 / 2 1 - 2^{-1} = 1/2 1 − 2 − 1 = 1/2 なので
∑ n ≥ 1 ( − 1 ) n − 1 n 2 = 1 2 ζ ( 2 ) = 1 2 ⋅ π 2 6 = π 2 12 = 0.822467 … \sum_{n\ge1}\frac{(-1)^{n-1}}{n^{2}} = \frac{1}{2}\zeta(2) = \frac{1}{2}\cdot\frac{\pi^{2}}{6} = \frac{\pi^{2}}{12} = 0.822467\ldots n ≥ 1 ∑ n 2 ( − 1 ) n − 1 = 2 1 ζ ( 2 ) = 2 1 ⋅ 6 π 2 = 12 π 2 = 0.822467 … です。
演習 7.2 標準
μ \mu μ をメビウス関数とします。μ ( n ) 2 \mu(n)^{2} μ ( n ) 2 は n n n が平方因子を持たないとき 1 1 1 、そうでないとき 0 0 0 です。Re s > 1 \operatorname{Re} s > 1 Re s > 1 において
∑ n = 1 ∞ μ ( n ) 2 n s = ζ ( s ) ζ ( 2 s ) \sum_{n=1}^{\infty}\frac{\mu(n)^{2}}{n^{s}} = \frac{\zeta(s)}{\zeta(2s)} n = 1 ∑ ∞ n s μ ( n ) 2 = ζ ( 2 s ) ζ ( s ) を示してください。
解答 σ > 1 \sigma > 1 σ > 1 とします。定理 3.1 の証明と同じ論法を使います。n n n が平方因子を持たないことは、n n n の素因数分解に各素数が高々 1 1 1 回しか現れないことと同値です。したがって有限積を展開すると
∏ p ≤ N ( 1 + p − s ) = ∑ n ∈ A N , μ ( n ) 2 = 1 n − s \prod_{p \le N}\left(1 + p^{-s}\right) = \sum_{n \in A_N,\ \mu(n)^2 = 1} n^{-s} p ≤ N ∏ ( 1 + p − s ) = n ∈ A N , μ ( n ) 2 = 1 ∑ n − s となります(各素数について「使わない」か「1 1 1 回だけ使う」かの二択で、素因数分解の一意性からこれらの n n n が重複なく現れます)。ここで A N A_N A N は素因数がすべて N N N 以下の自然数の集合です。n ≤ N n \le N n ≤ N なら必ず n ∈ A N n \in A_N n ∈ A N ですから、左辺と無限和の差は n > N n > N n > N の項だけからなり、その絶対値は ∑ n > N n − σ → 0 \sum_{n>N} n^{-\sigma} \to 0 ∑ n > N n − σ → 0 で抑えられます。よって N → ∞ N \to \infty N → ∞ として
∑ n ≥ 1 μ ( n ) 2 n s = ∏ p ( 1 + p − s ) . \sum_{n\ge1}\frac{\mu(n)^{2}}{n^{s}} = \prod_{p}\left(1+p^{-s}\right). n ≥ 1 ∑ n s μ ( n ) 2 = p ∏ ( 1 + p − s ) . (∑ n μ ( n ) 2 n − σ ≤ ζ ( σ ) < ∞ \sum_n \mu(n)^2 n^{-\sigma} \le \zeta(\sigma) < \infty ∑ n μ ( n ) 2 n − σ ≤ ζ ( σ ) < ∞ なので絶対収束は保証されます。)
次に恒等式 1 + u = 1 − u 2 1 − u 1 + u = \dfrac{1-u^{2}}{1-u} 1 + u = 1 − u 1 − u 2 (u = p − s u = p^{-s} u = p − s )を使うと
∏ p ( 1 + p − s ) = ∏ p 1 − p − 2 s 1 − p − s = ∏ p ( 1 − p − 2 s ) ∏ p ( 1 − p − s ) = ζ ( s ) ζ ( 2 s ) . \prod_p \left(1+p^{-s}\right) = \prod_p \frac{1 - p^{-2s}}{1-p^{-s}} = \frac{\prod_p\left(1-p^{-2s}\right)}{\prod_p\left(1-p^{-s}\right)} = \frac{\zeta(s)}{\zeta(2s)} . p ∏ ( 1 + p − s ) = p ∏ 1 − p − s 1 − p − 2 s = ∏ p ( 1 − p − s ) ∏ p ( 1 − p − 2 s ) = ζ ( 2 s ) ζ ( s ) . 最後の等号で 定理 3.1 を s s s と 2 s 2s 2 s の両方に適用しました(Re ( 2 s ) > 2 > 1 \operatorname{Re}(2s) > 2 > 1 Re ( 2 s ) > 2 > 1 なので適用可能で、系 3.2 より ζ ( 2 s ) ≠ 0 \zeta(2s) \ne 0 ζ ( 2 s ) = 0 なので割り算も正当です)。
演習 7.3 標準
定理 4.5 と ζ ( 4 ) = π 4 / 90 \zeta(4) = \pi^{4}/90 ζ ( 4 ) = π 4 /90 を使って ζ ( − 3 ) \zeta(-3) ζ ( − 3 ) を求めてください。
解答 関数等式
ζ ( s ) = 2 s π s − 1 sin ( π s 2 ) Γ ( 1 − s ) ζ ( 1 − s ) \zeta(s) = 2^{s}\pi^{s-1}\sin\!\left(\frac{\pi s}{2}\right)\Gamma(1-s)\zeta(1-s) ζ ( s ) = 2 s π s − 1 sin ( 2 π s ) Γ ( 1 − s ) ζ ( 1 − s ) に s = − 3 s = -3 s = − 3 を代入します。各因子を順に計算します。
2 − 3 = 1 8 2^{-3} = \dfrac18 2 − 3 = 8 1 、π − 3 − 1 = π − 4 \pi^{-3-1} = \pi^{-4} π − 3 − 1 = π − 4 。
sin ( − 3 π 2 ) = + 1 \sin\left(\dfrac{-3\pi}{2}\right) = +1 sin ( 2 − 3 π ) = + 1 (− 3 π / 2 -3\pi/2 − 3 π /2 は π / 2 \pi/2 π /2 と 2 π 2\pi 2 π の差なので値は 1 1 1 です)。
Γ ( 1 − ( − 3 ) ) = Γ ( 4 ) = 3 ! = 6 \Gamma(1-(-3)) = \Gamma(4) = 3! = 6 Γ ( 1 − ( − 3 )) = Γ ( 4 ) = 3 ! = 6 。
ζ ( 1 − ( − 3 ) ) = ζ ( 4 ) = π 4 90 \zeta(1-(-3)) = \zeta(4) = \dfrac{\pi^{4}}{90} ζ ( 1 − ( − 3 )) = ζ ( 4 ) = 90 π 4 。
したがって
ζ ( − 3 ) = 1 8 ⋅ 1 π 4 ⋅ 1 ⋅ 6 ⋅ π 4 90 = 6 8 ⋅ 90 = 1 120 . \zeta(-3) = \frac{1}{8}\cdot\frac{1}{\pi^{4}}\cdot 1 \cdot 6 \cdot \frac{\pi^{4}}{90} = \frac{6}{8 \cdot 90} = \frac{1}{120}. ζ ( − 3 ) = 8 1 ⋅ π 4 1 ⋅ 1 ⋅ 6 ⋅ 90 π 4 = 8 ⋅ 90 6 = 120 1 . π 4 \pi^{4} π 4 が約分されて有理数になる点に注意してください。一般に負の奇数点 ζ ( 1 − 2 m ) \zeta(1-2m) ζ ( 1 − 2 m ) はベルヌーイ数を用いて − B 2 m / ( 2 m ) -B_{2m}/(2m) − B 2 m / ( 2 m ) と書ける有理数になります。実際 m = 2 m=2 m = 2 では B 4 = − 1 / 30 B_4 = -1/30 B 4 = − 1/30 より − B 4 / 4 = 1 / 120 -B_4/4 = 1/120 − B 4 /4 = 1/120 で、上の計算と一致します。
演習 7.4 難
Re s > 1 \operatorname{Re} s > 1 Re s > 1 において
− ζ ′ ( s ) ζ ( s ) = ∑ n = 1 ∞ Λ ( n ) n s -\frac{\zeta'(s)}{\zeta(s)} = \sum_{n=1}^{\infty}\frac{\Lambda(n)}{n^{s}} − ζ ( s ) ζ ′ ( s ) = n = 1 ∑ ∞ n s Λ ( n ) が成り立つことを示してください(Λ \Lambda Λ は 定義 6.1 のフォン・マンゴルト関数)。これが 注意 6.2 の明示公式の出発点です。
解答 σ 0 > 1 \sigma_0 > 1 σ 0 > 1 を固定し、半平面 Re s ≥ σ 0 \operatorname{Re} s \ge \sigma_0 Re s ≥ σ 0 で考えます。系 3.2 より ζ ( s ) ≠ 0 \zeta(s)\ne0 ζ ( s ) = 0 なので ζ ′ / ζ \zeta'/\zeta ζ ′ / ζ はここで正則です。
命題 3.4 の証明と同じ展開により、Re s > 1 \operatorname{Re} s > 1 Re s > 1 で
log ζ ( s ) = ∑ p ∑ k = 1 ∞ p − k s k \log \zeta(s) = \sum_{p}\sum_{k=1}^{\infty}\frac{p^{-ks}}{k} log ζ ( s ) = p ∑ k = 1 ∑ ∞ k p − k s が成り立ちます(対数は ζ \zeta ζ の連続な分枝を取ります)。右辺の二重級数は Re s ≥ σ 0 \operatorname{Re} s \ge \sigma_0 Re s ≥ σ 0 で
∑ p ∑ k p − k σ 0 k ≤ ∑ p ∑ k p − k σ 0 ≤ 2 ∑ p p − σ 0 ≤ 2 ζ ( σ 0 ) < ∞ \sum_p \sum_k \frac{p^{-k\sigma_0}}{k} \le \sum_p\sum_k p^{-k\sigma_0} \le 2\sum_p p^{-\sigma_0} \le 2\zeta(\sigma_0) < \infty p ∑ k ∑ k p − k σ 0 ≤ p ∑ k ∑ p − k σ 0 ≤ 2 p ∑ p − σ 0 ≤ 2 ζ ( σ 0 ) < ∞ と一様に評価できるので、広義一様収束します。正則関数列が広義一様収束するとき項別微分できる(ワイヤシュトラスの定理)ので、d d s p − k s = − ( k log p ) p − k s \frac{d}{ds}p^{-ks} = -(k\log p)\,p^{-ks} d s d p − k s = − ( k log p ) p − k s を使って
ζ ′ ( s ) ζ ( s ) = d d s log ζ ( s ) = ∑ p ∑ k ≥ 1 − ( k log p ) p − k s k = − ∑ p ∑ k ≥ 1 ( log p ) p − k s . \frac{\zeta'(s)}{\zeta(s)} = \frac{d}{ds}\log\zeta(s) = \sum_p \sum_{k\ge1} \frac{-(k\log p)p^{-ks}}{k} = -\sum_p\sum_{k \ge 1}(\log p)\, p^{-ks}. ζ ( s ) ζ ′ ( s ) = d s d log ζ ( s ) = p ∑ k ≥ 1 ∑ k − ( k log p ) p − k s = − p ∑ k ≥ 1 ∑ ( log p ) p − k s . 符号を移すと
− ζ ′ ( s ) ζ ( s ) = ∑ p ∑ k ≥ 1 log p ( p k ) s . -\frac{\zeta'(s)}{\zeta(s)} = \sum_p \sum_{k\ge1}\frac{\log p}{\left(p^{k}\right)^{s}} . − ζ ( s ) ζ ′ ( s ) = p ∑ k ≥ 1 ∑ ( p k ) s log p . 右辺の和は「素数べき n = p k n = p^{k} n = p k 」を走り、そのとき係数は log p = Λ ( n ) \log p = \Lambda(n) log p = Λ ( n ) です。素数べきでない n n n には Λ ( n ) = 0 \Lambda(n)=0 Λ ( n ) = 0 を割り当てればよいので、右辺は ∑ n ≥ 1 Λ ( n ) n − s \sum_{n \ge 1}\Lambda(n)n^{-s} ∑ n ≥ 1 Λ ( n ) n − s に等しくなります。素因数分解の一意性より、各素数べき n n n を与える組 ( p , k ) (p,k) ( p , k ) は一意なので、重複はありません。
B. Riemann, “Über die Anzahl der Primzahlen unter einer gegebenen Grösse”, Monatsberichte der Berliner Akademie , 1859. — 関数等式・明示公式・リーマン予想が最初に現れた原論文。
H. M. Edwards, Riemann’s Zeta Function , Academic Press, 1974(Dover 版 2001)— 第 1 章がリーマンの原論文の逐条解説、第 3 章が明示公式。
E. C. Titchmarsh, The Theory of the Riemann Zeta-Function , 2nd ed. (revised by D. R. Heath-Brown), Oxford University Press, 1986. — 第 2 章に関数等式、第 14 章にリーマン予想と誤差項の同値性。
T. M. Apostol, Introduction to Analytic Number Theory , Springer, 1976. — 第 11〜13 章。オイラー積から素数定理までを学部レベルで丁寧に扱う。
H. von Koch, “Sur la distribution des nombres premiers”, Acta Mathematica 24 (1901), 159–182. — 定理 6.3 の原論文。
J. B. Conrey, “The Riemann Hypothesis”, Notices of the American Mathematical Society 50 (2003), 341–353. — 現代的な研究状況の概観。
明示公式の各項を波として読む方法を補足します。 注意 6.2 の和に現れる零点は共役対 ρ = 1 2 + i γ \rho = \tfrac12 + i\gamma ρ = 2 1 + iγ と ρ ˉ = 1 2 − i γ \bar\rho = \tfrac12 - i\gamma ρ ˉ = 2 1 − iγ で現れます(定理 4.5 の後で述べた実軸対称性によります)。リーマン予想を仮定してこの対の寄与をまとめると、x ρ = x 1 / 2 e i γ log x x^{\rho} = x^{1/2}e^{i\gamma\log x} x ρ = x 1/2 e iγ l o g x と x ρ ‾ = x 1 / 2 e − i γ log x \overline{x^{\rho}} = x^{1/2}e^{-i\gamma\log x} x ρ = x 1/2 e − iγ l o g x の和になるので
x ρ ρ + x ρ ˉ ρ ˉ = 2 Re x ρ ρ = 2 x ∣ ρ ∣ cos ( γ log x − arg ρ ) \frac{x^{\rho}}{\rho} + \frac{x^{\bar\rho}}{\bar\rho} = 2\operatorname{Re}\frac{x^{\rho}}{\rho} = \frac{2\sqrt{x}}{|\rho|}\cos\!\left(\gamma\log x - \arg\rho\right) ρ x ρ + ρ ˉ x ρ ˉ = 2 Re ρ x ρ = ∣ ρ ∣ 2 x cos ( γ log x − arg ρ )
という実数の式になります。これは変数 u = log x u = \log x u = log x について周期 2 π / γ 2\pi/\gamma 2 π / γ の余弦波で、振幅は 2 x / ∣ ρ ∣ 2\sqrt{x}/|\rho| 2 x /∣ ρ ∣ です。
つまり明示公式は「ψ ( x ) = x − ( 無数の余弦波の重ね合わせ ) \psi(x) = x - (\text{無数の余弦波の重ね合わせ}) ψ ( x ) = x − ( 無数の余弦波の重ね合わせ ) 」と読めます。素数の分布は、主要項 x x x に零点由来の波が干渉して階段状の形が作られる、という描像です。振幅は 2 x / ∣ ρ ∣ 2\sqrt{x}/|\rho| 2 x /∣ ρ ∣ なので、虚部 γ \gamma γ の小さい零点ほど大きく効きます。γ 1 = 14.13 \gamma_1 = 14.13 γ 1 = 14.13 の零点が最も強い波を出し、その周期は log x \log x log x について 2 π / 14.13 ≈ 0.445 2\pi/14.13 \approx 0.445 2 π /14.13 ≈ 0.445 、すなわち x x x が e 0.445 ≈ 1.56 e^{0.445} \approx 1.56 e 0.445 ≈ 1.56 倍になるごとに一周します。
もし臨界線から外れた零点 ρ = β + i γ \rho = \beta + i\gamma ρ = β + iγ (β > 1 / 2 \beta > 1/2 β > 1/2 )が存在すれば、その波の振幅は x β x^{\beta} x β に比例して他より速く育ち、いずれすべての波を圧倒します。リーマン予想とは「どの波も突出しない」、言い換えれば素数の分布が主要項のまわりで可能な限り整然と揺れている、という主張なのです。
定理 6.3 の証明に使う補題をまとめます。s = σ + i t s=\sigma+it s = σ + i t と書き、θ ( x ) = ∑ p ≤ x log p \theta(x)=\sum_{p\le x}\log p θ ( x ) = ∑ p ≤ x log p 、ψ ( x ) = ∑ n ≤ x Λ ( n ) \psi(x)=\sum_{n\le x}\Lambda(n) ψ ( x ) = ∑ n ≤ x Λ ( n ) (定義 6.1 )とします。零点は位数の分だけ重複して数えます。
補題群は二つの向きに分かれます。リーマン予想から誤差の評価を出す向き(補題 7.11 )には、ゼータ関数の大きさ、零点の局所的な個数、打ち切ったペロンの公式の補題を順に使います。逆向き、つまり誤差の評価から零点の位置を出す向きには 補題 7.12 を使います。最初の補題は両方に共通する準備です。
補題 7.5 (チェビシェフ関数の評価と対数微分 ) ✓
x ≥ 1 x\ge1 x ≥ 1 について θ ( x ) ≤ x log 4 \theta(x)\le x\log4 θ ( x ) ≤ x log 4 、0 ≤ ψ ( x ) − θ ( x ) ≤ 2 x log x 0\le\psi(x)-\theta(x)\le2\sqrt x\log x 0 ≤ ψ ( x ) − θ ( x ) ≤ 2 x log x 、ψ ( x ) ≤ 3 x \psi(x)\le3x ψ ( x ) ≤ 3 x 。
σ > 1 \sigma\gt1 σ > 1 で − ζ ′ ( s ) ζ ( s ) = ∑ n ≥ 1 Λ ( n ) n s = s ∫ 1 ∞ ψ ( x ) x − s − 1 d x -\dfrac{\zeta'(s)}{\zeta(s)}=\displaystyle\sum_{n\ge1}\frac{\Lambda(n)}{n^s}=s\int_1^\infty\psi(x)x^{-s-1}dx − ζ ( s ) ζ ′ ( s ) = n ≥ 1 ∑ n s Λ ( n ) = s ∫ 1 ∞ ψ ( x ) x − s − 1 d x で、0 ≤ ∑ n Λ ( n ) n − σ ≤ 3 σ σ − 1 0\le\sum_n\Lambda(n)n^{-\sigma}\le\dfrac{3\sigma}{\sigma-1} 0 ≤ ∑ n Λ ( n ) n − σ ≤ σ − 1 3 σ 。
証明(補題 7.5) 1. θ ( x ) ≤ x log 4 \theta(x)\le x\log4 θ ( x ) ≤ x log 4 は 命題 5.3[素数の魅力と素数定理] の証明(同記事 Appendix の第 1 段)で示されています。素数冪 p k ≤ x p^k\le x p k ≤ x を k k k ごとにまとめると ψ ( x ) = ∑ k ≥ 1 θ ( x 1 / k ) \psi(x)=\sum_{k\ge1}\theta(x^{1/k}) ψ ( x ) = ∑ k ≥ 1 θ ( x 1/ k ) で、k > log x / log 2 k\gt\log x/\log2 k > log x / log 2 では x 1 / k < 2 x^{1/k}\lt2 x 1/ k < 2 なので項は 0 0 0 です。k ≥ 2 k\ge2 k ≥ 2 の項は θ ( x ) ≤ x log 4 \theta(\sqrt x)\le\sqrt x\log4 θ ( x ) ≤ x log 4 以下なので、ψ ( x ) − θ ( x ) ≤ log x log 2 x log 4 = 2 x log x \psi(x)-\theta(x)\le\frac{\log x}{\log2}\sqrt x\log4=2\sqrt x\log x ψ ( x ) − θ ( x ) ≤ l o g 2 l o g x x log 4 = 2 x log x 。log x / x ≤ 2 / e \log x/\sqrt x\le2/e log x / x ≤ 2/ e から ψ ( x ) ≤ ( log 4 + 4 / e ) x ≤ 3 x \psi(x)\le(\log4+4/e)x\le3x ψ ( x ) ≤ ( log 4 + 4/ e ) x ≤ 3 x です。
2. 定理 4.2[素数の魅力と素数定理] で n = ∏ p e p n=\prod p^{e_p} n = ∏ p e p と分解すると、n n n の約数のうち素数冪であるものは p k p^k p k (1 ≤ k ≤ e p 1\le k\le e_p 1 ≤ k ≤ e p )なので ∑ d ∣ n Λ ( d ) = log n \sum_{d\mid n}\Lambda(d)=\log n ∑ d ∣ n Λ ( d ) = log n です。σ > 1 \sigma\gt1 σ > 1 では ∑ n − s \sum n^{-s} ∑ n − s と ∑ Λ ( n ) n − s \sum\Lambda(n)n^{-s} ∑ Λ ( n ) n − s (0 ≤ Λ ( n ) ≤ log n 0\le\Lambda(n)\le\log n 0 ≤ Λ ( n ) ≤ log n )が絶対収束するので、積を ( d , m ) (d,m) ( d , m ) を添字とする絶対収束する二重級数とみて n = d m n=dm n = d m ごとにまとめて(定理 7.11[級数と収束判定] の 2) ζ ( s ) ∑ d Λ ( d ) d − s = ∑ n ( log n ) n − s = − ζ ′ ( s ) \zeta(s)\sum_d\Lambda(d)d^{-s}=\sum_n(\log n)n^{-s}=-\zeta'(s) ζ ( s ) ∑ d Λ ( d ) d − s = ∑ n ( log n ) n − s = − ζ ′ ( s ) です(最後は部分和が σ ≥ σ 0 > 1 \sigma\ge\sigma_0\gt1 σ ≥ σ 0 > 1 で一様収束するので項別微分できること、定理 5.2[偏角の原理と Rouché の定理] )。ζ ( s ) ≠ 0 \zeta(s)\ne0 ζ ( s ) = 0 (系 3.2 )で割れば第 1 の等式です。次に n − s = s ∫ n ∞ x − s − 1 d x n^{-s}=s\int_n^\infty x^{-s-1}dx n − s = s ∫ n ∞ x − s − 1 d x を代入し、∑ n Λ ( n ) ∣ s ∣ ∫ n ∞ x − σ − 1 d x = ∣ s ∣ σ ∑ Λ ( n ) n − σ < ∞ \sum_n\Lambda(n)|s|\int_n^\infty x^{-\sigma-1}dx=\frac{|s|}\sigma\sum\Lambda(n)n^{-\sigma}\lt\infty ∑ n Λ ( n ) ∣ s ∣ ∫ n ∞ x − σ − 1 d x = σ ∣ s ∣ ∑ Λ ( n ) n − σ < ∞ なので和と積分を交換できて、∑ n Λ ( n ) n − s = s ∫ 1 ∞ ( ∑ n ≤ x Λ ( n ) ) x − s − 1 d x \sum_n\Lambda(n)n^{-s}=s\int_1^\infty\bigl(\sum_{n\le x}\Lambda(n)\bigr)x^{-s-1}dx ∑ n Λ ( n ) n − s = s ∫ 1 ∞ ( ∑ n ≤ x Λ ( n ) ) x − s − 1 d x です。s = σ s=\sigma s = σ とし 1 の ψ ( x ) ≤ 3 x \psi(x)\le3x ψ ( x ) ≤ 3 x を使えば最後の評価です。
∎
補題 7.6 (臨界帯の近くでのゼータ関数の大きさ ) ✓
P 2 ( x ) = 1 2 ( { x } 2 − { x } ) P_2(x)=\frac12(\{x\}^2-\{x\}) P 2 ( x ) = 2 1 ({ x } 2 − { x }) ({ x } \{x\} { x } は小数部分)とおく。σ > − 1 \sigma\gt-1 σ > − 1 、s ≠ 1 s\ne1 s = 1 で
ζ ( s ) = s s − 1 − 1 2 − s ( s + 1 ) ∫ 1 ∞ P 2 ( x ) x − s − 2 d x \zeta(s)=\frac{s}{s-1}-\frac12-s(s+1)\int_1^\infty P_2(x)\,x^{-s-2}\,dx ζ ( s ) = s − 1 s − 2 1 − s ( s + 1 ) ∫ 1 ∞ P 2 ( x ) x − s − 2 d x が成り立つ。とくに − 1 2 ≤ σ ≤ 3 -\frac12\le\sigma\le3 − 2 1 ≤ σ ≤ 3 、∣ t ∣ ≥ 2 |t|\ge2 ∣ t ∣ ≥ 2 で ∣ ζ ( s ) ∣ ≤ 3 ∣ t ∣ 2 |\zeta(s)|\le3|t|^2 ∣ ζ ( s ) ∣ ≤ 3∣ t ∣ 2 、σ ≥ 2 \sigma\ge2 σ ≥ 2 で ∣ ζ ( s ) ∣ ≤ 2 |\zeta(s)|\le2 ∣ ζ ( s ) ∣ ≤ 2 、また ∣ ζ ( 2 + i t ) ∣ ≥ c 0 : = 1 − ∑ n ≥ 2 n − 2 > 0.3 |\zeta(2+it)|\ge c_0:=1-\sum_{n\ge2}n^{-2}\gt0.3 ∣ ζ ( 2 + i t ) ∣ ≥ c 0 := 1 − ∑ n ≥ 2 n − 2 > 0.3 。
証明(補題 7.6) σ > 1 \sigma\gt1 σ > 1 とします。n ( n − s − ( n + 1 ) − s ) = n 1 − s − ( n + 1 ) 1 − s + ( n + 1 ) − s n(n^{-s}-(n+1)^{-s})=n^{1-s}-(n+1)^{1-s}+(n+1)^{-s} n ( n − s − ( n + 1 ) − s ) = n 1 − s − ( n + 1 ) 1 − s + ( n + 1 ) − s を n = 1 , … , N n=1,\ldots,N n = 1 , … , N で足して N → ∞ N\to\infty N → ∞ とすると、s ∫ 1 ∞ ⌊ x ⌋ x − s − 1 d x = ∑ n n ⋅ s ∫ n n + 1 x − s − 1 d x = ζ ( s ) s\int_1^\infty\lfloor x\rfloor x^{-s-1}dx=\sum_nn\cdot s\int_n^{n+1}x^{-s-1}dx=\zeta(s) s ∫ 1 ∞ ⌊ x ⌋ x − s − 1 d x = ∑ n n ⋅ s ∫ n n + 1 x − s − 1 d x = ζ ( s ) です。P 1 ( x ) = { x } − 1 2 P_1(x)=\{x\}-\frac12 P 1 ( x ) = { x } − 2 1 とおくと ⌊ x ⌋ = x − 1 2 − P 1 ( x ) \lfloor x\rfloor=x-\frac12-P_1(x) ⌊ x ⌋ = x − 2 1 − P 1 ( x ) で、s ∫ 1 ∞ x − s d x = s s − 1 s\int_1^\infty x^{-s}dx=\frac s{s-1} s ∫ 1 ∞ x − s d x = s − 1 s 、s ∫ 1 ∞ 1 2 x − s − 1 d x = 1 2 s\int_1^\infty\frac12x^{-s-1}dx=\frac12 s ∫ 1 ∞ 2 1 x − s − 1 d x = 2 1 なので ζ ( s ) = s s − 1 − 1 2 − s ∫ 1 ∞ P 1 ( x ) x − s − 1 d x \zeta(s)=\frac s{s-1}-\frac12-s\int_1^\infty P_1(x)x^{-s-1}dx ζ ( s ) = s − 1 s − 2 1 − s ∫ 1 ∞ P 1 ( x ) x − s − 1 d x 。P 2 P_2 P 2 は連続で P 2 ′ = P 1 P_2'=P_1 P 2 ′ = P 1 (整数点以外)、P 2 P_2 P 2 は整数点で 0 0 0 、− 1 8 ≤ P 2 ≤ 0 -\frac18\le P_2\le0 − 8 1 ≤ P 2 ≤ 0 なので、部分積分で ∫ 1 ∞ P 1 x − s − 1 d x = ( s + 1 ) ∫ 1 ∞ P 2 x − s − 2 d x \int_1^\infty P_1x^{-s-1}dx=(s+1)\int_1^\infty P_2x^{-s-2}dx ∫ 1 ∞ P 1 x − s − 1 d x = ( s + 1 ) ∫ 1 ∞ P 2 x − s − 2 d x を得て、σ > 1 \sigma\gt1 σ > 1 で等式が成り立ちます。
右辺の積分は、σ ≥ − 1 + δ \sigma\ge-1+\delta σ ≥ − 1 + δ で被積分関数の絶対値が 1 8 x − 1 − δ \frac18x^{-1-\delta} 8 1 x − 1 − δ 以下なので、∫ 1 N \int_1^N ∫ 1 N (Fubini と Morera の定理 定理 3.6[正則関数の強力な性質] で整関数)の局所一様な極限として σ > − 1 \sigma\gt-1 σ > − 1 で正則です(定理 5.2[偏角の原理と Rouché の定理] )。ζ \zeta ζ は C \mathbb{C} C 上の有理型関数(定理 4.5 )なので、一致の定理(定理 7.3[正則関数の強力な性質] )により等式は σ > − 1 \sigma\gt-1 σ > − 1 、s ≠ 1 s\ne1 s = 1 で成り立ちます。
評価: − 1 2 ≤ σ ≤ 3 -\frac12\le\sigma\le3 − 2 1 ≤ σ ≤ 3 、∣ t ∣ ≥ 2 |t|\ge2 ∣ t ∣ ≥ 2 では ∣ s ∣ ≤ 3 + ∣ t ∣ ≤ 5 2 ∣ t ∣ |s|\le3+|t|\le\frac52|t| ∣ s ∣ ≤ 3 + ∣ t ∣ ≤ 2 5 ∣ t ∣ 、∣ s + 1 ∣ ≤ 4 + ∣ t ∣ ≤ 3 ∣ t ∣ |s+1|\le4+|t|\le3|t| ∣ s + 1∣ ≤ 4 + ∣ t ∣ ≤ 3∣ t ∣ 、∣ s − 1 ∣ ≥ ∣ t ∣ |s-1|\ge|t| ∣ s − 1∣ ≥ ∣ t ∣ 、積分は 1 8 ( σ + 1 ) ≤ 1 4 \frac1{8(\sigma+1)}\le\frac14 8 ( σ + 1 ) 1 ≤ 4 1 以下なので ∣ ζ ( s ) ∣ ≤ 5 2 + 1 2 + 15 8 ∣ t ∣ 2 ≤ 3 ∣ t ∣ 2 |\zeta(s)|\le\frac52+\frac12+\frac{15}8|t|^2\le3|t|^2 ∣ ζ ( s ) ∣ ≤ 2 5 + 2 1 + 8 15 ∣ t ∣ 2 ≤ 3∣ t ∣ 2 。σ ≥ 2 \sigma\ge2 σ ≥ 2 では ∣ ζ ( s ) ∣ ≤ ∑ n − 2 < 2 |\zeta(s)|\le\sum n^{-2}\lt2 ∣ ζ ( s ) ∣ ≤ ∑ n − 2 < 2 、∣ ζ ( 2 + i t ) ∣ ≥ 1 − ∑ n ≥ 2 n − 2 |\zeta(2+it)|\ge1-\sum_{n\ge2}n^{-2} ∣ ζ ( 2 + i t ) ∣ ≥ 1 − ∑ n ≥ 2 n − 2 で、∑ n ≥ 2 n − 2 ≤ 1 4 + 1 9 + ∑ n ≥ 4 1 n ( n − 1 ) = 1 4 + 1 9 + 1 3 < 0.7 \sum_{n\ge2}n^{-2}\le\frac14+\frac19+\sum_{n\ge4}\frac1{n(n-1)}=\frac14+\frac19+\frac13\lt0.7 ∑ n ≥ 2 n − 2 ≤ 4 1 + 9 1 + ∑ n ≥ 4 n ( n − 1 ) 1 = 4 1 + 9 1 + 3 1 < 0.7 です。
∎
補題 7.7 (Borel–Carathéodory 型の評価 ) ✓
h h h を閉円板 ∣ z − z 0 ∣ ≤ R |z-z_0|\le R ∣ z − z 0 ∣ ≤ R を含む開集合で正則とし、h ( z 0 ) = 0 h(z_0)=0 h ( z 0 ) = 0 、円周 ∣ z − z 0 ∣ = R |z-z_0|=R ∣ z − z 0 ∣ = R 上で R e h ≤ A \mathrm{Re}\,h\le A Re h ≤ A とする。このとき ∣ z − z 0 ∣ < R |z-z_0|\lt R ∣ z − z 0 ∣ < R で ∣ h ′ ( z ) ∣ ≤ 2 A R ( R − ∣ z − z 0 ∣ ) 2 |h'(z)|\le\dfrac{2AR}{(R-|z-z_0|)^2} ∣ h ′ ( z ) ∣ ≤ ( R − ∣ z − z 0 ∣ ) 2 2 A R 。
証明(補題 7.7) z 0 = 0 z_0=0 z 0 = 0 とし、∣ z ∣ < R |z|\lt R ∣ z ∣ < R を固定します。導関数版の積分公式(定理 3.2[正則関数の強力な性質] )から h ′ ( z ) = 1 2 π i ∮ ∣ w ∣ = R h ( w ) ( w − z ) 2 d w h'(z)=\frac1{2\pi i}\oint_{|w|=R}\frac{h(w)}{(w-z)^2}dw h ′ ( z ) = 2 π i 1 ∮ ∣ w ∣ = R ( w − z ) 2 h ( w ) d w 。h = ∑ n a n w n h=\sum_na_nw^n h = ∑ n a n w n は ∣ w ∣ = R |w|=R ∣ w ∣ = R で一様収束し(定理 5.2[正則関数の強力な性質] )、円周上では h ( w ) ‾ = ∑ n a ˉ n R 2 n w − n \overline{h(w)}=\sum_n\bar a_nR^{2n}w^{-n} h ( w ) = ∑ n a ˉ n R 2 n w − n 、1 ( w − z ) 2 = ∑ k ≥ 0 ( k + 1 ) z k w − k − 2 \frac1{(w-z)^2}=\sum_{k\ge0}(k+1)z^kw^{-k-2} ( w − z ) 2 1 = ∑ k ≥ 0 ( k + 1 ) z k w − k − 2 なので、∮ h ( w ) ‾ ( w − z ) 2 d w \oint\frac{\overline{h(w)}}{(w-z)^2}dw ∮ ( w − z ) 2 h ( w ) d w は w − 1 w^{-1} w − 1 の項をもたない級数の項別積分で 0 0 0 です。∮ d w ( w − z ) 2 = 0 \oint\frac{dw}{(w-z)^2}=0 ∮ ( w − z ) 2 d w = 0 も同様です。よって u = R e h u=\mathrm{Re}\,h u = Re h について
h ′ ( z ) = 1 2 π i ∮ ∣ w ∣ = R 2 ( u ( w ) − A ) ( w − z ) 2 d w , ∣ h ′ ( z ) ∣ ≤ 2 R ( R − ∣ z ∣ ) 2 ⋅ 1 2 π ∫ 0 2 π ( A − u ( R e i ϕ ) ) d ϕ . h'(z)=\frac1{2\pi i}\oint_{|w|=R}\frac{2(u(w)-A)}{(w-z)^2}dw,\qquad|h'(z)|\le\frac{2R}{(R-|z|)^2}\cdot\frac1{2\pi}\int_0^{2\pi}\bigl(A-u(Re^{i\phi})\bigr)d\phi . h ′ ( z ) = 2 π i 1 ∮ ∣ w ∣ = R ( w − z ) 2 2 ( u ( w ) − A ) d w , ∣ h ′ ( z ) ∣ ≤ ( R − ∣ z ∣ ) 2 2 R ⋅ 2 π 1 ∫ 0 2 π ( A − u ( R e i ϕ ) ) d ϕ . Cauchy の積分公式(定理 6.1[コーシーの積分定理と積分公式] )の実部から 1 2 π ∫ 0 2 π u ( R e i ϕ ) d ϕ = u ( 0 ) = 0 \frac1{2\pi}\int_0^{2\pi}u(Re^{i\phi})d\phi=u(0)=0 2 π 1 ∫ 0 2 π u ( R e i ϕ ) d ϕ = u ( 0 ) = 0 なので、右辺は 2 A R ( R − ∣ z ∣ ) 2 \frac{2AR}{(R-|z|)^2} ( R − ∣ z ∣ ) 2 2 A R です。
∎
補題 7.8 (零点の近くでの対数微分(Landau の補題) ) ✓
f f f を閉円板 ∣ z − z 0 ∣ ≤ 5 2 |z-z_0|\le\frac52 ∣ z − z 0 ∣ ≤ 2 5 を含む開集合で正則とし、f ( z 0 ) ≠ 0 f(z_0)\ne0 f ( z 0 ) = 0 、M ≥ 1 M\ge1 M ≥ 1 について円板上で ∣ f ( z ) ∣ ≤ e M ∣ f ( z 0 ) ∣ |f(z)|\le e^M|f(z_0)| ∣ f ( z ) ∣ ≤ e M ∣ f ( z 0 ) ∣ とする。Z Z Z を ∣ ρ − z 0 ∣ ≤ 2 |\rho-z_0|\le2 ∣ ρ − z 0 ∣ ≤ 2 にある f f f の零点の多重集合とすると、(a) # Z ≤ M / log 5 4 \#Z\le M/\log\frac54 # Z ≤ M / log 4 5 、(b) ∣ z − z 0 ∣ ≤ 7 4 |z-z_0|\le\frac74 ∣ z − z 0 ∣ ≤ 4 7 で f ( z ) ≠ 0 f(z)\ne0 f ( z ) = 0 なら ∣ f ′ ( z ) f ( z ) − ∑ ρ ∈ Z 1 z − ρ ∣ ≤ 470 M \Bigl|\dfrac{f'(z)}{f(z)}-\displaystyle\sum_{\rho\in Z}\frac1{z-\rho}\Bigr|\le470M f ( z ) f ′ ( z ) − ρ ∈ Z ∑ z − ρ 1 ≤ 470 M 。
証明(補題 7.8) (a) Jensen の公式(定理 7.2[Jensen の公式と劣調和関数] )を半径 5 2 \frac52 2 5 で使うと、右辺の和の各項は非負で、∣ ρ − z 0 ∣ ≤ 2 |\rho-z_0|\le2 ∣ ρ − z 0 ∣ ≤ 2 の零点の項は log 5 4 \log\frac54 log 4 5 以上、左辺は log ∣ f ( z 0 ) ∣ + M \log|f(z_0)|+M log ∣ f ( z 0 ) ∣ + M 以下なので、# Z log 5 4 ≤ M \#Z\log\frac54\le M # Z log 4 5 ≤ M です。
(b) n = # Z n=\#Z n = # Z とし、g ( z ) = f ( z ) ∏ ρ ∈ Z ( z − ρ ) − 1 g(z)=f(z)\prod_{\rho\in Z}(z-\rho)^{-1} g ( z ) = f ( z ) ∏ ρ ∈ Z ( z − ρ ) − 1 とおくと、g g g は除去可能特異点を除いて同じ開集合で正則で、∣ z − z 0 ∣ ≤ 2 |z-z_0|\le2 ∣ z − z 0 ∣ ≤ 2 に零点をもちません。∣ z − z 0 ∣ = 5 2 |z-z_0|=\frac52 ∣ z − z 0 ∣ = 2 5 では ∣ z − ρ ∣ ≥ 1 2 |z-\rho|\ge\frac12 ∣ z − ρ ∣ ≥ 2 1 、∣ z 0 − ρ ∣ ≤ 2 |z_0-\rho|\le2 ∣ z 0 − ρ ∣ ≤ 2 なので ∣ g ( z ) g ( z 0 ) ∣ = ∣ f ( z ) f ( z 0 ) ∣ ∏ ρ ∣ z 0 − ρ ∣ ∣ z − ρ ∣ ≤ e M 4 n \bigl|\frac{g(z)}{g(z_0)}\bigr|=\bigl|\frac{f(z)}{f(z_0)}\bigr|\prod_\rho\frac{|z_0-\rho|}{|z-\rho|}\le e^M4^n g ( z 0 ) g ( z ) = f ( z 0 ) f ( z ) ∏ ρ ∣ z − ρ ∣ ∣ z 0 − ρ ∣ ≤ e M 4 n で、最大値原理(定理 7.7[正則関数の強力な性質] )から閉円板全体でも同じ評価が成り立ちます。g g g は ∣ z − z 0 ∣ ≤ 2 |z-z_0|\le2 ∣ z − z 0 ∣ ≤ 2 を含む開円板で零点をもたないので、g ′ / g g'/g g ′ / g の原始関数 h h h で h ( z 0 ) = 0 h(z_0)=0 h ( z 0 ) = 0 となるものが取れ(凸領域では Goursat の定理 定理 4.5[コーシーの積分定理と積分公式] から h ( z ) = ∫ [ z 0 , z ] g ′ / g h(z)=\int_{[z_0,z]}g'/g h ( z ) = ∫ [ z 0 , z ] g ′ / g が原始関数)、( g e − h ) ′ = 0 (ge^{-h})'=0 ( g e − h ) ′ = 0 から g = g ( z 0 ) e h g=g(z_0)e^h g = g ( z 0 ) e h 、R e h = log ∣ g / g ( z 0 ) ∣ ≤ M + n log 4 \mathrm{Re}\,h=\log|g/g(z_0)|\le M+n\log4 Re h = log ∣ g / g ( z 0 ) ∣ ≤ M + n log 4 です。補題 7.7 を半径 2 2 2 で使うと、∣ z − z 0 ∣ ≤ 7 4 |z-z_0|\le\frac74 ∣ z − z 0 ∣ ≤ 4 7 で
∣ f ′ f ( z ) − ∑ ρ ∈ Z 1 z − ρ ∣ = ∣ h ′ ( z ) ∣ ≤ 2 ⋅ 2 ( M + n log 4 ) ( 1 / 4 ) 2 ≤ 64 ( 1 + log 4 log ( 5 / 4 ) ) M ≤ 470 M . \Bigl|\frac{f'}{f}(z)-\sum_{\rho\in Z}\frac1{z-\rho}\Bigr|=|h'(z)|\le\frac{2\cdot2(M+n\log4)}{(1/4)^2}\le64\Bigl(1+\frac{\log4}{\log(5/4)}\Bigr)M\le470M . f f ′ ( z ) − ρ ∈ Z ∑ z − ρ 1 = ∣ h ′ ( z ) ∣ ≤ ( 1/4 ) 2 2 ⋅ 2 ( M + n log 4 ) ≤ 64 ( 1 + log ( 5/4 ) log 4 ) M ≤ 470 M . ∎
補題 7.9 (ゼータ関数の零点の局所的な個数 ) ✓
∣ t ∣ ≥ 5 |t|\ge5 ∣ t ∣ ≥ 5 とし、Z ( t ) Z(t) Z ( t ) を ∣ ρ − ( 2 + i t ) ∣ ≤ 2 |\rho-(2+it)|\le2 ∣ ρ − ( 2 + i t ) ∣ ≤ 2 にある ζ \zeta ζ の零点の多重集合とする。# Z ( t ) ≤ 18 log ∣ t ∣ \#Z(t)\le18\log|t| # Z ( t ) ≤ 18 log ∣ t ∣ で、∣ s − ( 2 + i t ) ∣ ≤ 7 4 |s-(2+it)|\le\frac74 ∣ s − ( 2 + i t ) ∣ ≤ 4 7 かつ ζ ( s ) ≠ 0 \zeta(s)\ne0 ζ ( s ) = 0 なら
∣ ζ ′ ( s ) ζ ( s ) − ∑ ρ ∈ Z ( t ) 1 s − ρ ∣ ≤ 1900 log ∣ t ∣ . \Bigl|\frac{\zeta'(s)}{\zeta(s)}-\sum_{\rho\in Z(t)}\frac1{s-\rho}\Bigr|\le1900\log|t| . ζ ( s ) ζ ′ ( s ) − ρ ∈ Z ( t ) ∑ s − ρ 1 ≤ 1900 log ∣ t ∣. とくに、零点 ρ = β + i γ \rho=\beta+i\gamma ρ = β + iγ で β ≥ 1 2 \beta\ge\frac12 β ≥ 2 1 、∣ γ − t ∣ ≤ 1 |\gamma-t|\le1 ∣ γ − t ∣ ≤ 1 を満たすものはすべて Z ( t ) Z(t) Z ( t ) に属する。
証明(補題 7.9) 閉円板 ∣ s − ( 2 + i t ) ∣ ≤ 5 2 |s-(2+it)|\le\frac52 ∣ s − ( 2 + i t ) ∣ ≤ 2 5 は s = 1 s=1 s = 1 を含まず、その上で − 1 2 ≤ σ ≤ 9 2 -\frac12\le\sigma\le\frac92 − 2 1 ≤ σ ≤ 2 9 、∣ I m s ∣ ∈ [ 2 , 3 2 ∣ t ∣ ] |\mathrm{Im}\,s|\in[2,\frac32|t|] ∣ Im s ∣ ∈ [ 2 , 2 3 ∣ t ∣ ] です。補題 7.6 から円板上で ∣ ζ ( s ) ∣ ≤ 3 ⋅ 9 4 t 2 ≤ 7 t 2 |\zeta(s)|\le3\cdot\frac94t^2\le7t^2 ∣ ζ ( s ) ∣ ≤ 3 ⋅ 4 9 t 2 ≤ 7 t 2 、∣ ζ ( 2 + i t ) ∣ ≥ c 0 > 0.3 |\zeta(2+it)|\ge c_0\gt0.3 ∣ ζ ( 2 + i t ) ∣ ≥ c 0 > 0.3 なので、M = log ( 24 t 2 ) ≤ 4 log ∣ t ∣ M=\log(24t^2)\le4\log|t| M = log ( 24 t 2 ) ≤ 4 log ∣ t ∣ (∣ t ∣ ≥ 5 |t|\ge5 ∣ t ∣ ≥ 5 なら 24 ≤ t 2 24\le t^2 24 ≤ t 2 )として 補題 7.8 が使えます。4 / log 5 4 ≤ 18 4/\log\frac54\le18 4/ log 4 5 ≤ 18 、470 ⋅ 4 ≤ 1900 470\cdot4\le1900 470 ⋅ 4 ≤ 1900 です。最後の主張は β ≤ 1 \beta\le1 β ≤ 1 (系 3.2 )から ∣ ρ − ( 2 + i t ) ∣ 2 ≤ ( 2 − 1 2 ) 2 + 1 < 4 |\rho-(2+it)|^2\le(2-\frac12)^2+1\lt4 ∣ ρ − ( 2 + i t ) ∣ 2 ≤ ( 2 − 2 1 ) 2 + 1 < 4 によります。
∎
補題 7.10 (打ち切った Perron の公式 ) ✓
y > 0 y\gt0 y > 0 、y ≠ 1 y\ne1 y = 1 、c > 0 c\gt0 c > 0 、T > 0 T\gt0 T > 0 について
∣ 1 2 π i ∫ c − i T c + i T y s s d s − 1 y > 1 ∣ ≤ y c π T ∣ log y ∣ . \Bigl|\frac1{2\pi i}\int_{c-iT}^{c+iT}\frac{y^s}{s}\,ds-\mathbf 1_{y\gt1}\Bigr|\le\frac{y^c}{\pi T|\log y|}. 2 π i 1 ∫ c − i T c + i T s y s d s − 1 y > 1 ≤ π T ∣ log y ∣ y c .
証明(補題 7.10) y > 1 y\gt1 y > 1 のとき、U > 0 U\gt0 U > 0 について頂点 c ± i T c\pm iT c ± i T 、− U ± i T -U\pm iT − U ± i T の長方形の周に留数定理(定理 5.1[留数定理と定積分への応用] )を使うと、内部の特異点は s = 0 s=0 s = 0 (留数 1 1 1 )だけです。水平な 2 辺の積分は、∣ y s / s ∣ ≤ y σ / T |y^s/s|\le y^\sigma/T ∣ y s / s ∣ ≤ y σ / T から合わせて 2 T ∫ − ∞ c y σ d σ = 2 y c T log y \frac2T\int_{-\infty}^cy^\sigma d\sigma=\frac{2y^c}{T\log y} T 2 ∫ − ∞ c y σ d σ = T l o g y 2 y c 以下、左の縦の辺は 2 T y − U U → 0 \frac{2Ty^{-U}}U\to0 U 2 T y − U → 0 (U → ∞ U\to\infty U → ∞ )なので、2 π 2\pi 2 π で割って主張です。y < 1 y\lt1 y < 1 のときは右側の長方形(頂点 c ± i T c\pm iT c ± i T 、U ± i T U\pm iT U ± i T )を使うと内部に特異点がなく、同様に 2 T ∫ c ∞ y σ d σ = 2 y c T ∣ log y ∣ \frac{2}{T}\int_c^\infty y^\sigma d\sigma=\frac{2y^c}{T|\log y|} T 2 ∫ c ∞ y σ d σ = T ∣ l o g y ∣ 2 y c で抑えられます。
∎
補題 7.11 (リーマン予想の下での ψ の誤差 ) ✓
リーマン予想(公理 5.3 )を仮定する。定数 C C C があって、x ≥ 2 x\ge2 x ≥ 2 で ∣ ψ ( x ) − x ∣ ≤ C x log 2 x |\psi(x)-x|\le C\sqrt x\log^2x ∣ ψ ( x ) − x ∣ ≤ C x log 2 x 。
証明(補題 7.11) ψ \psi ψ は [ N , N + 1 ) [N,N+1) [ N , N + 1 ) (N N N は整数)で定数なので、x = N + 1 2 x=N+\frac12 x = N + 2 1 (N ≥ 5 N\ge5 N ≥ 5 )の場合に示せば十分です(一般の x x x では ∣ x − ( ⌊ x ⌋ + 1 2 ) ∣ ≤ 1 2 |x-(\lfloor x\rfloor+\frac12)|\le\frac12 ∣ x − (⌊ x ⌋ + 2 1 ) ∣ ≤ 2 1 、小さい x x x は C C C で吸収)。c = 1 + 1 / log x c=1+1/\log x c = 1 + 1/ log x とおきます。
零点の位置. 定義 5.2 と仮定から、σ ≥ 1 4 \sigma\ge\frac14 σ ≥ 4 1 にある ζ \zeta ζ の零点はすべて臨界線上の ρ = 1 2 + i γ \rho=\frac12+i\gamma ρ = 2 1 + iγ です(σ > 1 \sigma\gt1 σ > 1 には零点がない)。ζ ( s ˉ ) = ζ ( s ) ‾ \zeta(\bar s)=\overline{\zeta(s)} ζ ( s ˉ ) = ζ ( s ) なので零点は実軸について対称です。
(1) Perron の公式. 線分 R e s = c \mathrm{Re}\,s=c Re s = c 上では ∑ n Λ ( n ) ∣ n − s ∣ ≤ ∑ Λ ( n ) n − c < ∞ \sum_n\Lambda(n)|n^{-s}|\le\sum\Lambda(n)n^{-c}\lt\infty ∑ n Λ ( n ) ∣ n − s ∣ ≤ ∑ Λ ( n ) n − c < ∞ なので、項別に積分して 補題 7.5 の 2 から
J ( x , T ) : = 1 2 π i ∫ c − i T c + i T ( − ζ ′ ζ ( s ) ) x s s d s = ∑ n Λ ( n ) 1 2 π i ∫ c − i T c + i T ( x / n ) s s d s . J(x,T):=\frac1{2\pi i}\int_{c-iT}^{c+iT}\Bigl(-\frac{\zeta'}\zeta(s)\Bigr)\frac{x^s}s\,ds=\sum_n\Lambda(n)\,\frac1{2\pi i}\int_{c-iT}^{c+iT}\frac{(x/n)^s}{s}ds . J ( x , T ) := 2 π i 1 ∫ c − i T c + i T ( − ζ ζ ′ ( s ) ) s x s d s = n ∑ Λ ( n ) 2 π i 1 ∫ c − i T c + i T s ( x / n ) s d s . 補題 7.10 から ∣ ψ ( x ) − J ( x , T ) ∣ ≤ ∑ n Λ ( n ) ( x / n ) c π T ∣ log ( x / n ) ∣ |\psi(x)-J(x,T)|\le\sum_n\Lambda(n)\frac{(x/n)^c}{\pi T|\log(x/n)|} ∣ ψ ( x ) − J ( x , T ) ∣ ≤ ∑ n Λ ( n ) π T ∣ l o g ( x / n ) ∣ ( x / n ) c です。n ≤ x / 2 n\le x/2 n ≤ x /2 または n ≥ 2 x n\ge2x n ≥ 2 x の項は ∣ log ( x / n ) ∣ ≥ log 2 |\log(x/n)|\ge\log2 ∣ log ( x / n ) ∣ ≥ log 2 、x c = e x x^c=ex x c = e x 、補題 7.5 の ∑ Λ ( n ) n − c ≤ 3 c log x \sum\Lambda(n)n^{-c}\le3c\log x ∑ Λ ( n ) n − c ≤ 3 c log x から合わせて 3 e c x log x π T log 2 ≤ 8 x log x T \frac{3ec\,x\log x}{\pi T\log2}\le\frac{8x\log x}{T} π T l o g 2 3 ec x l o g x ≤ T 8 x l o g x 以下です。x / 2 < n < 2 x x/2\lt n\lt2x x /2 < n < 2 x の項では Λ ( n ) ≤ log 2 x \Lambda(n)\le\log2x Λ ( n ) ≤ log 2 x 、( x / n ) c ≤ 4 (x/n)^c\le4 ( x / n ) c ≤ 4 、u = x / n ∈ ( 1 2 , 2 ) u=x/n\in(\frac12,2) u = x / n ∈ ( 2 1 , 2 ) について ∣ log u ∣ ≥ ∣ u − 1 ∣ / 2 ≥ ∣ x − n ∣ / ( 4 x ) |\log u|\ge|u-1|/2\ge|x-n|/(4x) ∣ log u ∣ ≥ ∣ u − 1∣/2 ≥ ∣ x − n ∣/ ( 4 x ) 、∣ x − n ∣ ≥ 1 2 |x-n|\ge\frac12 ∣ x − n ∣ ≥ 2 1 なので、合わせて 16 x log 2 x π T ∑ x / 2 < n < 2 x 1 ∣ x − n ∣ ≤ 16 x log 2 x π T ⋅ 2 ( 2 + log ( 4 x + 1 ) ) \frac{16x\log2x}{\pi T}\sum_{x/2\lt n\lt2x}\frac1{|x-n|}\le\frac{16x\log2x}{\pi T}\cdot2\bigl(2+\log(4x+1)\bigr) π T 16 x l o g 2 x ∑ x /2 < n < 2 x ∣ x − n ∣ 1 ≤ π T 16 x l o g 2 x ⋅ 2 ( 2 + log ( 4 x + 1 ) ) 以下です。よって ∣ ψ ( x ) − J ( x , T ) ∣ ≤ C 1 x log 2 x T |\psi(x)-J(x,T)|\le C_1\frac{x\log^2x}{T} ∣ ψ ( x ) − J ( x , T ) ∣ ≤ C 1 T x l o g 2 x です。
(2) 高さ T T T の選び方. 補題 7.9 を t = x + 1 2 t=x+\frac12 t = x + 2 1 で使うと、γ ∈ [ x , x + 1 ] \gamma\in[x,x+1] γ ∈ [ x , x + 1 ] の零点は Z ( x + 1 2 ) Z(x+\frac12) Z ( x + 2 1 ) に属し、その個数 K K K は 18 log ( x + 1 2 ) 18\log(x+\frac12) 18 log ( x + 2 1 ) 以下です。[ x , x + 1 ] [x,x+1] [ x , x + 1 ] を K + 1 K+1 K + 1 等分すると、内部に γ \gamma γ を含まない小区間があり、その中点を T T T とすると、すべての零点について ∣ T − γ ∣ ≥ 1 2 ( K + 1 ) ≥ 1 40 log x |T-\gamma|\ge\frac1{2(K+1)}\ge\frac1{40\log x} ∣ T − γ ∣ ≥ 2 ( K + 1 ) 1 ≥ 40 l o g x 1 (x ≥ 5 x\ge5 x ≥ 5 )。対称性から ∣ − T − γ ∣ |-T-\gamma| ∣ − T − γ ∣ も同じです。
(3) 積分路の移動. 頂点 c ± i T c\pm iT c ± i T 、1 4 ± i T \frac14\pm iT 4 1 ± i T の長方形の周に留数定理を使います。領域としては、∣ γ ∣ ∈ ( T , T + δ ) |\gamma|\in(T,T+\delta) ∣ γ ∣ ∈ ( T , T + δ ) の零点がないように小さい δ > 0 \delta\gt0 δ > 0 を取った凸領域 { 0 < σ < 3 , ∣ t ∣ < T + δ } \{0\lt\sigma\lt3,\ |t|\lt T+\delta\} { 0 < σ < 3 , ∣ t ∣ < T + δ } を取ります。周の上に零点はなく、内部の特異点は s = 1 s=1 s = 1 (定理 4.3 から ζ ( s ) = h ( s ) / ( s − 1 ) \zeta(s)=h(s)/(s-1) ζ ( s ) = h ( s ) / ( s − 1 ) 、h ( 1 ) = 1 h(1)=1 h ( 1 ) = 1 と書けて − ζ ′ / ζ = 1 s − 1 − h ′ / h -\zeta'/\zeta=\frac1{s-1}-h'/h − ζ ′ / ζ = s − 1 1 − h ′ / h なので、− ζ ′ ζ x s s -\frac{\zeta'}{\zeta}\frac{x^s}s − ζ ζ ′ s x s の留数は x x x )と ∣ γ ∣ < T |\gamma|\lt T ∣ γ ∣ < T の零点 ρ \rho ρ (補題 2.3[偏角の原理と Rouché の定理] から留数は − m ρ x ρ / ρ -m_\rho x^\rho/\rho − m ρ x ρ / ρ 、m ρ m_\rho m ρ は位数)です。よって
J ( x , T ) = x − ∑ ∣ γ ∣ < T x ρ ρ − 1 2 π i ( ∫ c + i T 1 4 + i T + ∫ 1 4 + i T 1 4 − i T + ∫ 1 4 − i T c − i T ) ( − ζ ′ ζ ( s ) ) x s s d s J(x,T)=x-\sum_{|\gamma|\lt T}\frac{x^\rho}{\rho}-\frac1{2\pi i}\Bigl(\int_{c+iT}^{\frac14+iT}+\int_{\frac14+iT}^{\frac14-iT}+\int_{\frac14-iT}^{c-iT}\Bigr)\Bigl(-\frac{\zeta'}{\zeta}(s)\Bigr)\frac{x^s}{s}ds J ( x , T ) = x − ∣ γ ∣ < T ∑ ρ x ρ − 2 π i 1 ( ∫ c + i T 4 1 + i T + ∫ 4 1 + i T 4 1 − i T + ∫ 4 1 − i T c − i T ) ( − ζ ζ ′ ( s ) ) s x s d s (和は重複込み)。
(4) 評価. 水平な辺 s = σ ± i T s=\sigma\pm iT s = σ ± i T (1 4 ≤ σ ≤ c ≤ 2 \frac14\le\sigma\le c\le2 4 1 ≤ σ ≤ c ≤ 2 )では ∣ s − ( 2 ± i T ) ∣ = 2 − σ ≤ 7 4 |s-(2\pm iT)|=2-\sigma\le\frac74 ∣ s − ( 2 ± i T ) ∣ = 2 − σ ≤ 4 7 、∣ s − ρ ∣ ≥ ∣ T ∓ γ ∣ ≥ 1 40 log x |s-\rho|\ge|T\mp\gamma|\ge\frac1{40\log x} ∣ s − ρ ∣ ≥ ∣ T ∓ γ ∣ ≥ 40 l o g x 1 なので、補題 7.9 から ∣ ζ ′ / ζ ( s ) ∣ ≤ 18 log T ⋅ 40 log x + 1900 log T ≤ C 2 log 2 x |\zeta'/\zeta(s)|\le18\log T\cdot40\log x+1900\log T\le C_2\log^2x ∣ ζ ′ / ζ ( s ) ∣ ≤ 18 log T ⋅ 40 log x + 1900 log T ≤ C 2 log 2 x (T ≤ x + 1 T\le x+1 T ≤ x + 1 )。∣ x s / s ∣ ≤ x σ / T |x^s/s|\le x^\sigma/T ∣ x s / s ∣ ≤ x σ / T 、∫ 1 / 4 c x σ d σ ≤ x c / log x = e x / log x \int_{1/4}^cx^\sigma d\sigma\le x^c/\log x=ex/\log x ∫ 1/4 c x σ d σ ≤ x c / log x = e x / log x なので、水平な辺の寄与は C 3 x log x / T C_3x\log x/T C 3 x log x / T 以下です。縦の辺 s = 1 4 + i t s=\frac14+it s = 4 1 + i t では、∣ t ∣ ≥ 5 |t|\ge5 ∣ t ∣ ≥ 5 なら ∣ s − ρ ∣ ≥ 1 4 |s-\rho|\ge\frac14 ∣ s − ρ ∣ ≥ 4 1 から ∣ ζ ′ / ζ ( s ) ∣ ≤ 4 ⋅ 18 log ∣ t ∣ + 1900 log ∣ t ∣ ≤ 2000 log ∣ t ∣ |\zeta'/\zeta(s)|\le4\cdot18\log|t|+1900\log|t|\le2000\log|t| ∣ ζ ′ / ζ ( s ) ∣ ≤ 4 ⋅ 18 log ∣ t ∣ + 1900 log ∣ t ∣ ≤ 2000 log ∣ t ∣ 、∣ t ∣ ≤ 5 |t|\le5 ∣ t ∣ ≤ 5 では零点も極もない有界閉線分上の連続関数として有界です。よって寄与は x 1 / 4 ( C 4 + ∫ 5 T 4000 log t t d t ) ≤ C 5 x 1 / 4 log 2 x x^{1/4}\bigl(C_4+\int_5^T\frac{4000\log t}{t}dt\bigr)\le C_5x^{1/4}\log^2x x 1/4 ( C 4 + ∫ 5 T t 4000 l o g t d t ) ≤ C 5 x 1/4 log 2 x です。零点の和は、∣ x ρ ∣ = x |x^\rho|=\sqrt x ∣ x ρ ∣ = x 、∣ ρ ∣ ≥ max ( 1 2 , ∣ γ ∣ ) |\rho|\ge\max(\frac12,|\gamma|) ∣ ρ ∣ ≥ max ( 2 1 , ∣ γ ∣ ) と、k ≥ 5 k\ge5 k ≥ 5 について γ ∈ [ k , k + 1 ] \gamma\in[k,k+1] γ ∈ [ k , k + 1 ] の零点が Z ( k + 1 2 ) Z(k+\frac12) Z ( k + 2 1 ) に属し 18 log ( k + 1 ) 18\log(k+1) 18 log ( k + 1 ) 個以下であること(負の側も対称)、∣ γ ∣ ≤ 5 |\gamma|\le5 ∣ γ ∣ ≤ 5 の零点が有限個であることから
∑ ∣ γ ∣ < T ∣ x ρ ∣ ∣ ρ ∣ ≤ x ( C 6 + 2 ∑ 5 ≤ k ≤ T 18 log ( k + 1 ) k ) ≤ C 7 x log 2 x . \sum_{|\gamma|\lt T}\frac{|x^\rho|}{|\rho|}\le\sqrt x\Bigl(C_6+2\sum_{5\le k\le T}\frac{18\log(k+1)}{k}\Bigr)\le C_7\sqrt x\log^2x . ∣ γ ∣ < T ∑ ∣ ρ ∣ ∣ x ρ ∣ ≤ x ( C 6 + 2 5 ≤ k ≤ T ∑ k 18 log ( k + 1 ) ) ≤ C 7 x log 2 x . (1)〜(4) を合わせ、x ≤ T ≤ x + 1 x\le T\le x+1 x ≤ T ≤ x + 1 を使えば ∣ ψ ( x ) − x ∣ ≤ C x log 2 x |\psi(x)-x|\le C\sqrt x\log^2x ∣ ψ ( x ) − x ∣ ≤ C x log 2 x です。
∎
補題 7.12 (誤差の評価から零点のない領域へ ) ✓
1 2 ≤ α < 1 \frac12\le\alpha\lt1 2 1 ≤ α < 1 とし、任意の ε > 0 \varepsilon\gt0 ε > 0 について定数 C ε C_\varepsilon C ε があって x ≥ 1 x\ge1 x ≥ 1 で ∣ ψ ( x ) − x ∣ ≤ C ε x α + ε |\psi(x)-x|\le C_\varepsilon x^{\alpha+\varepsilon} ∣ ψ ( x ) − x ∣ ≤ C ε x α + ε とする。このとき ζ \zeta ζ は R e s > α \mathrm{Re}\,s\gt\alpha Re s > α に零点をもたない。とくに α = 1 2 \alpha=\frac12 α = 2 1 ならリーマン予想が成り立つ。
証明(補題 7.12) 補題 7.5 の 2 と s ∫ 1 ∞ x ⋅ x − s − 1 d x = s s − 1 s\int_1^\infty x\cdot x^{-s-1}dx=\frac s{s-1} s ∫ 1 ∞ x ⋅ x − s − 1 d x = s − 1 s から、σ > 1 \sigma\gt1 σ > 1 で
− ζ ′ ( s ) ζ ( s ) − s s − 1 = s ∫ 1 ∞ ( ψ ( x ) − x ) x − s − 1 d x = : G ( s ) . -\frac{\zeta'(s)}{\zeta(s)}-\frac{s}{s-1}=s\int_1^\infty(\psi(x)-x)x^{-s-1}dx=:G(s). − ζ ( s ) ζ ′ ( s ) − s − 1 s = s ∫ 1 ∞ ( ψ ( x ) − x ) x − s − 1 d x =: G ( s ) . σ ≥ α + 2 ε \sigma\ge\alpha+2\varepsilon σ ≥ α + 2 ε では被積分関数の絶対値は C ε x − 1 − ε C_\varepsilon x^{-1-\varepsilon} C ε x − 1 − ε 以下なので、補題 7.6 の証明と同じく G G G は σ > α \sigma\gt\alpha σ > α で正則です。左辺は σ > α \sigma\gt\alpha σ > α で有理型で、その極を除いた(連結な)領域で G G G と一致します(一致の定理)。ζ ( ρ ) = 0 \zeta(\rho)=0 ζ ( ρ ) = 0 、R e ρ > α \mathrm{Re}\,\rho\gt\alpha Re ρ > α なら ρ ≠ 1 \rho\ne1 ρ = 1 で、補題 2.3[偏角の原理と Rouché の定理] から − ζ ′ / ζ -\zeta'/\zeta − ζ ′ / ζ は ρ \rho ρ で留数 − m ≠ 0 -m\ne0 − m = 0 の極をもち、s s − 1 \frac s{s-1} s − 1 s は ρ \rho ρ の近くで正則なので、左辺は ρ \rho ρ の近くで有界になりません。これは G G G が ρ \rho ρ で正則であることに矛盾します。
α = 1 2 \alpha=\frac12 α = 2 1 とします。非自明零点 ρ \rho ρ は 0 ≤ R e ρ ≤ 1 0\le\mathrm{Re}\,\rho\le1 0 ≤ Re ρ ≤ 1 を満たし(定義 5.2 )、いま示したことから R e ρ ≤ 1 2 \mathrm{Re}\,\rho\le\frac12 Re ρ ≤ 2 1 です。R e ρ < 1 2 \mathrm{Re}\,\rho\lt\frac12 Re ρ < 2 1 と仮定します。ρ ≠ 0 \rho\ne0 ρ = 0 (ζ ( 0 ) = − 1 2 \zeta(0)=-\frac12 ζ ( 0 ) = − 2 1 、命題 4.7 )なので ρ / 2 \rho/2 ρ /2 と ( 1 − ρ ) / 2 (1-\rho)/2 ( 1 − ρ ) /2 は Γ \Gamma Γ の極(0 , − 1 , − 2 , … 0,-1,-2,\ldots 0 , − 1 , − 2 , … 、定理 6.2[整関数とガンマ関数] )でなく、Γ \Gamma Γ は零点をもたない(定理 6.4[整関数とガンマ関数] )ので、ξ ( ρ ) = 0 \xi(\rho)=0 ξ ( ρ ) = 0 で、定理 4.5 の ξ ( 1 − ρ ) = ξ ( ρ ) = 0 \xi(1-\rho)=\xi(\rho)=0 ξ ( 1 − ρ ) = ξ ( ρ ) = 0 から ζ ( 1 − ρ ) = 0 \zeta(1-\rho)=0 ζ ( 1 − ρ ) = 0 です。ところが R e ( 1 − ρ ) > 1 2 \mathrm{Re}(1-\rho)\gt\frac12 Re ( 1 − ρ ) > 2 1 なので前半に矛盾します。よってすべての非自明零点は R e ρ = 1 2 \mathrm{Re}\,\rho=\frac12 Re ρ = 2 1 を満たします。
∎ 執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 5 回(Claude Code)・最終 2026-10-12