この記事ではオイラーの級数変換公式というものとそれによるゼータ関数の解析接続について解説していきます。
オイラーの級数変換とは以下のような変換公式のことを言います。
数列$a_n$に対し
$$b_n=\sum^n_{k=0}(-1)^j\binom nka_k$$
とおくと
$$\sum^\infty_{n=0}(-1)^na_n=\sum^\infty_{n=0}\frac{b_n}{2^{n+1}}$$
が成り立つ(収束性については後述)。
数列$a_n$に対し差分作用素$\D$およびその反復作用を
$$(\D a)_n=a_n-a_{n+1},\quad\D^ka=\D(\D^{k-1}a)$$
によって定めます。つまり
$$(\D^ka)_n=\sum^k_{j=0}(-1)^j\binom kja_{n+j}$$
が成り立つことに注意しましょう。
左辺が収束するとき
$$\sum^\infty_{n=0}(-1)^na_n
=\sum^{N-1}_{k=0}\frac{(\D^ka)_0}{2^{k+1}}
+\sum^\infty_{n=0}(-1)^n\frac{(\D^Na)_n}{2^N}$$
が成り立つ。
\begin{align}
2\sum^\infty_{n=0}(-1)^na_n
&=\l(a_0+\sum^\infty_{n=0}(-1)^{n+1}a_{n+1}\r)+\sum^\infty_{n=0}(-1)^na_n\\
&=a_0+\sum^\infty_{n=0}(-1)^n(a_n-a_{n+1})\\
&=(\D^0a)_0+\sum^\infty_{n=0}(-1)^n(\D a)_n
\end{align}
が成り立つので
$$\sum^\infty_{n=0}(-1)^n\frac{(\D^Na)_n}{2^N}=\frac2{2^{N+1}}\sum^\infty_{n=0}(-1)^n(\D^Na)_n$$
にも同様の変形をしていくことでわかる。
あとは"適当な条件下"で
$$\lim_{N\to\infty}\sum^\infty_{n=0}(-1)^n\frac{(\D^Na)_n}{2^N}=0$$
となることが示せるので主張が得られます。
上の方法では
$$\lim_{N\to\infty}\sum^\infty_{n=0}(-1)^n\frac{(\D^Na)_n}{2^N}=0$$
を示すのに少し面倒な手順が必要ですが、二重級数に対するフビニの定理を使えばかなり近道ができます。
二重数列$a_{m,n}$に対し常に
$$\sum^\infty_{m=0}\sum^\infty_{n=0}|a_{m,n}|
=\sum^\infty_{n=0}\sum^\infty_{m=0}|a_{m,n}|$$
が成り立ち、またこれが有限値を取る(つまり$\infty$ではない)ときは
$$\sum^\infty_{m=0}\sum^\infty_{n=0}a_{m,n}
=\sum^\infty_{n=0}\sum^\infty_{m=0}a_{m,n}$$
も成り立つ。
ここでこれは二重級数の順序交換といったワードとしてよく聞くトピックではありますが、無限和はルベーグ積分とみなせることから個人的な好みとしてこれをフビニの定理と言うことにしています。二重級数の順序交換やルベーグ積分のフビニの定理の証明については随所にあると思うのでここでは特に解説しません。
いまフビニの定理を用いるとオイラーの変換公式は簡単に示せます。
左辺が絶対収束するとき
$$\sum^\infty_{n=0}(-1)^na_n=\sum^\infty_{n=0}\frac{(\D^na)_0}{2^{n+1}}$$
が成り立ち、この右辺も絶対収束する。
$$\sum^\infty_{n=0}x^n=\frac1{1-x}$$
を$m$回微分することで
$$\sum^\infty_{n=m}\binom nmx^{n+1}=\l(\frac{x}{1-x}\r)^{m+1}$$
となるので$x=1/2$とすることで
$$\sum^\infty_{n=m}\binom nm\frac1{2^{n+1}}=1$$
が成り立つことに注意する。
いま仮定より
\begin{align}
\sum^\infty_{n=0}\l|\frac{(\D^na)_0}{2^{n+1}}\r|
&\leq\sum^\infty_{n=0}\frac1{2^{n+1}}\sum^n_{m=0}\binom nm|a_m|\\
&=\sum^\infty_{m=0}|a_m|\sum^\infty_{n=m}\binom nm\frac1{2^{n+1}}\\
&=\sum^\infty_{m=0}|a_m|<\infty
\end{align}
と諸々の収束性がわかるのでフビニの定理より
\begin{align}
\sum^\infty_{n=0}\frac{(\D^na)_0}{2^{n+1}}
&=\sum^\infty_{n=0}\frac1{2^{n+1}}\sum^n_{m=0}(-1)^m\binom nma_m\\
&=\sum^\infty_{m=0}(-1)^ma_m\sum^\infty_{n=m}\binom nm\frac1{2^{n+1}}\\
&=\sum^\infty_{m=0}(-1)^ma_m
\end{align}
を得る。
オイラーの変換公式はしばしば級数の収束を加速させるのに便利な公式の一つとなっています。
例えば$a_n$が指数関数的、つまりある$|\a|<1$を用いて$a_n\sim \a^n$とみなせるとき
$$(\D^ka)_n\sim (1-\a)^ka_n$$
のように評価できるので
$$\sum^\infty_{n=0}(-1)^na_n=\sum^\infty_{n=0}\frac{(\D^na)_0}{2^{n+1}}$$
の左辺の収束速度は$|\a|^n$程度であるのに対し、右辺の収束速度は$(|1-\a|/2)^n$程度となります。言うほど加速してないな。場合によっては遅くなりますね。
また例えば$a_n=f(n)$みたいな場合は
$$f(n+1)-f(n)=\int^1_0f'(n+x)dx$$
のように変形することで
$$(\D^ka)_n=(-1)^k\int^1_0\int^1_0\cdots\int^1_0f^{(k)}(n+x_1+x_2+\cdots+x_k)dx_1dx_2\cdots dx_k$$
と表せるのでこれを上手いこと評価できれば何かハッピーな感じがしませんか?
実際ゼータ関数のオイラー変換を考えることによって級数の収束範囲を広げることができ、それは$\Re(s)\leq1$への解析接続を与えることになります。
$\Re(s)>1$において
$$\z(s)=\frac1{1-2^{1-s}}\sum^\infty_{n=0}\frac{(-1)^n}{(n+1)^s}$$
が成り立つ。
ちなみにこの右辺の収束範囲からこれは$\Re(s)>0$への解析接続を与えています。
\begin{align}
(1-2^{1-s})\sum^\infty_{n=1}\frac1{n^s}
&=\sum^\infty_{n=1}\frac1{n^s}-2\sum^\infty_{n=1}\frac1{(2n)^s}\\
&=\sum^\infty_{n=1}\frac1{(2n+1)^s}-\sum^\infty_{n=1}\frac1{(2n)^s}\\
&=\sum^\infty_{n=0}\frac{(-1)^n}{(n+1)^s}
\end{align}
とわかる。
これをオイラー変換することで以下が得られます(この右辺の絶対収束範囲は$\Re(s)>1$であることに注意する)。
$\Re(s)>1$において
$$\z(s)=\frac1{1-2^{1-s}}\sum^\infty_{n=0}\frac{\D^n1^{-s}}{2^{n+1}}$$
が成り立つ。ただし
$$\D^k1^{-s}=(\D^k(n+1)^{-s})_0=\sum^k_{j=0}(-1)^j\binom kj(j+1)^{-s}$$
とした。
後はこの右辺の級数が$\C$全体で正則関数を定めることを確かめればこれがゼータ関数の解析接続を与えることになります。
ところでオイラー変換を途中でやめると補題2のように以下が成り立つのでした。
$\Re(s)>1$において
$$(1-2^{1-s})\z(s)=\sum^{N-1}_{k=0}\frac{\D^k1^{-s}}{2^{k+1}}
+\sum^\infty_{n=0}(-1)^n\frac{\D^N(n+1)^{-s}}{2^N}$$
が成り立つ。
実はこの右辺の級数は$\Re(s)>1-N$で絶対一様収束することが示せ、また補題2の証明からもわかるように
$$\sum^\infty_{n=0}(-1)^n\frac{\D^N(n+1)^{-s}}{2^N}
=\sum^\infty_{n=N}\frac{\D^n1^{-s}}{2^{n+1}}$$
はオイラー変換となっているのでこの右辺も$\Re(s)>1-N$で絶対一様収束する、つまり正則関数を定めることがわかります。
上での議論から
$$\D^kf(n)
=(-1)^k\int^1_0\int^1_0\cdots\int^1_0f^{(k)}(n+x_1+x_2+\cdots+x_k)dx_1dx_2\cdots dx_k$$
と表せたので
$$\D^Nn^{-s}
=(s)_N\int^1_0\int^1_0\cdots\int^1_0(n+x_1+x_2+\cdots+x_k)^{-s-N}dx_1dx_2\cdots dx_k$$
が成り立つ。ただし$(s)_N$はポッホハマー記号
$$(s)_N=s(s+1)\cdots(s+N-1)$$
とした。
そしてこれは$\Re(s)+N>0$において
$$|\D^Nn^{-s}|<\frac{|(s)_N|}{n^{\Re(s)+N}}$$
と評価できるので、$\Re(s)+N>1$において
$$\sum^\infty_{n=0}\l|\frac{\D^N(n+1)^{-s}}{2^N}\r|<\frac{|(s)_N|}{2^N}\z(\Re(s)+N)$$
と(広義)一様収束することがわかる。
以上により以下の公式が得られます。
$\C\setminus\{1\}$において
$$\z(s)
=\frac1{1-2^{1-s}}\sum^\infty_{n=0}\frac1{2^{n+1}}\sum^n_{m=0}(-1)^m\binom nm(m+1)^{-s}$$
が成り立つ。
上の公式はゼータ関数の解析的な性質、例えば関数等式を導出するなどには向いていませんが、特定の点における値を求めるには何かと役に立ちます。
例えば今回参考にした
INTEGERSの記事
ではWorpitzkyの公式
$$B_{m+1}=\frac{m+1}{2^{m+1}-1}\sum^\infty_{n=0}\frac1{2^{n+1}}\sum^n_{j=0}(-1)^j\binom nj(j+1)^m$$
から負の整数点における特殊値
$$\z(-m)=-\frac{B_{m+1}}{m+1}$$
を導出していたり、
tsujimotterさんの記事
ではゼータ関数を数値計算し、そのグラフを描くのに利用しています。
ここではそういった例の一つとして$\z(0),\z'(0)$の値を求めてみようと思います。今回使うのはオイラー変換された公式ではなく、補題6として示された以下の公式($N=1$)
$$(1-2^{1-s})\z(s)=\frac12+\frac12\sum^\infty_{n=1}(-1)^{n-1}\l(\frac1{n^s}-\frac1{(n+1)^s}\r)$$
を利用します。この右辺は$\Re(s)>-1$において一様収束します(絶対収束とは限らない)。
$$(1-2^{1-s})\z(s)=\sum^{N-1}_{k=0}\frac{\D^k1^{-s}}{2^{k+1}}
+\sum^\infty_{n=0}(-1)^n\frac{\D^N(n+1)^{-s}}{2^N}$$
の右辺は$\Re(s)>-N$において一様収束する。
\begin{align}
\sum^\infty_{n=0}(-1)^n\D^N(n+1)^{-s}
&=\sum^\infty_{n=0}(\D^N(2n+1)^{-s}-\D^N(2n+2)^{-s})\\
&=\sum^\infty_{n=0}\D^{N+1}(2n+1)^{-s}
\end{align}
と変形すると、この右辺は
\begin{align}
\sum^\infty_{n=0}|\D^{N+1}(2n+1)^{-s}|
&\leq\sum^\infty_{n=0}\frac{|(s)_{N+1}|}{(2n+1)^{\Re(s)+N+1}}\\
&\leq|(s)_{N+1}|\z(\Re(s)+N+1)
\end{align}
と評価できるので$\Re(s)>-N$において一様収束することがわかる。
$\eta(s)=(1-2^{1-s})\z(s)$とおくと
$$\eta'(0)=\frac12\log\frac\pi2$$
が成り立つ。
$\Re(s)>-1$において
$$\eta(s)=\frac12+\frac12\sum^\infty_{n=1}(-1)^{n-1}\l(\frac1{n^s}-\frac1{(n+1)^s}\r)$$
が成り立つので、これを微分すると
$$\eta'(s)
=-\frac12\sum^\infty_{n=1}(-1)^{n-1}\l(\frac{\log n}{n^s}-\frac{\log(n+1)}{(n+1)^s}\r)$$
つまり
\begin{align}
\eta'(0)
&=-\frac12\sum^\infty_{n=1}(-1)^{n-1}(\log n-\log(n+1))\\
&=-\frac12\sum^\infty_{n=1}(\log(2n-1)-2\log2n+\log(2n+1))\\
&=-\frac12\log\l(\prod^\infty_{n=1}\frac{(2n-1)(2n+1)}{(2n)^2}\r)\\
&=-\frac12\log\l(\prod^\infty_{n=1}\l(1-\frac1{(2n)^2}\r)\r)
\end{align}
と表せる。
ここでウォリスの公式
$$\prod^\infty_{n=1}\l(1-\frac1{(2n)^2}\r)=\frac{\sin\frac\pi2}{\frac\pi2}=\frac2\pi$$
に注意すると
$$\eta'(0)=\frac12\log\frac\pi2$$
を得る。
$$\z(0)=-\frac12,\quad
\z'(0)=-\frac12\log2\pi$$
が成り立つ。
$\z(0)$については
$$(1-2^{1-s})\z(s)=\frac12+\frac12\sum^\infty_{n=1}(-1)^{n-1}\l(\frac1{n^s}-\frac1{(n+1)^s}\r)$$
において$s=0$とすることで$-\z(0)=1/2$が得られる。
また
\begin{align}
\eta'(s)&=(2^{1-s}\log2)\z(s)+(1-2^{1-s})\z'(s)\\
\eta'(0)&=(2\log2)\z(0)-\z'(0)
\end{align}
に注意すると
\begin{align}
\z'(0)
&=(2\log2)\z(0)-\eta'(0)\\
&=-\log2-\frac12\log\frac\pi2\\
&=-\frac12\log2\pi
\end{align}
を得る。