この記事は多重ガンマ関数の入門的な記事です。ゼータ関数や解析接続のお気持ち程度の理解があると役に立つかもしれません。今回は初回ということで多重ガンマの導入まではいかずに、Hurwitzのゼータ関数とレルヒの公式までを扱います。
内容には誤りがある可能性があります。留意してください。
素朴なHurwitzのゼータ関数は、$Re(s)>1$の複素数に対して次の式で定義されます。
$$ ζ(s,x):= \sum_{n=0}^{∞} (n+x)^{-s} $$
$x=1$とするとリーマンゼータ関数となるので、この関数はリーマンゼータ関数の一般化になっているわけですね。
この関数は$Re(s)>1$に対して次の積分表示を持ちます。
$$ ζ(s,x)= \frac{1}{Γ(s)} \int_{0}^{\infty} \frac{e^{-xt}}{1-e^{-t}} t^{s-1}dt $$
ガンマ関数については
この記事
かWikipediaなどを参照してください。
ちなみに、この積分表示から被積分関数の$t\rightarrow 0$での漸近的な振る舞いが$O(t^{s-2}$)であることがわかるので、ディリクレ級数表示と積分表示それぞれの収束範囲が先ほど言った$Re(s)>1$であることがわかります。証明は簡単なので済ませてしまいましょう。
$$
ζ(s,x)Γ(s)= \sum_{n=0}^{\infty} (n+x)^{-s}\cdot \int_{0}^{\infty} e^{-t}t^{s-1}dt
$$
$t=(n+x)u$と置換して、
$$
=\sum_{n=0}^{\infty} (n+x)^{-s}\cdot \int_{0}^{\infty} e^{-(n+x)u}(n+x)^{s-1}u^{s-1}du \cdot (n+x)
$$
$$
=\sum_{n=0}^{\infty} \int_{0}^{\infty} e^{-(n+x)t}t^{s-1}dt
$$
見栄えのため積分変数をtに戻しました。
適当な収束性の元無限和と積分を交換できるので、
$$
= \int_{0}^{\infty} e^{-xt} \sum_{n=0}^{\infty} e^{-nt}t^{s-1}dt
$$
$$
= \int_{0}^{\infty} \frac{e^{-xt}}{1-e^{-t}} t^{s-1}dt
$$
ただしここで等比級数の和公式を用いた。(Q.E.D)
ちなみに上記のような変換は一般のディリクレ級数にも可能で、メリン変換と呼ばれています。
$$
ζ_A(s):= \sum_{a_n\in A}a_n^{-s}
$$
$$
\psi(t):= \sum_{a_n\in A}e^{-a_nt}
$$
に対して、
$$
ζ_A(s)= \frac{1}{Γ(s)} \int_{0}^{\infty} \psi(t)t^{s-1}dt
$$
証明は先ほどと同じ流れです。
先ほどのHurwitzのゼータ関数は$Re(s)>1$でしか定義されていませんでしたが、実は$s=1$を除く全ての複素数に解析接続できます。早速やっていきましょう。
$ \psi(t,x)= \frac{e^{-xt}}{1-e^{-t}} $とします。
まず積分表示を3つに分解します。
$$
ζ(s,x)= \frac{A(s)}{Γ(s)} + \frac{B(s)}{Γ(s)} + \frac{C(s)}{Γ(s)}
$$
それぞれ
$$
A(s)= \int_{1}^{\infty} \psi(t,x)t^{s-1}dt
$$
$$
B(s)= \int_{0}^{1} (\sum_{n=-1}^{m} a_n(x)t^n)t^{s-1}dt
$$
$$
C(s)= \int_{0}^{1} \{\psi(t,x)-\sum_{n=-1}^{m} a_n(x)t^n\}t^{s-1}dt
$$
ここで$a_n(x)$というのは$\psi(t,x)$の$t=0$でのローラン展開の係数で、
$$
\psi(t,x)= \frac{e^{-xt}}{1-e^{-t}} = \sum_{n=-1}^{\infty} a_n(x)t^n
$$
と定義されます。
要するにローラン展開の第m項までをくくり出しているわけです。
さて、先ほど分解したそれぞれの関数の収束範囲を見てみましょう。
まず$A(s)$は積分範囲内に極はなく、$t\rightarrow \infty$で$e^{-xt}$が支配的なので急減少し全ての複素数$s$に対して収束します。
$B(s)$は簡単な積分なので計算すると
$$
B(s)= \sum_{n=-1}^{m} \frac{a_n(x)}{n+s}
$$
となるので、これは$s=1,0,-1,...,-m$に一位の極を持つだけの全平面で有理型関数になりましたね。
$C(s)$は少し面倒ですが、$t\rightarrow 0$の評価を考えます。
ローラン展開の定義より$ \psi(t,x)- \sum_{n=-1}^{m} a_n(x)t^n=O(t^{m+1}) $なので、
被積分関数の$t\rightarrow 0$でのオーダーは$O(t^{m+s})$となります。よって積分は$s>-m-1$で収束することがわかります。ここで$m$は任意の自然数なので、Hurwitzのゼータ関数は複素平面全域に解析接続できました。
せっかくなので特殊値を計算してみましょう。
公式
0を含む自然数$n$に対し、
$$
ζ(-n,x)=(-1)^nn!a_n(x)
$$
$$
ζ(-n,x)=\lim_{s \rightarrow -n} \frac{A(s)}{Γ(s)} + \frac{B(s)}{Γ(s)} + \frac{C(s)}{Γ(s)}
$$
ここで、ガンマ関数は0を含む負の整数点で一位の極を持つので、極を持たない$A(s)$と$C(s)$は0になります。
よって$B(s)/Γ(s)$を考えればいいですが、$B(s)$も先ほどと同じ理由で$\frac{a_n(x)}{s+n}$以外の項は全て0になります。よって、
$$
\lim_{s\rightarrow -n}a_n(x)\cdot \frac{1}{Γ(s)(s+n)}
$$
$$
=\lim_{s\rightarrow -n}a_n(x)\cdot \frac{s(s+1)...(s+n-1)}{(s+n)Γ(s+n)}
$$
$$
= a_n(x)\cdot\frac{(-n)\cdot(-n+1)...(-1)}{1}
$$
$$
=(-1)^nn!a_n(x)
$$
ただし、ガンマ関数の関数等式を用いた。
特に、$ a_0(x)= \frac{1}{2} -x $,$ a_1(x)= \frac{1}{12} - \frac{x}{2} + \frac{x^2}{2} $(これは$\psi(t,x)$をベルヌーイ数を使って頑張って計算できます)なので、次の式が成り立ちます。
$$
ζ(0,x)= \frac{1}{2} -x
$$
$$
ζ(-1,x)=- \frac{1}{12} + \frac{x}{2} - \frac{x^2}{2}
$$
これらの式を使って自然数の正規化和が$ー\frac1{12}$というような話もできるのですが、本題からズレるので今回は扱いません。
Hurwitzゼータ関数が用意できたことで、やっと準備が整いました。レルヒの公式とは次のゼータ関数とガンマ関数を対応付ける、多重ガンマにも関わるとても重要な公式です。
$$ \left. \frac{\partial}{\partial s} ζ(s,x) \right|_{s=0}=log \frac{Γ(x)}{\sqrt{2\pi}} $$
解析接続されたHurwitzゼータ関数の$s=0$での導関数に驚くことにガンマ関数が現れました。この式は次のように書くこともできます。
$$ \Rprod{n=0}{\infty}{(n+x)=\frac{\sqrt{2\pi}}{Γ(x)}} $$
この式の左辺は正規化積といって、普通の積ではありません(下が繋がっていて総乗記号と違うことに注意)
正規化積とは$s=0$で正則に解析接続できる全てのディリクレ級数に対して定義できるもので、次のように定義されます。
$ ζ_A(s):= \sum_{a_n\in A}a_n^{-s} $ とする。
この時、
$ \Rprod{a_n\in A}{}{a_n:=\left.exp( -\frac{\partial}{\partial s}ζ_A(s)\right|_{s=0} )} $
なぜこのように定義されるかというと、ディリクレ級数表示を形式的に微分して0を代入したものに-1をかけてexpすると、形式的にa_nに対する積が現れますが(これは簡単に確認できる)、ディリクレ級数表示がs=0で成り立たない場合はその式はもちろん発散します。しかしs=0まで解析接続できている場合にはその値にも意味が生まれると期待できそうなので、なら新しい記号を定義して扱ってやろうじゃないか、みたいなモチベーションです。
この正規化積、積と呼ばれるだけあっていくつかのそれっぽい性質を備えています。
&&&fml 正規化積の性質
I,Jを集合Aの分割とします。
このとき、
これらの性質はそれぞれゼータ関数を二つに分割したり定数をかけて微分してあげるだけで簡単に示せます。まるで本当の積のようですね。特に1の性質を使うと正規化積の中から有限項を抜き出して普通の積として書くことができます。
ちなみに、当たり前ですが例えば$ζ_A(s)=\sum_{n=1}^5n^{-s}$のようにs=0でも収束するような関数にすると正規化積は普通の積と一致します。
話が逸れてしまいましたが、レルヒの公式
$$
\Rprod{n=0}{\infty}{(n+x)=\frac{\sqrt{2\pi}}{Γ(x)}}
$$
にx=1を代入すると、
$$
\Rprod{n=1}{\infty}{n=\sqrt{2\pi}}
$$
という結果が得られます。
自然数の無限積はもちろん発散してしまいますが、正規化積を用いると一定の意味付けが行えるということです。
ところで、逆にガンマ関数が正規化積を用いて表されていると見ることもできて、正規化積を見るとガンマ関数の逆数の零点が一目瞭然です。
また、レルヒの公式の正規化積でxが自然数だと考えると、左辺は形式的にxから$\infty$までの無限積だと思うことができて、右辺は形式的に$1$から$\infty$までの積$\sqrt{2\pi}$を$1$から$x-1$までの積$Γ(x)$で割っていると考えられるので辻褄が合います。面白いですね。
話がめちゃくちゃ長引いてしまいましたが、早速レルヒの公式の証明を参考文献にある現代三角関数論に沿って三段階に分けて行います。
$f(x)=\left. \frac{\partial}{\partial s} ζ(s,x) \right|_{s=0}-logΓ(x)$として、
(1)$f''(x)=0$を示す。これによって$f(x)=ax+b$がわかる。
(2)$f(x+1)=f(x)$を示す。これによって$f(x)$=bがわかる。
(3)$f(\frac12)$=$-\frac12log2\pi$を示す。
まず先に3つの補題を証明します。
&&&lem xでの偏微分
$$
\frac{\partial}{\partial x}ζ(s,x)=-sζ(s+1,x)
$$
$$ \frac{d^2}{dx^2} logΓ(x)=ζ(2,x) $$
$$ ζ(s,x+1)=ζ(s,x)-x^{-s} $$
$$
\frac{\partial}{\partial x}\sum_{n=0}^{\infty}(n+x)^{-s}=-s \sum_{n=0}^{\infty} (n+x)^{-s-1}
$$
$$
=-sζ(s+1,x)
$$
$$
\frac{1}{Γ(x)} =xe^{γx} \prod_{n=1}^{\infty} \left(1+ \frac{x}{n}\right)e^{- \frac{x}{n} }
$$
なので、(詳しくは参考文献の前に書いた記事をみてください)
$$
-logΓ(x)=logx+γx+ \left\{\sum_{n=1}^{\infty} log( 1+ \frac{x}{n}) - \frac{x}{n} \right\}
$$
である。これをxで二階微分すると、
$$
\frac{d}{dx} -logΓ(x)= \frac{1}{x} +γ+ \sum_{n=1}^{\infty} \left( \frac{1}{n+x}- \frac{1}{n} \right)
$$
$$
\frac{d^2}{dx^2} -logΓ(x)=- \frac{1}{x^2} - \sum_{n=1}^{\infty} \frac{1}{(n+x)^2}
$$
$$
=- \sum_{n=0}^{\infty} (n+x)^{-2}
$$
$$
=-ζ(2,x)
$$
$$
ζ(s,x+1)= \sum_{n=0}^{\infty} (n+x+1)^{-s}
$$
$$
= \sum_{n=1}^{\infty} (n+x)^{-s}=ζ(s,x)-x^{-s}
$$
(1)の証明を始めます。
$$
f''(x)= \left.\frac{\partial ^3}{\partial x^2\partial s} ζ(s,x)\right|_{s=0}- \frac{d^2}{dx^2} logΓ(x)
$$
補題2と3より
$$
= \left. \frac{\partial }{\partial s}\bigr(s(s+1)ζ(s+2,x)\bigr)\right|_{s=0}-ζ(2,x)
$$
$$
=ζ(2,x)-ζ(2,x)=0
$$
よって(1)が示された。
(2)の証明
$$
f(x+1)=\left. \frac{\partial}{\partial s} ζ(s,x+1) \right|_{s=0}-logΓ(x+1)
$$
補題3とガンマ関数の周期性より
$$
=\left. \frac{\partial}{\partial s} \bigl(ζ(s,x)-x^{-s}\bigr) \right|_{s=0}-logΓ(x)-logx
$$
$$
=\left. \frac{\partial}{\partial s} ζ(s,x)\right|_{s=0}+logx\cdot x^{-0}-logΓ(x)-logx
$$
$$
=f(x)
$$
(3)の証明
$$
f(x)=\left. \frac{\partial}{\partial s} ζ(s, \frac{1}{2} ) \right|_{s=0}-logΓ( \frac{1}{2} )
$$
まずゼータ関数の方から計算します。
$$
ζ(s, \frac{1}{2} )= \sum_{n=0}^{\infty} \bigr(n+ \frac{1}{2} \bigr)^{-s}=2^s \sum_{n=0}^{\infty} (2n+1)^{-s}
$$
$$
=2^s\Bigr\{ \sum_{n=1}^{\infty} n^{-s}- \sum_{n=1}^{\infty} (2n)^{-s}\Bigr\}
$$
$$
=2^s(1-2^{-s})ζ(s)=(2^s-1)ζ(s)
$$
これを微分すると、
$$
\left. \frac{\partial}{\partial s} ζ(s, \frac{1}{2}) \right|_{s=0}=(1-1)ζ'(0)+log2ζ(0)
$$
$$
= -\frac{1}{2} log2
$$
ただし、先に計算したHurwitzゼータ関数の特殊値にs=0,x=1を代入したものを用いた。
また、$ Γ( \frac{1}{2} )=\sqrt\pi $なので、(証明は略しますが調べればいくつも出てきます)
$$ f( \frac{1}{2})=- \frac{1}{2}log2- \frac{1}{2} log\pi
$$
$$
=- \frac{1}{2} log2\pi
$$
以上の結果から、
$$
\left. \frac{\partial}{\partial s} ζ(s,x) \right|_{s=0}-logΓ(x)=- \frac{1}{2} log2\pi
$$
$$
\left. \frac{\partial}{\partial s} ζ(s,x) \right|_{s=0}=log \frac{Γ(x)}{\sqrt{2\pi}}
$$
を得るので、レルヒの公式を証明することができました。
モチベーション次第ですが、次回にこのレルヒの公式からガンマ関数を拡張した多重ガンマ関数、さらにそれを拡張したBarnes-Milnor Type(BM型)多重ガンマ関数の解説を行いたいと思います。では。