リーマンゼータ関数の特殊値

リーマンゼータ関数の特殊値（リーマンゼータかんすうのとくしゅち、英: Particular values of the Riemann zeta function）とは、数学におけるリーマンゼータ関数（英: Riemann zeta function）に整数を代入した際の値のことをいう。これはリーマンゼータ値（英: Riemann zeta value）とも呼ばれる。 ゼータ関数は複素解析に頻繁に登場する特殊関数であるが、解析的整数論においても重要な関数である。ゼータ関数は、実部が 1 より真に大きい複素数 s と自然数 n に対して、 ζ ( s ) = ∑ n = 1 ∞ 1 n s = 1 + 1 2 s + 1 3 s + 1 4 s + ⋯ {\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}=1+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+{\frac {1}{4^{s}}}+\cdots } で定義される関数 ζ のことをいい、例えば s = 2 とすると、 ζ ( 2 ) = ∑ n = 1 ∞ 1 n 2 = 1 + 1 2 2 + 1 3 2 + 1 4 2 + ⋯ {\displaystyle \zeta (2)=\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}=1+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+\cdots } のような級数が提供される。特に、整数引数に対してゼータ関数がとる値についてはこの例も含めすべて実数値をもち、さらに数値計算に効率のよい公式が存在する。この記事では、これらの公式を値の表とともに列挙し、その微分と整数引数でのゼータ関数からなる級数も記述する。 ゼータ関数は、s = 1 における一位の極を除き解析接続によって複素平面全体に拡張される。しかしながら、上の定義式は解析接続された s に対しては無効であり、直に計算を試みると対応する和が発散する。例えば、ゼータ関数において s = −1 のとき、 ζ ( − 1 ) = − 1 12 {\displaystyle \zeta (-1)=-{\frac {1}{12}}} となるが、これを上の定義式で計算すると、 ζ ( − 1 ) = lim N → ∞ ∑ n = 1 N n = lim N → ∞ N ( N + 1 ) 2 = ∞ {\displaystyle \zeta (-1)=\lim _{N\to \infty }\sum _{n=1}^{N}n=\lim _{N\to \infty }{\frac {N\,(N+1)}{2}}=\infty } となって発散級数となる。以下に列挙するゼータ関数の特殊値は負の偶数に対する特殊値も含み、これは恒等的に ζ (s) = 0 であり、いわゆる自明な零点となる。 1644年、イタリアのピエトロ・メンゴリ（イタリア語版、ドイツ語版）によって以下の問題が提起された。この問題は、解決に挑んだ数学者の多くがバーゼルの生まれであったことから、バーゼル問題と呼ばれる。 バーゼル問題 ― 以下の級数： ∑ n = 1 ∞ 1 n 2 = 1 + 1 2 2 + 1 3 2 + 1 4 2 + ⋯ {\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n^{2}}}=1+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+\cdots } は収束するか。収束するならばその値はいくつか。 バーゼル問題は、スイスのレオンハルト・オイラーによって初めて解決された。オイラーは、三角関数のテイラー級数およびその無限乗積の x2 の項の展開係数を比較することで、 − 1 3 ! = − 1 π 2 ∑ n = 1 ∞ 1 n 2 {\displaystyle -{\frac {1}{3!}}=-{\frac {1}{\pi ^{2}}}\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}} となることから、 ζ ( 2 ) = ∑ n = 1 ∞ 1 n 2 = π 2 6 {\displaystyle \zeta (2)=\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}={\frac {\pi ^{2}}{6}}} が成り立つことを示した。さらにオイラーの研究はバーゼル問題にとどまることはなく、より一般の場合の研究に努め、任意の自然数 n に対して、 ζ ( 2 n ) = ( − 1 ) n + 1 ( 2 π ) 2 n B 2 n 2 ( 2 n ) ! {\displaystyle \zeta (2n)={\frac {(-1)^{n+1}\,(2\pi )^{2n}\,B_{2n}}{2\,(2n)!}}} が成り立つことも示した。ただし、ここで B2n は 2n 番目のベルヌーイ数である。 正の偶数に対する特殊値の証明 — 余接関数の指数関数による定義から、 π z cot ⁡ π z = i π z e i π z + e − i π z e i π z − e − i π z {\displaystyle \pi z\cot \pi z=i\pi z\,{\frac {e^{i\pi z}+e^{-i\pi z}}{e^{i\pi z}-e^{-i\pi z}}}} である。ただし、ここで i は虚数単位である。ここで右辺は、 i π z e i π z + e − i π z e i π z − e − i π z = i π z 2 i π z e 2 i π z − 1 {\displaystyle i\pi z\,{\frac {e^{i\pi z}+e^{-i\pi z}}{e^{i\pi z}-e^{-i\pi z}}}=i\pi z\,{\frac {2i\pi z}{e^{2i\pi z}-1}}} であり、ベルヌーイ数の定義から、 π z cot ⁡ π z = i π z + ∑ n = 0 ∞ ( 2 i π z ) n B n n ! {\displaystyle \pi z\cot \pi z=i\pi z+\sum _{n=0}^{\infty }{\frac {(2i\pi z)^{n}\,B_{n}}{n!}}} となるので、ベルヌーイ数の特殊値を利用し、 π z cot ⁡ π z = 1 + ∑ n = 1 ∞ ( 2 i π z ) 2 n B 2 n ( 2 n ) ! {\displaystyle \pi z\cot \pi z=1+\sum _{n=1}^{\infty }{\frac {(2i\pi z)^{2n}\,B_{2n}}{(2n)!}}} 一方、正弦関数の無限乗積の対数は、 log ⁡ sin ⁡ π z = log ⁡ π z + ∑ k = 1 ∞ log ( 1 − z 2 k 2 ) {\displaystyle \log \sin \pi z=\log \pi z+\sum _{k=1}^{\infty }\log \!{\biggl (}1-{\frac {z^{2}}{k^{2}}}{\biggr )}} であり、両辺を微分して z を乗じると、 π z cot ⁡ π z = 1 − ∑ k = 1 ∞ 1 1 − z 2 / k 2 2 z 2 k 2 {\displaystyle \pi z\cot \pi z=1-\sum _{k=1}^{\infty }{\frac {1}{1-z^{2}\!/k^{2}}}{\frac {2z^{2}}{k^{2}}}} となるが、これをテイラー展開して整理すると、 π z cot ⁡ π z = 1 − 2 ∑ n = 1 ∞ ζ ( 2 n ) z 2 n {\displaystyle \pi z\cot \pi z=1-2\sum _{n=1}^{\infty }\zeta (2n)\,z^{2n}} それぞれで導いた πz cot πz を比較して、 ζ ( 2 n ) = ( − 1 ) n + 1 ( 2 π ) 2 n B 2 n 2 ( 2 n ) ! {\displaystyle \zeta (2n)={\frac {(-1)^{n+1}\,(2\pi )^{2n}\,B_{2n}}{2\,(2n)!}}} この公式により、正の偶数に対する特殊値を容易く計算することができる。しかるに n = 1 から小さい順に n = 10 まで計算してみると、 ζ ( 2 ) = ∑ n = 1 ∞ 1 n 2 = π 2 6 {\displaystyle \zeta (2)=\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}={\frac {\pi ^{2}}{6}}} ζ ( 4 ) = ∑ n = 1 ∞ 1 n 4 = π 4 90 {\displaystyle \zeta (4)=\sum _{n=1}^{\infty }{\frac {1}{n^{4}}}={\frac {\pi ^{4}}{90}}} ζ ( 6 ) = ∑ n = 1 ∞ 1 n 6 = π 6 945 {\displaystyle \zeta (6)=\sum _{n=1}^{\infty }{\frac {1}{n^{6}}}={\frac {\pi ^{6}}{945}}} ζ ( 8 ) = ∑ n = 1 ∞ 1 n 8 = π 8 9450 {\displaystyle \zeta (8)=\sum _{n=1}^{\infty }{\frac {1}{n^{8}}}={\frac {\pi ^{8}}{9450}}} ζ ( 10 ) = ∑ n = 1 ∞ 1 n 10 = π 10 93555 {\displaystyle \zeta (10)=\sum _{n=1}^{\infty }{\frac {1}{n^{10}}}={\frac {\pi ^{10}}{93555}}} ζ ( 12 ) = ∑ n = 1 ∞ 1 n 12 = 691 π 12 638512875 {\displaystyle \zeta (12)=\sum _{n=1}^{\infty }{\frac {1}{n^{12}}}={\frac {691\pi ^{12}}{638512875}}} ζ ( 14 ) = ∑ n = 1 ∞ 1 n 14 = 2 π 14 18243225 {\displaystyle \zeta (14)=\sum _{n=1}^{\infty }{\frac {1}{n^{14}}}={\frac {2\pi ^{14}}{18243225}}} ζ ( 16 ) = ∑ n = 1 ∞ 1 n 16 = 3617 π 10 325641566250 {\displaystyle \zeta (16)=\sum _{n=1}^{\infty }{\frac {1}{n^{16}}}={\frac {3617\pi ^{10}}{325641566250}}} ζ ( 18 ) = ∑ n = 1 ∞ 1 n 18 = 43867 π 10 38979295480125 {\displaystyle \zeta (18)=\sum _{n=1}^{\infty }{\frac {1}{n^{18}}}={\frac {43867\pi ^{10}}{38979295480125}}} ζ ( 20 ) = ∑ n = 1 ∞ 1 n 20 = 173611 π 10 1531329465290625 {\displaystyle \zeta (20)=\sum _{n=1}^{\infty }{\frac {1}{n^{20}}}={\frac {173611\pi ^{10}}{1531329465290625}}} となる。またその近似値は、以下の表に示す通りである。 この表からもわかるように、ゼータ関数は s → ∞ の極限で ζ (s) → 1 である。すなわち、 lim s → ∞ ζ ( s ) = 1 {\displaystyle \lim _{s\to \infty }\zeta (s)=1} である。また、自然数 n に対して、 a n ζ ( 2 n ) = π 2 n b n {\displaystyle a_{n}\,\zeta (2n)=\pi ^{2n}\,b_{n}} を満たすように an と bn を定める。ただし、ここで an と bn は任意の自然数 n に対して常に自然数をとるものとする。すると、このとき an と bn の n = 1 から n = 20 までの挙動は以下の表に示す通りである。 さらに cn = bn/an と定めると、偶数に対する特殊値はより簡単に、 ζ ( 2 n ) = ∑ m = 1 ∞ 1 m 2 n = c n π 2 n {\displaystyle \zeta (2n)=\sum _{m=1}^{\infty }{\frac {1}{m^{2n}}}=c_{n}\,\pi ^{2n}} とかくことができる。するとこのとき、 c n = ∑ k = 1 n − 1 ( − 1 ) k − 1 c n − k ( 2 k + 1 ) ! + ( − 1 ) n + 1 n ( 2 n + 1 ) ! {\displaystyle c_{n}=\sum _{k=1}^{n-1}{\frac {(-1)^{k-1}\,c_{n-k}}{(2k+1)!}}+{\frac {(-1)^{n+1}\,n}{(2n+1)!}}} なる漸化式が存在することがわかる。この漸化式は、ベルヌーイ数を効率的に求める漸化式に基づいている。また、特殊値の係数ではなく、ゼータ関数についての漸化式も存在する。余接関数の微分： d d z cot ⁡ z = − 1 − cot 2 ⁡ z {\displaystyle {\frac {\rm {d}}{{\rm {d}}z}}\cot z=-1-\cot ^{2}z} およびその部分分数分解による表現： cot ⁡ z = 1 z − 2 ∑ n = 1 ∞ ζ ( 2 n ) z 2 n − 1 π 2 n {\displaystyle \cot z={\frac {1}{z}}-2\sum _{n=1}^{\infty }{\frac {\zeta (2n)\,z^{2n-1}}{\pi ^{2n}}}} を用いれば、 ζ ( 2 n ) = 2 2 n + 1 ∑ k = 1 n − 1 ζ ( 2 k ) ζ ( 2 n − 2 k ) {\displaystyle \zeta (2n)={\frac {2}{2n+1}}\sum _{k=1}^{n-1}\zeta (2k)\,\zeta (2n-2k)} が容易に導かれる。ただし、ここで n > 1 である。 ゼータ関数は Re s > 1 なる複素数 s に対して定義される関数であるが、その定義式に s = 1 を代入すると、 ζ ( 1 ) = ∑ n = 1 ∞ 1 n = 1 + 1 2 + 1 3 + 1 4 + ⋯ {\displaystyle \zeta (1)=\sum _{n=1}^{\infty }{\frac {1}{n}}=1+{\frac {1}{2}}+{\frac {1}{3}}+{\frac {1}{4}}+\cdots } となって調和級数に一致する。調和級数は、古くにおいては収束すると考えられていたが、今日においては発散することが知られている。しかしこれはコーシーの主値は存在し、 lim ε → 0 ζ ( 1 + ε ) + ζ ( 1 − ε ) 2 = γ {\displaystyle \lim _{\varepsilon \to 0}{\frac {\zeta (1+\varepsilon )+\zeta (1-\varepsilon )}{2}}=\gamma } である。ただし、ここで γ はオイラーの定数である。 また、正の偶数に対する特殊値はベルヌーイ数を用いる形で一般化されたが、正の奇数に対する特殊値は簡潔な形で表すことができないことが知られている。例えば、ゼータ関数に s = 3 を代入した実数 ζ (3) はアペリーの定数として知られ、様々な積分表示や級数表示が発見されているものの、簡単な形で表すことができない。また ζ (3) は無理数であることがわかっている。この主張をアペリーの定理という。また、正の偶数に対する特殊値が常に無理数となることはその一般化された公式を見れば一目瞭然である一方、正の奇数に対する特殊値がすべて無理数であるかどうかは現在もまだわかっていないが、すべて無理数ではないかと予想されている。以下の表にその近似値を示す。 アペリーの定数をはじめとした正の奇数に対する特殊値には様々な積分表示や級数表示が与えられており、それらを計算する場合はゼータ関数の定義式を利用するのではなく、別の収束速度の速い公式を利用することが多い。 ζ ( 3 ) = 2 π 2 7 [ 3 2 − log ⁡ π + ∑ n = 1 ∞ ζ ( 2 n ) 2 2 n − 1 n ( 2 n + 1 ) ( 2 n + 2 ) ] {\displaystyle \zeta (3)={\frac {2\pi ^{2}}{7}}{\Biggl [}{\frac {3}{2}}-\log \pi +\sum _{n=1}^{\infty }{\frac {\zeta (2n)}{2^{2n-1}\,n\,(2n+1)(2n+2)}}{\Biggr ]}} ζ ( 3 ) = π 2 7 [ 1 − 4 ∑ n = 1 ∞ ζ ( 2 n ) 2 2 n ( 2 n + 1 ) ( 2 n + 2 ) ] {\displaystyle \zeta (3)={\frac {\pi ^{2}}{7}}{\Biggl [}1-4\sum _{n=1}^{\infty }{\frac {\zeta (2n)}{2^{2n}\,(2n+1)(2n+2)}}{\Biggr ]}} ζ ( 3 ) = 7 π 3 180 − 2 ∑ n = 1 ∞ 1 n 3 ( e 2 n π − 1 ) {\displaystyle \zeta (3)={\frac {7\pi ^{3}}{180}}-2\sum _{n=1}^{\infty }{\frac {1}{n^{3}\,(e^{2n\pi }-1)}}} ζ ( 3 ) = 5 2 ∑ n = 1 ∞ ( − 1 ) n − 1 ( n ! ) 2 n 3 ( 2 n ) ! {\displaystyle \zeta (3)={\frac {5}{2}}\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}\,(n!)^{2}}{n^{3}\,(2n)!}}} ζ ( 3 ) = 1 4 ∑ n = 1 ∞ ( n ! ) 4 ( 30 n − 11 ) n 3 [ ( 2 n ) ! ] 2 ( 2 n − 1 ) {\displaystyle \zeta (3)={\frac {1}{4}}\sum _{n=1}^{\infty }{\frac {(n!)^{4}\,(30n-11)}{n^{3}\,[(2n)!]^{2}\,(2n-1)}}} ζ ( 3 ) = 1 2 ∑ n = 1 ∞ ( − 1 ) n − 1 ( n ! ) 10 ( 205 n 2 − 160 n + 32 ) n 5 [ ( 2 n ) ! ] 5 {\displaystyle \zeta (3)={\frac {1}{2}}\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}\,(n!)^{10}\,(205n^{2}-160n+32)}{n^{5}\,[(2n)!]^{5}}}} ζ ( 3 ) = 14 ∑ n = 1 ∞ 1 n 3 sinh ⁡ n π − 11 2 ∑ n = 1 ∞ 1 n 3 ( e 2 n π − 1 ) − 7 2 ∑ n = 1 ∞ 1 n 3 ( e 2 n π + 1 ) {\displaystyle \zeta (3)=14\sum _{n=1}^{\infty }{\frac {1}{n^{3}\sinh n\pi }}-{\frac {11}{2}}\sum _{n=1}^{\infty }{\frac {1}{n^{3}\,(e^{2n\pi }-1)}}-{\frac {7}{2}}\sum _{n=1}^{\infty }{\frac {1}{n^{3}\,(e^{2n\pi }+1)}}} ζ ( 3 ) = 1 2 ∫ 0 ∞ x 2 e x − 1 d x {\displaystyle \zeta (3)={\frac {1}{2}}\int _{0}^{\infty }\!{\frac {x^{2}}{e^{x}-1}}\,{\rm {d}}x} ζ ( 3 ) = 2 π 2 log ⁡ 2 7 + 16 7 ∫ 0 π 2 x log ⁡ sin ⁡ x d x {\displaystyle \zeta (3)={\frac {2\pi ^{2}\log 2}{7}}+{\frac {16}{7}}\int _{0}^{\frac {\pi }{2}}x\log \sin x\,{\rm {d}}x} ζ ( 3 ) = − 1 2 ∫ 0 1 ∫ 0 1 log ⁡ x y 1 − x y d x d y {\displaystyle \zeta (3)=-{\frac {1}{2}}\int _{0}^{1}\!\!\int _{0}^{1}\!{\frac {\log xy}{1-xy}}\,{\rm {d}}x\,{\rm {d}}y} ζ ( 3 ) = ∫ 0 1 ∫ 0 1 ∫ 0 1 d x d y d z 1 − x y z {\displaystyle \zeta (3)=\int _{0}^{1}\!\!\int _{0}^{1}\!\!\int _{0}^{1}\!{\frac {{\rm {d}}x\,{\rm {d}}y\,{\rm {d}}z}{1-xyz}}} ζ ( 3 ) = Re ⁡ ∫ 0 1 ∫ 0 1 ∫ 0 1 d x d y d z f ( f ( f ( x ) y ) z ) w h e r e f ( x ) = 1 + 3 i 2 x {\displaystyle \zeta (3)=\operatorname {Re} \int _{0}^{1}\!\!\int _{0}^{1}\!\!\int _{0}^{1}\!{\frac {{\rm {d}}x\,{\rm {d}}y\,{\rm {d}}z}{f(f(f(x)\,y)\,z)}}\quad {\rm {where}}\,f(x)={\frac {1+{\sqrt {3}}i}{2}}\,x} ζ ( 5 ) = π 5 294 − 72 35 ∑ n = 1 ∞ 1 n 5 ( e 2 n π − 1 ) − 2 35 ∑ n = 1 ∞ 1 n 5 ( e 2 n π + 1 ) {\displaystyle \zeta (5)={\frac {\pi ^{5}}{294}}-{\frac {72}{35}}\sum _{n=1}^{\infty }{\frac {1}{n^{5}\,(e^{2n\pi }-1)}}-{\frac {2}{35}}\sum _{n=1}^{\infty }{\frac {1}{n^{5}\,(e^{2n\pi }+1)}}} ζ ( 5 ) = 12 ∑ n = 1 ∞ 1 n 5 sinh ⁡ n π − 39 20 ∑ n = 1 ∞ 1 n 5 ( e 2 n π − 1 ) + 1 20 ∑ n = 1 ∞ 1 n 5 ( e 2 n π + 1 ) {\displaystyle \zeta (5)=12\sum _{n=1}^{\infty }{\frac {1}{n^{5}\sinh n\pi }}-{\frac {39}{20}}\sum _{n=1}^{\infty }{\frac {1}{n^{5}\,(e^{2n\pi }-1)}}+{\frac {1}{20}}\sum _{n=1}^{\infty }{\frac {1}{n^{5}\,(e^{2n\pi }+1)}}} ζ ( 5 ) = 2 π 4 3 ( 2 5 − 1 ) [ log ⁡ 2 + ∑ n = 1 ∞ ( 2 n + 5 ) ζ ( 2 n ) 2 2 n ( n + 2 ) ( 2 n + 3 ) ] {\displaystyle \zeta (5)={\frac {2\pi ^{4}}{3\,(2^{5}-1)}}{\Biggl [}\log 2+\sum _{n=1}^{\infty }{\frac {(2n+5)\,\zeta (2n)}{2^{2n}\,(n+2)(2n+3)}}{\Biggr ]}} ζ ( 5 ) = 2 π 4 31 [ ∑ n = 1 ∞ ζ ( 2 n ) 2 2 n ( n + 2 ) ( 2 n + 1 ) − 3 ∑ n = 1 ∞ ζ ( 3 ) 2 2 n ( n + 1 ) ( 2 n + 3 ) ] {\displaystyle \zeta (5)={\frac {2\pi ^{4}}{31}}{\Biggl [}\sum _{n=1}^{\infty }{\frac {\zeta (2n)}{2^{2n}(n+2)(2n+1)}}-3\sum _{n=1}^{\infty }{\frac {\zeta (3)}{2^{2n}\,(n+1)(2n+3)}}{\Biggr ]}} ζ ( 5 ) = 1 24 ∫ 0 ∞ x 4 e x − 1 d x {\displaystyle \zeta (5)={\frac {1}{24}}\int _{0}^{\infty }\!{\frac {x^{4}}{e^{x}-1}}\,{\rm {d}}x} 正の奇数 n に対して、 S ± ( s ) = ∑ n = 1 ∞ 1 n s ( e 2 n π ± 1 ) {\displaystyle S_{\pm }(s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}\,(e^{2n\pi }\pm 1)}}} なる級数を定めるとき、ζ (3) や ζ (5) で見られたような一連の級数は次の形で定式化される。 A n ζ ( n ) − B n π n + C n S − ( n ) + D n S + ( n ) = 0 {\displaystyle A_{n}\,\zeta (n)-B_{n}\,\pi ^{n}+C_{n}\,S_{-}(n)+D_{n}\,S_{+}(n)=0} ただし、ここで An 、Bn 、Cn および Dn は、任意の正の奇数 n に対して常に自然数をとるものとする。ここでの Bn はベルヌーイ数とは異なる。すると、このとき An 、Bn 、Cn および Dn の n = 3 から n = 19 までの挙動は以下の表に示す通りである。 これらの整数はベルヌーイ数の和として表現することができる。任意の整数引数に対するゼータ関数の高速計算アルゴリズムは、アナトリー・カラツバ（英語版）によって与えられている。 ゼータ関数の定義式は、 ζ ( s ) = ∑ n = 1 ∞ 1 n s = 1 + 1 2 s + 1 3 s + 1 4 s + ⋯ {\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}=1+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+{\frac {1}{4^{s}}}+\cdots } であったが、このままでは負の整数に対する特殊値の計算を実行することができない。しかしながらゼータ関数には、 ζ ( s ) = − Γ ( 1 − s ) 2 π i ∮ C ( − z ) s − 1 e z − 1 d z {\displaystyle \zeta (s)=-{\frac {\varGamma (1-s)}{2\pi i}}\oint _{C}\!{\frac {(-z)^{s-1}}{e^{z}-1}}\,{\rm {d}}z} なる複素平面全体で定義された関数が存在する。この積分を利用することで、任意の自然数 n に対して、 ζ ( − n ) = ( − 1 ) n B n + 1 n + 1 {\displaystyle \zeta (-n)={\frac {(-1)^{n}\,B_{n+1}}{n+1}}} が成り立つことがわかる。 負の整数に対するリーマンゼータ値の証明 — ハンケルによるリーマンゼータ関数の積分表示： ζ ( s ) = − Γ ( 1 − s ) 2 π i ∮ C ( − z ) s − 1 e z − 1 d z {\displaystyle \zeta (s)=-{\frac {\varGamma (1-s)}{2\pi i}}\oint _{C}\!{\frac {(-z)^{s-1}}{e^{z}-1}}\,{\rm {d}}z} から始める。ただし、ここで s は自然数ではない複素数で、積分路 C はハンケルの積分路（英語版）である。ここで n を自然数として s = − n とすると、 ζ ( − n ) = i n ! 2 π ∮ C 1 ( − z ) n + 2 − z e z − 1 d z {\displaystyle \zeta (-n)={\frac {in!}{2\pi }}\oint _{C}\!{\frac {1}{(-z)^{n+2}}}{\frac {-z}{e^{z}-1}}\,{\rm {d}}z} となる。留数定理に基づいてこの複素積分を実行すれば、 ζ ( − n ) = − i n ! 2 π 2 π i Res ( 1 ( − z ) n + 2 − z e z − 1 , 0 ) {\displaystyle \zeta (-n)=-{\frac {in!}{2\pi }}\,2\pi i\operatorname {Res} \!{\biggl (}{\frac {1}{(-z)^{n+2}}}{\frac {-z}{e^{z}-1}},\,0{\biggr )}} となるのでベルヌーイ数の定義から、 ζ ( − n ) = n ! ( − 1 ) n Res ( ∑ k = 1 ∞ B k k ! z k − n − 2 , 0 ) {\displaystyle \zeta (-n)=n!\,(-1)^{n}\operatorname {Res} \!{\Biggl (}\sum _{k=1}^{\infty }{\frac {B_{k}}{k!}}\,z^{k-n-2},\,0{\Biggr )}} ここで、右辺の留数は、 Res ( ∑ k = 1 ∞ B k k ! z k − n − 2 , 0 ) = B n + 1 ( n + 1 ) ! {\displaystyle \operatorname {Res} \!{\Biggl (}\sum _{k=1}^{\infty }{\frac {B_{k}}{k!}}\,z^{k-n-2},\,0{\Biggr )}={\frac {B_{n+1}}{(n+1)!}}} であるから、 ζ ( − n ) = ( − 1 ) n B n + 1 n + 1 {\displaystyle \zeta (-n)={\frac {(-1)^{n}\,B_{n+1}}{n+1}}} この公式を利用することで、負の整数に対する特殊値を計算することができる。一般に負の偶数に対しては、 ζ ( − 2 n ) = 0 {\displaystyle \zeta (-2n)=0} が恒等的に成り立つ。これを自明な零点という。また負の奇数については n = 1 から小さい順に n = 9 まで計算してみると、 ζ ( − 1 ) = − 1 12 {\displaystyle \zeta (-1)=-{\frac {1}{12}}} ζ ( − 3 ) = 1 120 {\displaystyle \zeta (-3)={\frac {1}{120}}} ζ ( − 5 ) = − 1 252 {\displaystyle \zeta (-5)=-{\frac {1}{252}}} ζ ( − 7 ) = 1 240 {\displaystyle \zeta (-7)={\frac {1}{240}}} ζ ( − 9 ) = − 1 132 {\displaystyle \zeta (-9)=-{\frac {1}{132}}} ζ ( − 11 ) = 691 32760 {\displaystyle \zeta (-11)={\frac {691}{32760}}} ζ ( − 13 ) = − 1 12 {\displaystyle \zeta (-13)=-{\frac {1}{12}}} ζ ( − 15 ) = 3617 8160 {\displaystyle \zeta (-15)={\frac {3617}{8160}}} ζ ( − 17 ) = − 43867 14364 {\displaystyle \zeta (-17)=-{\frac {43867}{14364}}} ζ ( − 19 ) = 174611 6600 {\displaystyle \zeta (-19)={\frac {174611}{6600}}} となる。特に ζ (−1) はラマヌジャン総和法（英語版）に関連する。 ゼータ関数の負の偶数での微分係数は、 d ζ d s | s = − 2 n = ( − 1 ) n ( 2 n ) ! ζ ( 2 n + 1 ) 2 ( 2 π ) 2 n {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=-2n}={\frac {(-1)^{n}\,(2n)!\,\zeta (2n+1)}{2\,(2\pi )^{2n}}}} である。これはゼータ関数の一様収束性から項別に微分すれば簡単に示すことができる。この公式を用いて特殊値を計算すると、 d ζ d s | s = − 2 = − ζ ( 3 ) 4 π 2 {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=-2}=-{\frac {\zeta (3)}{4\pi ^{2}}}} d ζ d s | s = − 4 = 3 ζ ( 5 ) 4 π 4 {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=-4}={\frac {3\,\zeta (5)}{4\pi ^{4}}}} d ζ d s | s = − 6 = − 45 ζ ( 7 ) 8 π 6 {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=-6}=-{\frac {45\,\zeta (7)}{8\pi ^{6}}}} d ζ d s | s = − 8 = 315 ζ ( 9 ) 4 π 8 {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=-8}={\frac {315\,\zeta (9)}{4\pi ^{8}}}} となる。またこれ以外にも、 d ζ d s | s = − 1 = 1 12 − log ⁡ A {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=-1}={\frac {1}{12}}-\log A} d ζ d s | s = 0 = − 1 2 log ⁡ 2 π {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=0}=-{\frac {1}{2}}\log 2\pi } d ζ d s | s = 2 = π 2 6 ( γ + log ⁡ 2 − 12 log ⁡ A + log ⁡ π ) {\displaystyle {\frac {{\rm {d}}\zeta }{{\rm {d}}s}}{\Bigg \vert }_{s=2}={\frac {\pi ^{2}}{6}}\,(\gamma +\log 2-12\log A+\log \pi )} なる特殊値が存在する。ただし、ここで A はグレイシャー・キンケリンの定数、γ はオイラーの定数である。また、これらの特殊値の近似値は以下の表に示す通りである。