たけのこ赤軍の自由帳

反復積分とNona ReevesとPrinceとRED SPIDER

I Wanna Be Your Modular

この記事はこの曲を聴きながら読むのがオススメです:
Nona Reeves - Sweet Survivor

NONA REEVES 『Sweet Survivor』




本記事は 日曜数学 Advent Calendar 2017 - Adventar 23日目の記事です。22日目の記事はmattyuu氏による p進距離はなぜ特別か? - mattyuuの数学ネタ集 でした。



世の中には「モジュラー形式」というものがあります。

よく知られているものだと、フェルマーの最終定理を証明する際にワイルズがこれと楕円曲線との関連を示したというところです。

これを黒川信重先生による"絶対数学"の舞台に乗せてできたものが Absolute Modular Form、「絶対保型形式」です。定義は少々面倒なので、今日はすこし具体例を眺めてみましょう。

1900年ごろ、バーンズはガンマ関数の多重化というアイデアを考えつきました―――これこそが黒川先生らの考える多重三角関数論の源流です。

これはレルヒの公式の一般化によって定義されるものであり、通常のガンマ関数と同じように無限積で書き下すことができます。

\begin{eqnarray}\displaystyle\Gamma(w)^{-1}&=&\exp\left(\left.\frac{\partial }{\partial s}\zeta(s,w)\right|_{s=0}\right)\\&=&\rho_1(1)we^{\gamma w}\prod_{0\neq{n}\geq{0}}^{} \left(1+\frac{w}{n}\right)e^{-\frac{w}{n}}\\\Gamma_r(w,{\bf \omega})&=&\exp\left(\left.\frac{\partial}{\partial s}\zeta_r(s,w,\omega)\right|_{s=0}\right)\\&=&\rho_r({\bf \omega})we^{P_{\omega}(w)}\prod_{{\bf 0}\neq{\bf n}\geq{\bf 0}}^{}\left(1+\frac{w}{{\bf n}\cdot{\bf \omega}}\right)\exp\left(\sum_{k=1}^{r} \left(-\frac{w}{\bf n\cdot{\bf \omega}}\right)^k\right) \end{eqnarray}

ここで \begin{eqnarray}\displaystyle\zeta_r(s,w,\omega)=\sum_{{\bf n}\geq{\bf 0}}^{} ({\bf n}\cdot\omega+w)^{-s},\end{eqnarray}

P_{\omega}(w)\omega によってきまる wr多項式で、 {\bf 0}\neq{\bf n}\geq{\bf 0} というのは {\bf n}=(n_1,\cdots n_r)\in\mathbb{Z}^r に対して n_i\geq{0}(i=1,2,3,\cdots) かつ {\bf n}\neq{(0,\cdots 0)} であることをいいます。

\omega はボールド体ではない(TeXで打てなかった)ですが、複素数 r 個の組 (\omega_1,\cdots,\omega_r) だと理解してください。

{\bf n}\cdot{\omega}=n_1\omega_1+n_2\omega_2+\cdots+n_r\omega_r です。




さてこれらを比較してみると、まず \rho_1(1) が気になりますね。 \Gamma_r はもともと \Gamma を一般化したものですから、 \rho_1(1) を一般化して \rho_r(\omega) が出てきたのだと推測できます。

o-v-e-r-h-e-a-t.hatenablog.com

で紹介したレルヒの公式からすると、 \rho_1(1)=\sqrt{2\pi} でなくてはならないはずです。そして、この値はまた

\begin{eqnarray}\displaystyle\sqrt{2\pi}=\prod_{0\neq{n}\geq{0}}^{} n\cdot{1}\end{eqnarray}

でもありました。複素数 r 個の組 \omega は「周期」と呼ばれ、多重ガンマ関数 \Gamma_r の補助的な変数になっています(メインの変数はもちろん w )。つまり、一般化を逆にたどると \Gamma_1(w,1)=\Gamma(w) であるということですね。

ここまでくると、 \rho_r(\omega) をどう定義すればいいかはもう見当がつきます。もちろん

\begin{eqnarray}\displaystyle\rho_r(\omega):=\prod_{{\bf 0}\neq{\bf n}\geq{\bf 0}}^{} {\bf n}\cdot\omega\end{eqnarray}

ということですね。ところでこれは正規積なので、書き直してみると

\begin{eqnarray}\displaystyle\zeta_r^{\star} (s,\omega)=\sum_{{\bf 0}\neq{\bf n}\geq{\bf 0}}^{} ({\bf n}\cdot\omega)^{-s}\end{eqnarray}

として

\begin{eqnarray}\displaystyle\rho_r(\omega):=\exp\left(\left.\frac{\partial}{\partial s}\zeta_r^{\star} (s,\omega)\right|_{s=0}\right)\end{eqnarray}

とも書けますね。ということで、こいつが今回扱う「絶対保型形式の一例」です。本記事の目標は r=2 における特殊な例、 \rho_2(1,i) を明示的に求めることとしましょう。




それを示すためにはたくさんの道具が必要なのですが、そのうちの一つをここで定義します:

\begin{eqnarray}\displaystyle\zeta_4^M(s):=\sum_{m\geq{1},n\geq{0}}^{} (m+ni)^{-s}\end{eqnarray}

これは私が導入した関数で、私が勝手に4次元モノイドゼータと呼んでいるものです。詳細については2018/1/7に数学カフェで行われる講演の資料(終了後公開予定)を参照してください。


[補題1(ポアソン和公式)]

\begin{eqnarray}\displaystyle\sum_{n\in\mathbb{Z}}^{} f(n)=\sum_{m\in\mathbb{Z}}^{} \widehat{f}(m)\end{eqnarray}

[証明]

 f:\mathbb{R}\rightarrow\mathbb{R} に対してのフーリエ変換の定義

\begin{eqnarray}\displaystyle\widehat{f}(y):=\int_{-\infty}^{\infty} f(x)e^{-2\pi ixy} dx\end{eqnarray}

を思い出しておきます。

\displaystyle h(x)=\sum_{n\in{\mathbb{Z}}}^{} f(x+n)

としたとき、 h(x) は明らかに周期 1 なのでフーリエ展開

h(x)=\sum_{n\in{\mathbb{Z}}}^{} a_ne^{2\pi inx}

をもつ、というのがわかります。

\begin{eqnarray}a_n&=&\int_0^1 h(x)e^{-2\pi inx}dx\\&=&\sum_{m\in{\mathbb{Z}}}^{} \int_0^1 f(x+m)e^{-2\pi inx}dx\end{eqnarray}

ですが、 y=x+m で変数変換すると積分範囲は [0,1] から [m,m+1] となり、 \displaystyle\frac{dy}{dx}=1 より

\begin{eqnarray}a_n&=&\sum_{m\in\mathbb{Z}}\int_m^{m+1} f(y)e^{-2\pi in(y-m)}dy\\&=&\int_{-\infty}^{\infty} f(y)e^{-2\pi iny}dy\\&=&\widehat{f}(n).\end{eqnarray}

ゆえに

\begin{eqnarray}h(x)&=&\sum_{n\in\mathbb{Z}}^{} f(x+n)\\&=&\sum_{n\in\mathbb{Z}}^{} \widehat{f}(n)e^{2\pi inx}.\end{eqnarray}

となります。両辺で x=0 とすると目的の公式を得られます。[証明終]




[補題2(リプシッツ和公式)]

\begin{eqnarray}\displaystyle\sum_{m\in\mathbb{Z}}(m+z)^{-s}=\frac{(-2\pi i)^s}{\Gamma(s)}\sum_{m=1}^{\infty} m^{s-1}e^{2\pi imz}\end{eqnarray}

[証明] ※日本語で証明している本・文献が見つからなかったのでクノップとロビンスの論文"Easy proofs of Riemann's functional equation for \zeta(s) and of Lipschitz summation"のTheorem 1の証明を参考にしています。

f(x)=\begin{cases}
    x^{s-1}e^{2\pi ixz} & (x>0) \\
    0 & (x\leq{0})
  \end{cases}

として、ポアソン和公式を適用します:

\begin{eqnarray}\displaystyle\sum_{m=1}^{\infty} m^{s-1}e^{2\pi imz}&=&\sum_{n\in\mathbb{Z}}^{} f(n)\\&=&\sum_{m\in\mathbb{Z}}^{} \widehat{f}(-m)\\&=&\sum_{m\in\mathbb{Z}}^{} \int_{0}^{\infty} x^{s-1}e^{2\pi ixz}e^{2\pi imx}dx\\&=&\sum_{m\in\mathbb{Z}}^{} \int_{0}^{\infty} x^{s-1}e^{-(-2\pi i(m+z))x}dx\\&=&\sum_{m\in\mathbb{Z}}^{} (-2\pi i(m+z))^{-s}\Gamma(s)\end{eqnarray}

となって、両辺を整理すると証明できます。[証明終]


さてこれを使うことで、 \zeta_4^M の解析接続(のような何か)を行うことが可能です。自然対数の偏角(-\pi,\pi] にとるものとして、

\begin{eqnarray}\displaystyle\zeta_4^M(s)&=&\frac{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\zeta_4^M(s)\\&=&\frac{1+i^{-s}+(-1)^{-s}+(-i)^{-s}}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\zeta_4^M(s)\\&=&\frac{1}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\sum_{m,n\in\mathbb{Z}\atop{(m,n)\neq{(0,0)}}}^{} (m+ni)^{-s}\\&=&\frac{1}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\sum_{m,n\in\mathbb{Z}\atop{(m,n)\neq{(0,0)}}}^{} (m+ni)^{-s}\\&=&\frac{1}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\left(\sum_{m,n\in\mathbb{Z}\atop{n\neq{0}}}^{} (m+ni)^{-s}+\sum_{0\neq\in\mathbb{Z}}^{} m^{-s}\right)\\&=&\frac{1}{(1+e^{-\frac{i\pi}{2}}+e^{-i\pi s}+e^{\frac{i\pi}{2}})}\left( (1+(-1)^{-s})\sum_{m\in\mathbb{Z},n\geq{1}}^{} (m+ni)^{-s}+(1+(-1)^{-s})\zeta(s)\right)\\&=&\frac{1+e^{-i\pi s}}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\left(\zeta(s)+\sum_{m\in\mathbb{Z}}^{} \sum_{n=1}^{\infty} (m+ni)^{-s}\right)\\&=&\frac{1+e^{-i\pi s}}{(1+e^{-\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{i\pi s}{2}})}\left(\zeta(s)+\sum_{n=1}^{\infty} \frac{(-2\pi i)^s}{\Gamma(s)}\sum_{m=1}^{\infty} m^{s-1}e^{-2\pi mn}\right)\end{eqnarray}

となるので、和の変数を mn から n におきかえて整理すると
\begin{eqnarray}\displaystyle \zeta_4^M(s)&=&\frac{1+e^{i\pi s}}{(1+e^{\frac{i\pi s}{2}}+e^{-i\pi s}+e^{\frac{3i\pi s}{2}})}\left(\zeta(s)+\frac{(-2\pi i)^s}{\Gamma(s)}\sum_{n=1}^{\infty} \sigma_{s-1}(n)e^{-2\pi n}\right)\end{eqnarray}

を得ます。これは通常発散するはずの値も求めることができて、今回使うのは s=0 の場合です(一般の負の整数に対しても明示的に求めることが可能です。1/7数学カフェで詳細をお話します)。簡単に計算できて

\begin{eqnarray}\displaystyle\zeta_4^M(0)=-\frac{1}{4}\end{eqnarray}

となります。この値を後で使います(っていうかこの記事で \zeta_4^M を導入したのはこれを使うためだけです)。


これはまだ準備段階の更に最初の方で、ここから延々と計算のキツイ補題が続いていくのですが、一つの記事に収めるのは到底不可能なのでここではいくつかをFactとして認めることにします(あまりにも記述量が膨大になってしまうことに投稿一日前に気づきました。ごめんなさい。今度書きます)。

とりあえず、明解にするために一番大きな定理をここで出してしまいます:



[主定理(新谷、1980)]

\begin{eqnarray}\displaystyle\rho_2(1,z)\rho_2(1,-z)=(2\pi)^{\frac{3}{2}}z^{-\frac{1}{2}}\eta(z)\exp\left(\pi i\left(\frac{1}{4}+\frac{1}{12z}\right)\right)\end{eqnarray}

ここで\begin{eqnarray}\displaystyle\eta(z):=e^{\frac{\pi iz}{12}}\prod_{n=1}^{\infty} (1-e^{2\pi inz})\end{eqnarray}


まずこれを使って目標である \rho_2(1,i) を先に求めてしまいます(膨大な記述、というのはこの定理の証明のことです): z=i として、

\begin{eqnarray}\displaystyle\rho_2(1,i)\rho_2(1,-i)&=&(2\pi)^{\frac{3}{2}}e^{-\frac{i\pi}{4}}\eta(i)e^{\frac{i\pi}{4}+\frac{\pi}{12}}\\&=&(2\pi)^{\frac{3}{2}}\eta(i)e^{\frac{\pi}{12}}\end{eqnarray}

となり、それなりに簡単に書けました。

次に、また記号を改めた関数を導入します:

\begin{eqnarray}\displaystyle\zeta^{\star}(s,a):=\sum_{m,n\geq{0}\atop{(m,n)\neq{(0,0)}}}^{} (m+na)^{-s}\end{eqnarray}

これを使うと、今問題にしている \rho_2(1,i), \rho_2(1,-i) が次のように書けます:

\begin{eqnarray}\displaystyle\rho_2(1,i)&=&\exp\left(-\left.\frac{d}{ds}\zeta^{\star}(s,i)\right|_{s=0}\right)\\\rho_2(1,-i)&=&\exp\left(-\left.\frac{d}{ds}\zeta^{\star}(s,-i)\right|_{s=0}\right)\end{eqnarray}

そして、

\begin{eqnarray}\displaystyle\zeta^{\star}(s,-i)&=&\sum_{m,n\geq{0}\atop{(m,n)\neq{(0,0)}}}^{} (m-in)^{-s}\\&=&(-i)^{-s}\sum_{m,n\geq{0}\atop{(m,n)\neq{(0,0)}}}^{} (n+im)^{-s}\\&=&e^{\frac{i\pi s}{2}}\zeta^{\star}(s,i)\end{eqnarray}

となることより、

\begin{eqnarray}\displaystyle\rho_2(1,i)\rho_2(1,-i)&=&\exp\left(-\left.\frac{d}{ds}\left(\zeta^{\star}(s,i)+\zeta^{\star}(s,-i)\right)\right|_{s=0}\right)\\&=&\exp\left(-\left.\frac{d}{ds}(1+e^{\frac{i\pi s}{2}})\zeta^{\star}(s,i)\right|_{s=0}\right)\\&=&\exp\left(-\left.\frac{i\pi}{2}e^{\frac{i\pi s}{2}}\zeta^{\star}(s,i)+(1+e^{\frac{i\pi s}{2}})\left(\frac{d}{ds}\zeta^{\star}(s,i)\right)\right|_{s=0}\right)\end{eqnarray}

というのが導かれます。また、

\begin{eqnarray}\displaystyle\zeta^{\star}(s,i)&=&\left(\sum_{m\geq{1},n\geq{0}}^{} (m+ni)^{-s}\right)+\left(\sum_{n=1}^{\infty} (ni)^{-s}\right)&=&\zeta_M^4+e^{-\frac{i\pi s}{2}}\zeta(s)\end{eqnarray}

なのでもちろん \displaystyle \zeta^{\star}(0)=-\frac{1}{4}-\frac{1}{2}=-\frac{3}{4} です。ここから

\begin{eqnarray}\displaystyle\rho_2(1,i)\rho_2(1,-i)&=&e^{\frac{3i\pi}{8}}\exp\left(-\left.2\left(\frac{d}{ds}\zeta^{\star}(s,i)\right)\right|_{s=0}\right)\\&=&e^{\frac{3i\pi}{8}}\rho_2(1,i)^2\end{eqnarray}

がわかって、さっき導いた \rho_2(1,i), \rho_2(1,-i) の式と合わせて

\begin{eqnarray}\displaystyle e^{\frac{3i\pi}{8}}\rho_2(1,i)^2=(2\pi)^{\frac{3}{2}}\eta(i)e^{\frac{\pi}{12}}\end{eqnarray}

が出ます。 \rho_2(1,i) を中心に整理すると

\begin{eqnarray}\displaystyle \rho_2(1,i)=(2\pi)^{\frac{3}{4}}\sqrt{\eta(i)}e^{\frac{\pi}{24}-\frac{3i\pi}{16}}\end{eqnarray}

となりますね。まだ残っている \eta(i) が気になりますが、これは

o-v-e-r-h-e-a-t.hatenablog.com

の定理1と保型形式間に成り立つ関係式

\begin{eqnarray}\displaystyle\eta(z)^{24}=\frac{E_4^3-E_6^2}{1728}\end{eqnarray}

からすぐに値を求めることができます(この関係式については

tsujimotter.hatenablog.com

も参照)。実際には \displaystyle\eta(i)=\frac{\Gamma\left(\frac{1}{4}\right)}{2\pi^{\frac{3}{4}}} なので、結局

\begin{eqnarray}\displaystyle\rho_2(1,i)=2^{\frac{1}{4}}\pi^{\frac{3}{8}}\sqrt{\Gamma\left(\frac{1}{4}\right)}e^{\frac{\pi}{24}-\frac{3i\pi}{16}}\end{eqnarray}

となります。



さて、一番面倒な主定理の証明ですが、長すぎるのでかなり天下り的になってしまいました。ご了承ください(後ろから読んでいくのを推奨します)。先に記号をいくつか定義します:

\begin{eqnarray}\displaystyle\zeta_2(s,w,z):&=&\sum_{m,n\geq{0}}^{} (m+nz+w)^{-s}\\\Gamma^{\star}(w,z):&=&\exp\left(\left.\frac{\partial}{\partial s}\zeta_2(s,w,z)\right|_{s=0}\right)\\\Gamma_2(w,z):&=&\rho_2(1,z)\Gamma^{\star}(w,z)\\LG(w):&=&\frac{1}{2\pi i}\int_{I(\lambda,\infty)}^{} \frac{e^{-wt}}{1-e^{-t}}\frac{\log t}{t^2}dt-\frac{\gamma-\pi i}{2}B_2(w)\end{eqnarray}

積分I(\lambda,\infty)(-\infty,\infty) から \mathbb{C} への写像 \phi:(-\infty,\infty)\rightarrow\mathbb{C} として

I(\lambda,\infty):\phi(u)=\begin{cases}

u & (u<-\lambda) \\
\lambda\exp\left(\pi i\frac{u+\lambda}{\lambda}\right) & (-\lambda\geq{u}\geq{\lambda}) \\
u & (u>\lambda)
\end{cases}

として定義されます。イメージとしては無限の長さを持つ丸底フラスコを右に倒した感じです。

B_n(w) はベルヌーイ多項式です。定義は B_1=\frac{1}{2} になるほう。


以下の補題3,4,5,6をひとまず認めることにします:



[補題3]


(1)
\displaystyle LG(w)=\frac{B_2(w)}{2}-\left.\frac{\partial}{\partial s}\zeta(s,w)\right|_{s=-1}

(2)
\displaystyle \frac{d}{dw}LG(w)=\log\frac{\sqrt{2\pi}}{\Gamma(w)}



[補題4]


(1)

\begin{eqnarray}\displaystyle LG(n)=-\frac{n^2}{2}\log n+\frac{3}{4}n^2-B_1(n\log n-n)-\frac{B_2}{2}\log n+c+O(n^{-1})\end{eqnarray}

(2)

\begin{eqnarray}\displaystyle LG(w+nz)=-\frac{B_2(w+nz)}{2}\log nz+\frac{3}{4}n^2z^2+nzB_1(w)+O(n^{-1})\end{eqnarray}



[補題5]


\begin{eqnarray}\displaystyle\log\Gamma^{\star}(w,z)=\frac{1}{z}LG(w)-B_1\log\frac{\Gamma(w)}{\sqrt{2\pi}}-\frac{B_2}{2}\psi(w)+O(w^{-1})\end{eqnarray}



[補題6]


(1)

\begin{eqnarray}\displaystyle\log\frac{\Gamma(w+nz)}{\sqrt{2\pi}}=\left(w+nz-\frac{1}{2}\right)\log n-n+\frac{\sqrt{2\pi}}+O(n^{-1})\end{eqnarray}

(2)

\begin{eqnarray}\displaystyle\psi(n)=\log n+O(n^{-1})\end{eqnarray}


補題6はスターリングの公式なので簡単に証明できるはずですが、補題3,4,5はどれも面倒です(3は比較的簡単。

補題7以降は証明を与えます。



[補題7]

(1)
\begin{eqnarray}\displaystyle\frac{\Gamma_2(w+z,z)}{\Gamma_2(w,z)}=\frac{\sqrt{2\pi}}{\Gamma(w)}\end{eqnarray}

(2)
\begin{eqnarray}\displaystyle\Gamma_2(1,z)=\sqrt{\frac{2\pi}{z}}\end{eqnarray}

[証明]

(1)

\begin{eqnarray}\displaystyle\frac{\Gamma_2(w+z,z)}{\Gamma_2(w,z)}&=&\frac{\Gamma^{\star}(w+z,z)}{\Gamma^{\star}(w,z)}\\&=&\exp\left(\left.\frac{\partial}{\partial s}\zeta_2(s,w+z,z)-\zeta_2(s,w,z)\right|_{s=0}\right)\end{eqnarray}

\begin{eqnarray}\displaystyle\zeta_2(s,w+z,z)-\zeta_2(s,w,z)&=&\sum_{m,n\geq{0}}^{} \left( (m+nz+w+z)^{-s}-(m+nz+w)^{-s}\right)\\&=&\left(\sum_{m\geq{1},n\geq{0}}^{} (m+nz+w)^{-s}\right)-\left(\sum_{m,n\geq{0}}^{} (m+nz+w)^{-s}\right)\\&=&-\sum_{m=0}^{\infty} (m+w)^{-s}\\&=&-\zeta(s,w)\end{eqnarray}

あとはレルヒの公式でおわり。

(2)

o-v-e-r-h-e-a-t.hatenablog.com

より

\begin{eqnarray}\displaystyle\frac{\Gamma_2(w+1,z)}{\Gamma_2(w,z)}=\frac{\sqrt{2\pi}}{\Gamma\left(\frac{w}{z}\right)}z^{\frac{1}{2}-\frac{w}{z}}\end{eqnarray}

がいえるので

\begin{eqnarray}\displaystyle\frac{\Gamma_2(w,z)}{\Gamma\left(\frac{w}{z}\right)}&=&\frac{\frac{1}{w}+O(1)}{\frac{z}{w}+O(1)}\\&=&\frac{1+O(w)}{z+O(w)}\end{eqnarray}

より、上の式で w\rightarrow{0} として \displaystyle\Gamma_2(1,z)=\sqrt{2\pi}\frac{1}{z}z^{\frac{1}{2}}=\sqrt{\frac{2\pi}{z}} となります。

[証明終]




[補題8]


\displaystyle\begin{align*}& \log\Gamma^{\star}(w,z)-\log\frac{\Gamma(w)}{\sqrt{2\pi}}-\sum_{m=1}^{n-1} \left(\log\frac{\Gamma(w+mz)}{\sqrt{2\pi}}-\left(mz+w-\frac{1}{2}\right)+mz-\frac{B_2(w)}{2mz}\right)\\& =\log\Gamma^{\star}(w+nz,z)+\left(w-\frac{1}{2}\right)\log\Gamma(n)+z\sum_{m=1}^{n-1} \left( (m\log m-m)+\frac{n(n-1)}{2}\right)z\log z+\frac{B_2(w)}{2z}\sum_{m=1}^{n-1} \frac{1}{m}\end{align*}

[証明]

補題7より

\displaystyle \log\Gamma^{\star}(w,z)-\log\Gamma^{\star}(w+z,z)=\log\frac{\Gamma(w)}{\sqrt{2\pi}}.

これを繰り返し使って

\begin{eqnarray}\displaystyle\log\Gamma^{\star}(w,z)&=&\log\Gamma^{\star}(w+nz,z)+\sum_{m=0}^{n-1} \frac{\Gamma(w+mz)}{\sqrt{2\pi}}\\&=&\log\frac{\Gamma(w)}{\sqrt{2\pi}}+\log\Gamma^{\star}(w+nz,z)\\&+&\sum_{m=1}^{n-1} \left(\log\frac{\Gamma(w+mz)}{\sqrt{2\pi}}-\left(mz+w-\frac{1}{2}\right)\log mz+mz-\frac{B_2(w)}{2mz}\right)\\&+&\sum_{m=1}^{n-1} \left(\left(mz+w-\frac{1}{2}\right)\log mz-mz+\frac{B_2(w)}{2mz}\right)\end{eqnarray}

となって、

\displaystyle\begin{align*}& \log\Gamma^{\star}(w,z)-\log\frac{\Gamma(w)}{\sqrt{2\pi}}-\sum_{m=1}^{n-1} \left(\log\frac{\Gamma(w+mz)}{\sqrt{2\pi}}-\left(mz+w-\frac{1}{2}\right)+mz-\frac{B_2(w)}{2mz}\right)\\& =\log\Gamma^{\star}(w+nz,z)+\frac{B_2(w)}{2z}\left(\sum_{m=1}^{n-1} \frac{1}{m}\right)+z\left(\sum_{m=1}^{n-1} (m\log m-m)\right)+\left(w-\frac{1}{2}\right)\log ( (n-1)!z)+z\sum_{m=1}^{n-1} mz\log z\\& =\log\Gamma^{\star}(w+nz,z)+\left(w-\frac{1}{2}\right)\log\Gamma(n)+z\sum_{m=1}^{n-1} \left( (m\log m-m)+\frac{n(n-1)}{2}\right)z\log z+\frac{B_2(w)}{2z}\sum_{m=1}^{n-1} \frac{1}{m}.\end{align*}

[証明終]



[補題9]

\displaystyle\log\Gamma^{\star}(w,z)-\log\frac{\Gamma(w)}{\sqrt{2\pi}}-\sum_{m=1}^{\infty} \left(\log\frac{\Gamma(w+mz)}{\sqrt{2\pi}}-\left(mz+w-\frac{1}{2}\right)+mz-\frac{B_2(w)}{2mz}\right)
\displaystyle =\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+LG(1)z-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma-\frac{B_2(w)}{2z}\log z+\left(w-\frac{1}{2}\right)\frac{\log z}{2}-\frac{z\log z}{12}

[証明]

補題8に補題4,5,6を代入する:

\begin{align*}\displaystyle & \log\Gamma^{\star}(w,z)-\log\frac{\Gamma(w)}{\sqrt{2\pi}}-\sum_{m=1}^{\infty} \left(\log\frac{\Gamma(w+mz)}{\sqrt{2\pi}}-\left(mz+w-\frac{1}{2}\right)+mz-\frac{B_2(w)}{2mz}\right)\\& =\frac{1}{z}LG(w+nz)+\frac{1}{2}\log\frac{\Gamma(w+nz)}{\sqrt{2\pi}}-\frac{z}{12}\psi(w+nz)+O(n^{-1})\\& +\left(w-\frac{1}{2}\right)\left(\left(n-\frac{1}{2}\right)\log n-n+\log\sqrt{2\pi}+O(n^{-1})\right)\\& +zLG(1)-zLG(n)+\frac{n(n-1)}{2}z\log z-\left(w-\frac{1}{2}\right)(n-1)\log z\\& +\frac{B_2(w)}{2z}\left(\log n+\gamma+O(n^{-1})\right)\\& =\frac{1}{z}\left(-\frac{1}{2}B_2(w+nz)\log nz+\frac{3}{4}n^2z^2+nzB_1(w)+O(n^{-1})\right)\\& +\frac{1}{2}\left(\left(nz+w-\frac{1}{2}\right)\log nz-nz+O(n^{-1})\right)-\frac{z}{12}(\log nz+O(n^{-1}))\\& +\left(w-\frac{1}{2}\right)\left(\left(n-\frac{1}{2}\right)\log n-n+\log\sqrt{2\pi}+O(n^{-1})\right)\\& +zLG(1)-z\left(-\frac{n^2}{2}\log n+\frac{3}{4}n^2+\frac{1}{2}(n\log n-n)-\frac{1}{12}\log n+O(n^{-1})\right)\\& +\frac{n(n-1)}{2}z\log z+\left(w-\frac{1}{2}\right)(n-1)\log z+\frac{B_2(w)}{2z}\log n+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\& =-\frac{z}{2}B_2(w+nz)\log nz+\frac{3}{4}n^2z+nw-\frac{n}{2}+\frac{1}{2}\left(nz+w-\frac{1}{2}\right)\log nz-\frac{nz}{2}-\frac{z\log nz}{12}\\& +\left(w-\frac{1}{2}\right)\left(n-\frac{1}{2}\right)\log n-\left(w-\frac{1}{2}\right)n+\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)+\frac{n^2z}{2}\log n-\frac{3}{4}n^2z\\& -\frac{z}{2}(n\log n-n)+\frac{z}{12}\log n+\frac{n(n-1)}{2}z\log z+\left(w-\frac{1}{2}\right)(n-1)\log z+\frac{B_2(w)}{2z}\log n+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\& =\left(-\frac{1}{2z}B_2(w+nz)+\frac{1}{2}nz+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}\right)\log nz+\left(\left(w-\frac{1}{2}\right)\left(n-\frac{1}{2}\right)+\frac{n^2z}{2}+\frac{12}{z}+\frac{B_2(w)}{2z}\right)\log n\\& +\left(w-\frac{1}{2}-\frac{z}{2}-\left(w-\frac{1}{2}\right)\right)n-\frac{z}{2}(n\log n-n)+\frac{n(n-1)}{2}z\log z\\& +\left(w-\frac{1}{2}\right)n\log z+\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\& =\left(-\frac{1}{2z}B_2(w+nz)+\frac{1}{2}nz+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}\right)\log nz+\left(\left(w-\frac{1}{2}\right)\left(n-\frac{1}{2}\right)+\frac{n^2z}{2}+\frac{12}{z}+\frac{B_2(w)}{2z}\right)\log n\\& +\left(w-\frac{1}{2}-\frac{z}{2}-\left(w-\frac{1}{2}\right)\right)n-\frac{nz}{2}\log nz+\frac{n^2}{2}z\log z\\& +\left(w-\frac{1}{2}\right)n\log z+\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\& =\left(-\frac{1}{2z}B_2(w+nz)+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}+\frac{n^2z}{12}\right)\log z\\& +\left(-\frac{1}{2z}B_2(w+nz)+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}+\left(w-\frac{1}{2}\right)\left(n-\frac{1}{2}\right)+\frac{n^2z}{2}+\frac{z}{12}+\frac{B_2(w)}{2z}\right)\log n\\& +\left(w-\frac{1}{2}\right)n\log z+\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\& =\left(-\frac{1}{2z}B_2(w+nz)+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}+\frac{n^2z}{12}\right)\log z\\& +\left(-wn-\frac{n^2z}{2}+\frac{n}{2}+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}+wn-\frac{n}{2}-\frac{w}{2}+\frac{1}{4}+\frac{n^2z}{2}+\frac{z}{12}\right)\log n\\& +\left(w-\frac{1}{2}\right)n\log z+\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\& =\left(-\frac{w^2}{2z}-nw-\frac{n^2z}{2}+\frac{w}{2z}+\frac{n}{2}-\frac{1}{12z}+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}+\frac{n^2z}{12}+wn-\frac{n}{2}\right)\log z\\& +\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\ &=\left(-\frac{B_2(w)}{2z}+\frac{w}{2}-\frac{1}{4}-\frac{z}{12}\right)\log z+\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+zLG(1)-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma+O(n^{-1})\\ &\overset{n\rightarrow\infty}{\rightarrow}\left(w-\frac{1}{2}\right)\log\sqrt{2\pi}+LG(1)z-\left(w-\frac{1}{2}\right)\log z+\frac{B_2(w)}{2z}\gamma-\frac{B_2(w)}{2z}\log z+\left(w-\frac{1}{2}\right)\frac{\log z}{2}-\frac{z\log z}{12}.\end{align*}

[証明終]




[補題10]

\begin{eqnarray}\displaystyle\rho_2(z)=(2\pi)^{\frac{3}{4}}\exp\left(-\frac{\gamma}{12z}-\frac{12}{z}+z\zeta'(-1)+\left(\frac{z}{12}-\frac{1}{4}+\frac{1}{12z}\right)\log z\right)\\\times\prod_{n=1}^{\infty} \frac{\sqrt{2\pi}}{\Gamma(1+nz)}\exp\left(\frac{1}{12nz}+\left(\frac{1}{2}+nz\right)\log nz-nz\right)\end{eqnarray}



[証明]

補題9で w=1 とすると

\begin{eqnarray}\displaystyle\log\Gamma^{\star}(1,z)-\log\frac{1}{\sqrt{2\pi}}-\sum_{m=1}^{\infty} \left(\log\frac{\Gamma(mz+1)}{\sqrt{2\pi}}-\left(mz+\frac{1}{2}\right)\log mz+mz-\frac{1}{12mz}\right)\\=\frac{\log 2\pi}{4}+LG(1)z-\frac{\log z}{2}+\frac{\gamma}{12z}-\frac{\log z}{12z}+\frac{\log z}{4}-\frac{z\log z}{12}.\end{eqnarray}

\displaystyle LG(1)=\frac{1}{12}-\zeta'(-1) より

\begin{eqnarray}\displaystyle\log\Gamma^{\star}(1,z)=-\frac{\log 2\pi}{4}+\frac{z}{12}-\zeta'(-1)z-\frac{\log z}{2}+\frac{\gamma}{12z}-\frac{\log z}{12z}+\frac{\log z}{4}-\frac{z\log z}{12}\\+\sum_{m=1}^{\infty} \left(\log\frac{\Gamma(mz+1)}{\sqrt{2\pi}}-\left(mz+\frac{1}{2}\right)\log mz+mz-\frac{1}{12mz}\right).\end{eqnarray}

補題7(2)より \displaystyle \Gamma_2(1,z)=\sqrt{\frac{2\pi}{z}}なので\Gamma^{\star} の定義より \displaystyle \Gamma^{\star}(1,z)=\frac{1}{\rho_2(1,z)}\sqrt{\frac{2\pi}{z}} となる、よって \displaystyle \log\Gamma^{\star}(1,z)=\frac{\log 2\pi}{2}-\frac{\log z}{2}-\log\rho_2(1,z) となって

\begin{eqnarray}\displaystyle -\log\rho_2(1,z)=-\frac{3\log 2\pi}{4}+\frac{z}{12}-\zeta'(-1)z+\frac{\gamma}{12z}-\frac{\log z}{12z}+\frac{\log z}{4}-\frac{z\log z}{12}\\+\sum_{m=1}^{\infty} \left(\log\frac{\Gamma(mz+1)}{\sqrt{2\pi}}-\left(mz+\frac{1}{2}\right)\log mz+mz-\frac{1}{12mz}\right).\end{eqnarray}

[証明終]

12/24の記事は盛田みずすまし氏の記事です、お楽しみに。