複素数 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 を含めない)とします。
Definition 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 が級数表示の限界です。
Theorem 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 の意味であり、この極限は存在する。
Proof(Theorem 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 とします。絶対収束する級数は有限個であれば自由に掛け合わせて項別に並べ替えられるので
∏ 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 以下であるような自然数」の集合です。二つ目の等号で素因数分解の一意性 (Theorem 4.2[Primes and the Prime Number Theorem] )を使いました。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 ζ が右半平面で消えないことが出ます。
Corollary 3.2 (右半平面での非零性 )
Re s > 1 \operatorname{Re} s > 1 Re s > 1 ならば ζ ( s ) ≠ 0 \zeta(s) \neq 0 ζ ( s ) = 0 である。
Proof(Corollary 3.2) σ > 1 \sigma > 1 σ > 1 とします。∑ n ∣ n − s ∣ \sum_n |n^{-s}| ∑ n ∣ n − s ∣ が収束するので並べ替えが自由にでき、任意の素数 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 です。
∎ Proof(Corollary 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 しかないと仮定します。Theorem 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 ) → ∞ に矛盾します。
∎ 同じ道具をもう少し押すと、ユークリッドの証明(ユークリッドの定理(Theorem 3.2)[Primes and the Prime Number Theorem] )では届かない「素数の密度」の情報が出ます。
Proposition 3.4 (素数の逆数和の発散 )
∑ p : 素数 1 p = ∞ \displaystyle\sum_{p:\ \text{素数}} \frac{1}{p} = \infty p : 素数 ∑ p 1 = ∞ である。
Proof(Proposition 3.4) s > 1 s > 1 s > 1 を実数とします。Theorem 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 を得ます(二重級数はすべて正項なので順序交換は自由です)。剰余項 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. Corollary 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 は発散する。これは「素数は平方数よりずっと多い」という定量的な情報で、素数の無限性よりはるかに強い主張です。オイラー積が素数の分布 を掴んでいることの最初の証拠と言えます(この発散は Theorem 5.2[Primes and the Prime Number Theorem] としても扱いました)。
Example 3.5 (バーゼル問題と「二数が互いに素である確率」 )
バーゼル問題の答え ζ ( 2 ) = π 2 / 6 = 1.644934 … \zeta(2) = \pi^{2}/6 = 1.644934\ldots ζ ( 2 ) = π 2 /6 = 1.644934 … を認めます。Theorem 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 ポイント未満でした。素数の「独立性」という素朴な直観が、オイラー積という等式によって正確な形を与えられています。
級数表示 Definition 2.1 は Re s > 1 \operatorname{Re} s > 1 Re s > 1 でしか意味を持ちません。しかし素数の情報が濃く現れるのは、むしろその外側です。そこで解析接続 という操作が要ります。
Definition 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 ζ に一致し、より広い範囲で収束する式を一つ見つければ十分です。交代級数がその役を果たします。
Lemma 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 の各コンパクト部分集合上で一様収束し、そこで正則関数を定める。
Proof(Lemma 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 ) 自身も同じ値に収束します。
∎ Theorem 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. Proof(Theorem 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 ) と定めます。Lemma 4.2 より分子は Re s > 0 \operatorname{Re} s > 0 Re s > 0 で正則、分母は整関数ですから、右辺は分母の零点を除いて正則で、そこでは有理型です。第 1 段よりこの関数は Re s > 1 \operatorname{Re} s > 1 Re s > 1 で元の ζ \zeta ζ に一致します。一意性は Definition 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. ∎ 同じ手口をさらに進める(あるいはテータ関数のメリン変換を使う)と、ζ \zeta ζ は複素平面全体へ接続され、次の顕著な対称性が現れます。
Theorem 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 ) を満たす。
ξ ( 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 で成り立ち、一致の定理で全平面へ伝わります。したがって零点は実軸に関しても対称 です。零点は上下左右の二重の対称性を持って配置されているわけです。
Example 4.7 (関数等式から特殊値を計算する )
s = − 1 s = -1 s = − 1 を Theorem 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 、また Theorem 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 を除外する根拠になります。
Proposition 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 ζ の零点はこれらに限る。
Proof(Proposition 5.1) σ = Re s < 0 \sigma = \operatorname{Re} s < 0 σ = Re s < 0 とし、Theorem 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 なので、Corollary 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 になります。
∎ Definition 5.2 (非自明零点と臨界帯 )
ζ \zeta ζ の零点のうち s = − 2 , − 4 , − 6 , … s = -2, -4, -6, \ldots s = − 2 , − 4 , − 6 , … 以外のものを非自明零点 と呼ぶ。Corollary 3.2 (Re s > 1 \operatorname{Re} s > 1 Re s > 1 に零点なし)と Proposition 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 非自明零点 (実軸に関して対称)
臨界帯と零点の配置(虚部の目盛りは模式的) Axiom 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 である。
Example 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 の重みを付けたチェビシェフ関数を使うと式が最も簡単になります。
Definition 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 ) とおく。また対数積分(Definition 6.1[Primes and the Prime Number Theorem] )を主値積分
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 の程度に揃います。リーマン予想とは「誤差項が可能な限り小さい」という主張なのです。
flowchart TD
A["オイラー積: zeta(s) = prod (1 - p^-s)^-1"] --> B["対数微分: -zeta'/zeta(s) = sum Lambda(n) n^-s"]
B --> C["ペロンの公式で複素積分に直す"]
C --> D["留数定理: s=1 の極 と 非自明零点 rho の寄与を拾う"]
D --> E["明示公式: psi(x) = x - sum_rho x^rho / rho - ..."]
E --> F["|x^rho| = x^(Re rho): 誤差の指数 = 零点の実部"]
F --> G["リーマン予想 (Re rho = 1/2) は誤差 O(sqrt(x)) と同値"] ゼータ関数から素数分布へ至る論理の流れ この直観を厳密にしたのがフォン・コッホの定理です。
Theorem 6.3 (フォン・コッホの同値定理 )
次の二つは同値である。
リーマン予想(Axiom 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 + ε ) 」とも同値である。
Example 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 よりむしろ緩やかです。Theorem 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 の位置は今も特定されておらず、上界(スキューズ数と呼ばれる系列の評価)が知られているだけです。有限の数値実験がどれほど説得的に見えても証明にはならない、という数論の教訓がここにあります。
リーマン予想が素数以外にも波及することを示すため、同値な言明をいくつか挙げます。見た目の分野がまったく違うのに、すべて臨界線上の零点という一点に帰着します。
同値な言明 内容 出典 素数計数の誤差 π ( 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 ) \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 文として書き下すための道具にもなります(Example 8.2[ゲーデルの不完全性定理] )。
さらに、この構図は ζ \zeta ζ に限りません。数論に現れる多くの対象——ディリクレの L L L 関数、代数体のデデキント・ゼータ関数、そして 楕円曲線とモジュラー形式 で扱う楕円曲線の L L L 関数(Definition 4.4[楕円曲線とモジュラー形式] のトレース a p a_p a p から作られます)——にも同型の関数等式と「一般化リーマン予想」があります。同じ形の予想が広い範囲で立つこと自体が、零点の位置が偶然ではないことを示唆しています。リーマン予想は、素数分布という一つの問題であると同時に、L L L 関数という一族全体の設計図についての問いでもあるのです。
Exercise 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 の値を求めてください。
Solution 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 … です。
Exercise 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 ) を示してください。
Solution σ > 1 \sigma > 1 σ > 1 とします。Theorem 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 ) . 最後の等号で Theorem 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 なので適用可能で、Corollary 3.2 より ζ ( 2 s ) ≠ 0 \zeta(2s) \ne 0 ζ ( 2 s ) = 0 なので割り算も正当です)。
Exercise 7.3 標準
Theorem 4.5 と ζ ( 4 ) = π 4 / 90 \zeta(4) = \pi^{4}/90 ζ ( 4 ) = π 4 /90 を使って ζ ( − 3 ) \zeta(-3) ζ ( − 3 ) を求めてください。
Solution 関数等式
ζ ( 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 で、上の計算と一致します。
Exercise 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 Λ は Definition 6.1 のフォン・マンゴルト関数)。これが Remark 6.2 の明示公式の出発点です。
Solution σ 0 > 1 \sigma_0 > 1 σ 0 > 1 を固定し、半平面 Re s ≥ σ 0 \operatorname{Re} s \ge \sigma_0 Re s ≥ σ 0 で考えます。Corollary 3.2 より ζ ( s ) ≠ 0 \zeta(s)\ne0 ζ ( s ) = 0 なので ζ ′ / ζ \zeta'/\zeta ζ ′ / ζ はここで正則です。
Proposition 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. — Theorem 6.3 の原論文。
J. B. Conrey, “The Riemann Hypothesis”, Notices of the American Mathematical Society 50 (2003), 341–353. — 現代的な研究状況の概観。
明示公式の各項を波として読む方法を補足します。 Remark 6.2 の和に現れる零点は共役対 ρ = 1 2 + i γ \rho = \tfrac12 + i\gamma ρ = 2 1 + iγ と ρ ˉ = 1 2 − i γ \bar\rho = \tfrac12 - i\gamma ρ ˉ = 2 1 − iγ で現れます(Theorem 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 β に比例して他より速く育ち、いずれすべての波を圧倒します。リーマン予想とは「どの波も突出しない」、言い換えれば素数の分布が主要項のまわりで可能な限り整然と揺れている、という主張なのです。