ベルヌーイ数

ベルヌーイ数 (ベルヌーイすう、英: Bernoulli number) は数論における基本的な係数を与える数列であり、もともと、連続する整数のべき乗和を定式化する際の展開係数として1713年にヤコブ・ベルヌーイが著書 Ars Conjectandi (推測術) にて導入したことからこの名称がついた。ベルヌーイ数は、べき乗和の展開係数にとどまらず、級数展開の係数や剰余項、リーマンゼータ関数においても登場する。また、ベルヌーイ数はすべてが有理数である。

定義

ベルヌーイ数 B_n は下に示すマクローリン展開 (テイラー展開) の展開係数として定義される。

[math] f(x) = \frac{x}{e^x-1} = \sum_{n=0}^{\infty} \frac{B_n}{n!} x^n. [/math]

この定義から、関数 x/(e^x − 1) を繰り返し微分していけばベルヌーイ数を得ることができるが、そのような手段でベルヌーイ数を得るのは容易ではない。 ベルヌーイ数を計算するには、マクローリン展開ではなく、次の漸化式を用いる。この漸化式から、ベルヌーイ数がすべて有理数であることがわかる。

[math] B_0 = 1,\quad{}B_n = -{1 \over n+1}\sum^{n-1}_{k=0}{n+1 \choose k}B_k. [/math]

ここで、[math]\left({\textstyle{n+1\atop k}}\right)[/math] は二項係数である。

上の漸化式を用いて、ベルヌーイ数を第 29 項までを算出すると下表のようになる。この表には、ベルヌーイ数が有理数であるので、分子と分母を記載している。

n	分子	分母		n	分子	分母		n	分子	分母
0	1	1		10	5	66		20	−174 611	330
1	−1	2		11	0	—		21	0	—
2	1	6		12	−691	2 730		22	854 513	138
3	0	—		13	0	—		23	0	—
4	−1	30		14	7	6		24	−236 364 091	2 730
5	0	—		15	0	—		25	0	—
6	1	42		16	−3 617	510		26	8 553 103	6
7	0	—		17	0	—		27	0	—
8	−1	30		18	43 867	798		28	−23 749 461 029	870
9	0	—		19	0	—		29	0	—

ベルヌーイ数の漸化式は、上記の関数 f(x) = x/(e^x − 1) の逆数をテイラー展開し、その 2 つの積が 1 になることから導出できる。その漸化式は厳密な計算には有用であるが、n が大きくなると途中の式の値が非常に大きくなるため、浮動小数点数を使って計算する場合、精度が著しく悪くなる計算として知られている。

奇数番目のベルヌーイ数は B₁ 以外はすべて 0 であり、偶数番目は B₀ を除いて正の数と負の数が交互に並ぶ。 ベルヌーイ数の第 3 項以降の奇数項が 0 となることは、正接関数 (tangent) のマクローリン展開から証明できる。

ベルヌーイ数の一般項

第 2 種スターリング数との関係から、次のようなベルヌーイ数の一般項を算出する公式が存在する。

[math] B_n = \sum_{j=0}^n (-1)^j\,j^n\sum_{m=j}^n \frac{1}{m+1}{m\choose j}. [/math]

この公式は、総和記号が二重になっているため、上に示した漸化式ほど手軽にベルヌーイ数を計算する公式ではない。

漸近的性質

ベルヌーイ数とリーマンゼータ関数の関係から、

[math]B_{2n} = (-1)^{n+1}\frac {2(2n)!} {(2\pi)^{2n}} \left( 1 + \frac{1}{2^{2n}} + \frac{1}{3^{2n}} + \frac{1}{4^{2n}} + \dotsb \right) [/math]

が成り立つ。従ってスターリングの公式から、n → ∞ のとき、

[math] |B_{2 n}| \sim 4 \sqrt{\pi n} \left(\frac{n}{ \pi e} \right)^{2n} [/math]

が成り立つ。

ベルヌーイ数を用いた級数展開

ベルヌーイ数は、いくつかの双曲線関数と三角関数の級数展開における展開係数となる。 ベルヌーイ数を展開係数とする関数とそのローラン級数による表現を挙げる。 まず、余接関数 (cotangent) のローラン級数展開は次のようになる。

[math] \begin{align} \coth z &= \frac1z + \sum_{k=1}^{\infty} \frac{2^{2k}B_{2k}}{(2k)!}z^{2k-1},\\ \cot z &= \frac1z + \sum_{k=1}^{\infty} (-1)^k \frac{2^{2k}B_{2k}}{(2k)!}z^{2k-1}. \end{align} [/math]

第 1 の関係式は、ベルヌーイ数が x/(e^x − 1) の展開係数であることを利用して数式変形すれば得られる。 第 2 の関係式は cot z = -i coth (-iz) であることを利用すれば、第 1 の関係式から導き出される。これらの級数の収束半径は テンプレート:Mabs < π である。 次に正接関数 (tangent) のローラン級数展開は次のようになる。

[math] \tan z = \sum_{k=1}^{\infty} \frac{(-1)^k\,(2^{2k}-4^{2k})\,B_{2k}}{(2k)!}z^{2k-1}. [/math]

この関係式は、tan z = cot z − 2cot 2z を利用して余接関数のローラン級数展開を変形すれば導出できる。 なお、この級数の収束半径は テンプレート:Mabs < π/2 である。この正接関数のローラン級数展開の展開係数による数列はタンジェント数と呼ばれる。 一方、余割関数 (cosecant) は次のようにローラン級数展開される。

[math] \csc z = \frac{1}{\sin z} = \frac{1}{z} +\sum_{k=1}^\infin \frac{(-1)^k\,(2-2^{2k})\,B_{2k}}{(2k)!}z^{2k-1}. [/math]

この関係式は、csc 2z = (tan z + cot z)/2 を利用すれば導出できる。 なお、この級数の収束半径は テンプレート:Mabs < π である。

べき乗和による導入

ベルヌーイ数は、もともと、連続する整数のべき乗和を定式化する際に、展開係数として導入された。 現代の表記法によって書くならば、定式化するべき乗和とは、

[math] S_k(n) \equiv \sum_{j=0}^{n-1} j^k = 0^k + 1^k + 2^k + \cdots + (n-1)^k [/math]

なる総和である。この総和は、ベルヌーイ数を用いて、

[math] S_k(n) = \frac{1}{k+1} \sum_{j=0}^k {k+1\choose j}B_j\,n^{k-j+1} [/math]

のように書くことができる。 ベルヌーイ数の漸化式は、べき乗和を定式化した際の考察から得られる。 さらに、ベルヌーイ数の指数型母関数が x/(e^x − 1) となることから、その母関数を現在ではベルヌーイ数の定義とする。

ヤコブ・ベルヌーイは彼の著書『推測術』でベルヌーイ数を導入した際、べき乗和を上に書いたような 0 から n − 1 にわたる和でなく、1 から n にわたる和:

[math] \hat{S}_k(n) \equiv \sum_{j=1}^{n} j^k = 1^k + 2^k + 3^k + \cdots + n^k [/math]

として扱っていた。 ベルヌーイは、その著書で整数のべき乗 n^c の和を計算する公式として、次の数式を記している。

[math] \begin{align} \int n^c & = {} \frac{1}{c+1}n^{c+1} + \frac{1}2n^c + \frac{c}2 An^{c-1} + \frac{c.c-1.c-2}{2.3.4}Bn^{c-3}\\ & \quad\quad\quad{} + \frac{c.c-1.c-2.c-3.c-4}{2.3.4.5.6}Cn^{c-5} \\ & \quad\quad\quad{} + \frac{c.c-1.c-2.c-3.c-4.c-5.c-6}{2.3.4.5.6.7.8}Dn^{c-7}+ \ldots\ldots \end{align} [/math]

この数式に記載されている展開係数 [math]A, B, C, D, \ldots[/math] がベルヌーイ数 ([math]B_2[/math] 以降) である。ベルヌーイが記した数式は、

[math] \hat{S}_k(n) = \frac{1}{k+1} \sum_{j=0}^k {k+1\choose j}\hat{B}_j\,n^{k-j+1} [/math]

に相当する。この数式に用いた展開係数 [math]\hat{B}_j[/math] は、

[math] \hat{B}_1 = 1/2\,(= -B_1),\quad \hat{B}_j = B_j\quad(j\neq1) [/math]

のように、[math]j\neq1[/math] においてベルヌーイ数と一致する。一部の文献では [math]B_j[/math] の代わりに [math]\hat{B}_j[/math] をベルヌーイ数と呼んでいる。

一方、日本ではベルヌーイとほぼ同時期に関孝和がべき乗和を定式化し、ベルヌーイ数を発見していた。 そのため、ベルヌーイ数を関・ベルヌーイ数と書いている文献もある。

一般ベルヌーイ数

一般ベルヌーイ数は代数的数で、ベルヌーイ数がリーマンゼータ函数の特殊値に関連する方法と同じ方法で、ディリクレの L-関数の特殊値に関連して定義される。

χ を mod f のディリクレ指標とすると、一般ベルヌーイ数 B_k,χ は、

[math]\sum_{a=1}^f\chi(a)\frac{te^{at}}{e^{ft}-1}=\sum_{k=0}^\infty B_{k,\chi}\frac{t^k}{k!}[/math]

により定義される。B_1,1 = 1/2 を除き、任意のディリクレ指標 χ に対し、χ(−1) ≠ (−1)^k であれば、B_k,χ = 0 である。

正でない整数におけるリーマンゼータ関数の値とベルヌーイ数の間の関係を一般化し、全ての整数 k ≥ 1 に対し、

[math]L(1-k,\chi)=-\frac{B_{k,\chi}}{k}[/math]

が成り立つ。ここに L(s, χ) は χ のディリクレの L-関数である。

関連項目

• ヤコブ・ベルヌーイ
• 関孝和
• ファウルハーバーの公式
• ベルヌーイ多項式
• ゼータ関数
• 三角関数