7
現代数学解説
文献あり

ラマヌジャン・佐藤級数を理解したい(その2)

246
0
$$\newcommand{a}[0]{\alpha} \newcommand{Aut}[0]{\operatorname{Aut}} \newcommand{b}[0]{\beta} \newcommand{C}[0]{\mathbb{C}} \newcommand{d}[0]{\delta} \newcommand{dis}[0]{\displaystyle} \newcommand{e}[0]{\varepsilon} \newcommand{F}[4]{{}_2F_1\left(\begin{matrix}#1,#2\\#3\end{matrix};#4\right)} \newcommand{farc}[2]{\frac{#1}{#2}} \newcommand{G}[0]{\Gamma} \newcommand{g}[0]{\gamma} \newcommand{Gal}[0]{\operatorname{Gal}} \newcommand{H}[0]{\mathbb{H}} \newcommand{id}[0]{\operatorname{id}} \newcommand{Im}[0]{\operatorname{Im}} \newcommand{Ker}[0]{\operatorname{Ker}} \newcommand{l}[0]{\left} \newcommand{L}[0]{\Lambda} \newcommand{la}[0]{\lambda} \newcommand{La}[0]{\Lambda} \newcommand{Li}[0]{\operatorname{Li}} \newcommand{li}[0]{\operatorname{li}} \newcommand{M}[4]{\begin{pmatrix}#1& #2\\#3& #4\end{pmatrix}} \newcommand{N}[0]{\mathbb{N}} \newcommand{o}[0]{\omega} \newcommand{O}[0]{\Omega} \newcommand{ol}[1]{\overline{#1}} \newcommand{ord}[0]{\operatorname{ord}} \newcommand{P}[0]{\mathfrak{P}} \newcommand{p}[0]{\mathfrak{p}} \newcommand{q}[0]{\mathfrak{q}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{r}[0]{\right} \newcommand{R}[0]{\mathbb{R}} \newcommand{Re}[0]{\operatorname{Re}} \newcommand{s}[0]{\sigma} \newcommand{t}[0]{\theta} \newcommand{ul}[1]{\underline{#1}} \newcommand{vp}[0]{\varphi} \newcommand{vt}[0]{\vartheta} \newcommand{Z}[0]{\mathbb{Z}} \newcommand{z}[0]{\zeta} \newcommand{ZZ}[1]{\mathbb{Z}/#1\mathbb{Z}} \newcommand{ZZt}[1]{(\mathbb{Z}/#1\mathbb{Z})^\times} $$

はじめに

 この記事では前回の記事に引き続きラマヌジャン・佐藤級数についての論文Chan, Cooper(2012)を読んでいきます。

Clausenの公式の類似

 Chan-Chan-Liuの手法では
$$A_n=(-1)^n\sum^n_{k=0}\binom nk^2\binom{2n-2k}{n-k}\binom{2k}k$$
のような数列に対し直接
$$Z=\sum^\infty_{n=0}A_nX^n$$
という級数を考えていたが、超幾何関数型の円周率公式はまず
$$z=\sum^\infty_{n=0}\frac{(\frac1s)_n(1-\frac1s)_n}{(1)_n^2}x^n$$
という級数を考えてからClausenの公式
$$\l(\sum^\infty_{n=0}\frac{(\frac1s)_n(1-\frac1s)_n}{(1)_n^2}x^n\r)^2 =\sum^\infty_{n=0}\frac{(\frac12)_n(\frac1s)_n(1-\frac1s)_n}{(1)_n^3}(4x(1-x))^n$$
を用いることで導出できるのであった。
 Chan-Cooperの手法ではこれの次のような一般化を考えている。

Clausenの公式の類似

 $a,b,c$を整数とし、数列$\{s_n\},\{t_n\}$$s_{-1}=t_{-1}=0,s_0=t_0=1$および
\begin{align*} (n+1)^2s_{n+1}&=(an^2+an+b)s_n+cn^2s_{n-1}\\ (n+1)^3t_{n+1}&=-(2n+1)(an^2+an+a-2b)t_n-(4c+a^2)n^3t_{n-1} \end{align*}
によって定める。このとき
\begin{align*} (1+cx^2)\l(\sum^\infty_{n=0}s_nx^n\r)^2 &=\sum^\infty_{n=0}\binom{2n}ns_n\l(\frac{x(1-ax-cx^2)}{(1+cx^2)^2}\r)^n\\ (1-ax-cx^2)\l(\sum^\infty_{n=0}s_nx^n\r)^2 &=\sum^\infty_{n=0}t_n\l(\frac x{1-ax-cx^2}\r)^n \end{align*}
が成り立つ。

Clausenの公式との関係

 上に挙げたClausenの公式はこの第一式の$c=0$の場合として得られる。実際
$$s_n=\frac{(\frac1s)_n(1-\frac1s)_n}{(1)_n^2}$$
とおいたとき
$$(n+1)^2s_n=\l(\frac1s+n\r)\l(1-\frac1s+n\r)s_n=(n^2+n+\frac1s\l(1-\frac1s\r))s_n$$
が成り立つので上の定理が適用でき($a,b,c$が整数でなくても上の定理は成り立つ)、
$$\farc{(\frac12)_n}{(1)_n}=\frac{(2n)!}{2^{2n}(n!)^2}=\binom{2n}n\frac1{4^n}$$
に注意すると
\begin{align*} \l(\sum^\infty_{n=0}s_nx^n\r)^2 &=\sum^\infty_{n=0}\binom{2n}ns_n\l(x(1-x)\r)^n\\ &=\sum^\infty_{n=0}\frac{(\frac12)_n}{(1)_n}s_n\l(4x(1-x)\r)^n \end{align*}
が得られる。
 ちなみに下では$b$を整数にするために適当な$N_s$を取って$s'_n=N_s^ns_n$と置くわけだが、具体的には
$$\l(\frac1m\r)_n\l(\frac2m\r)_n\cdots\bigg(\frac mm\bigg)_n=\frac{(mn)!}{m^{mn}}$$
が成り立つことに注意すると
\begin{alignat*}{2} s=2:s_n&=\frac1{(2^{2n})^2}\frac{((2n)!)^2}{(n!)^4} &&=\frac1{16^n}\binom{2n}n^2\\ s=3:s_n&=\frac1{3^{3n}}\frac{(3n)!}{(n!)^3} &&=\frac1{27^n}\binom{3n}n\binom{2n}n\\ s=4:s_n&=\frac{2^{2n}}{4^{4n}}\frac{(4n)!}{(n!)^2(2n)!} &&=\frac1{64^n}\binom{4n}{2n}\binom{2n}n\\ s=6:s_n&=\frac{2^{2n}3^{3n}}{6^{6n}}\frac{(6n)!}{n!(2n)!(3n)!} &&=\frac1{432^n}\binom{6n}{3n}\binom{3n}n \end{alignat*}
と表せるので
$$N_2=16,\quad N_3=27,\quad N_4=64,\quad N_6=432$$
と取ることになる。

 やることとしては両辺が同じ微分方程式と同じ初期値を満たすことを確認するだけである。いま
$$z=\sum^\infty_{n=0}s_nx^n$$
とおいたときこれの満たす微分方程式は$\vt=x\frac{d}{dx}$を用いて
$$\vt^2z=x(a\vt^2+a\vt+b)z+cx^2(\vt+1)^2z$$
と表せ、これを展開すると
$$x(1-ax-cx^2)z''+(1-2ax-3cx^2)z'-(b+cx)z=0$$
となることに注意する。
 またこれは$x=0,\infty$および$1-ax-cx^2=0$を確定特異点に持つFuchs型微分方程式であることにも注意したい。このことの意味についてはこの記事の1,2節で解説している。
 以降の議論は単純に計算するだけであるが、手計算でどうにかするのは非常に面倒なので省略する。

 このような公式の発見的な方法についても中々興味深い理論が広がっているようだが、まだ私の理解が浅いためここで解説できる段階にはないのが悲しいところ。いつかちゃんと理解してみたい。

$$X=\frac{x(1-ax-cx^2)}{(1+cx^2)^2},\quad Y=\frac x{1-ax-cx^2}$$
とおいたとき、上のような数列$s_n,t_n$に対し
\begin{align*} \sum^\infty_{n=0}t_nY^{n+\frac12}&=\sum^\infty_{n=0}\binom{2n}ns_nX^{n+\frac12}\\ \sum^\infty_{n=0}t_n(n+\frac12)Y^{n+\frac12} &=\sqrt{1-4aX-16cX^2}\sum^\infty_{n=0}\binom{2n}ns_n(n+\frac12)X^{n+\frac12} \end{align*}
が成り立つ。

 第一式については
$$\sqrt{\frac XY}=\frac{1-ax-cx^2}{1+cx^2}$$
に注意すると上の定理からわかる。第二式については
$$X=\frac{Y}{(1+aY)^2+4cY^2}$$
および
$$1-4aX-16cX^2=\l(\frac{1-(a^2+4c)Y^2}{(1+aY)^2+4cY^2}\r)^2$$
に注意して第一式を$Y$について微分することでわかる。

級数の収束性について

 いま上の定義だけだと級数
$$\sum^\infty_{n=0}s_nx^n,\quad\sum^\infty_{n=0}t_nx^n$$
$x=0$の近傍で収束するのか、つまり$s_n,t_n$は急速に増大してしまわないのか、という疑問がある。
 これについては漸化式に対するポアンカレやペロンの定理として知られる以下の結果が役に立つ。

 数列$\{x_n\}$についての漸化式
$$x_{n+k}+\sum^{k-1}_{j=0}p^{(j)}_nx_{n+j}=0$$
であって、その係数$p_n^{(j)}$が極限値$p_n^{(j)}\to p^{(j)}$を持つようなものをポアンカレ型の漸化式という。
 またポアンカレ型の漸化式に対して定まる方程式
$$\la^k+\sum^{k-1}_{j=0}p^{(j)}\la^j=0$$
のことを特性方程式と言う。

Perronの第二定理

 ポアンカレ型の漸化式
$$x_{n+k}+\sum^{k-1}_{j=0}p^{(j)}_nx_{n+j}=0$$
に対し、その解$\{x_n\}$は次のいずれかを満たす。
・ある$N$より先で常に$x_n=0$が成り立つ。
・特性方程式のある解$\la$が存在して$\dis\limsup_{n\to\infty}\sqrt[n]{|x_n|}=|\la|$が成り立つ。

 この結果によると$s_n,t_n$が満たす漸化式についての特性方程式はそれぞれ
$$\la^2-a\la-c=0,\quad\la^2+2a\la+(a^2+4c)=0$$
であり、このことから
\begin{align*} \limsup_{n\to\infty}\sqrt[n]{|s_n|}&\leq\max\l|\frac{a\pm\sqrt{a^2+4c}}2\r|=:M_s\\ \limsup_{n\to\infty}\sqrt[n]{|t_n|}&\leq\max|a\pm2\sqrt{-c}|=:M_t \end{align*}
と評価できる。したがってコーシーの収束判定法より
$$\sum^\infty_{n=0}s_nx^n,\quad\sum^\infty_{n=0}t_nx^n$$
はそれぞれ$|x|<1/M_s,1/M_t$で収束することが保証されるわけである。
 また
$$\limsup_{n\to\infty}\sqrt[n]{\binom{2n}n}=4$$
が成り立つことにも注意したい。

余談

 ちなみにペロンの定理の一般化として次のような定理が知られている。

Coffmanの定理

 $A$$k\times k$行列、$f_n(y):\R^k\to\R^k$を連続関数とし$y_n\in\R^k$についての漸化式
$$y_{n+1}=Ay_n+f_n(y_n)$$
を考える。このとき$f_n$
$$\frac{\|f_n(y)\|}{\|y\|}\to0\quad\big((n,y)\to(\infty,0)\big)$$
を満たすとすると
・ある$N$より先で$y_n\neq0$
$y_n\to0\quad(n\to\infty)$
なる解$y_n$に対しある$A$の固有値$\la$が存在して
$$\lim_{n\to\infty}\sqrt[n]{\|y_n\|}=|\la|<1$$
が成り立つ。

 これは
$$\frac{dy}{dt}=Ay+g(t,y)\qquad\l(\lim_{(t,y)\to(\infty,0)}\frac{\|g(t,y)\|}{\|y\|}=0\r)$$
なる微分方程式を考えたとき、$A$の任意の固有値$\la$に対し$\|y(t)\|\approx|e^{\la t}|$なる漸近挙動を持つ解$y(t)$が取れる、という感じの有名な事実の漸化式への類似となっている。
 これも中々興味深い話題ではあるのだが、足を踏み入れるには少し深そうな気がするのでまた気が向いたときに勉強しようと思う。

$s_n,t_n$の明示形

 あとはある整数$a,b,c$によって定まる$s_n$を用いて
$$z(\tau)=\sum^\infty_{n=0}s_nx(\tau)^n$$
と表せるような保型形式$z(\tau)$と保型関数$x(\tau)$の組を持ってくれば、上の公式とChan-Chan-Liuの手法から
\begin{align*} \frac1\pi &=\sum^\infty_{n=0}\binom{2n}ns_n(An+B)X^n\\ &=\sum^\infty_{n=0}t_n(A'n+B')Y^n\\ \end{align*}
という円周率公式が構成できるわけである。
 この記事では一旦具体的な$z(\tau)$$x(\tau)$を見ていく前に、どのような$a,b,c$に対し$s_n$$t_n$はどのように表せるのかを見ていこう。Chan-Cooperでは以下の表が提示されている($\ell$は対応するラマヌジャン・佐藤級数のレベルを表す。具体的には次回の記事で解説する)。

\begin{array}{cccc}\hline \ell&(a,b,c)&s_n&t_n\\\hline 1&(432,60,0)&\dis\binom{6n}{3n}\binom{3n}n &\dis\sum^n_{k=0}s_n\binom{2k}k\binom{n+k}{n-k}(-a)^{n-k}\\ 2&(64,12,0)&\dis\binom{4n}{2n}\binom{2n}n &\dis\sum^n_{k=0}s_n\binom{2k}k\binom{n+k}{n-k}(-a)^{n-k}\\ 3&(27,6,0)&\dis\binom{3n}n\binom{2n}n &\dis\sum^n_{k=0}s_n\binom{2k}k\binom{n+k}{n-k}(-a)^{n-k}\\ 4&(16,4,0)&\dis\binom{2n}n^2 &\dis(-1)^n\sum^n_{k=0}\binom{2k}k^2\binom{2n-2k}{n-k}\\ 5&(11,3,1)&\dis\sum^n_{k=0}\binom nk^2\binom{n+k}k &\dis\sum^n_{k=0}(-1)^{n-k}\binom nk^3\binom{4n-5k}{3n}\\ 6_A&(-17,-6,-72)&\dis\sum_{j\leq k}\binom nk\binom kj^3(-8)^{n-k} &\dis\sum^n_{k=0}\binom nk^2\binom{n+k}k^2\\ 6_B&(10,3,9)&\dis\sum^n_{k=0}\binom nk^2\binom{2k}k &\dis(-1)^n\sum^n_{k=0}\binom nk^2\binom{2n-2k}{n-k}\binom{2k}k\\ 6_C&(7,2,8)&\dis\sum^n_{k=0}\binom nk^3 &\dis\sum^n_{k=0}\binom{n+k}n\binom n{3k}\binom{3k}{2k}\binom{2k}k(-3)^{n-3k}\\ 8&(12,4,-32)&\dis\sum^n_{k=0}\binom n{2k}\binom{2k}k^24^{n-2k} &\dis(-1)^n\sum^n_{k=0}\binom nk^2\binom{2k}n^2\\ 9&(-9,-3,-27)&\dis\sum^n_{k=0}\binom nk\binom{n-k}k\binom{n-2k}k(-3)^{n-3k} &\dis\sum_{j\leq k}\binom nk^2\binom nj\binom kj\binom{k+j}n\\\hline \end{array}

 ラマヌジャンの考えた円周率公式
$$\frac1\pi=\sum^\infty_{n=0}\frac{(\frac12)_n(\frac1s)_n(1-\frac1s)_n}{(1)_n^3}\frac{An+B}{C^n}\quad(s=2,3,4,6)$$
はこの$\ell=1,2,3,4$$s_n$についてのものであり、Chan-Chan-Liuの発見した円周率公式
$$\frac1\pi=\frac{\sqrt3}8\sum^\infty_{n=0}(\sum^n_{k=0}\binom nk^2\binom{2n-2k}{n-k}\binom{2k}k)\frac{5n+1}{64^n}$$
はこの$\ell=6_B$$t_n$についてのものとなる。
 また佐藤猛 氏の発表した円周率公式
$$\frac{\sqrt{15}}\pi=6\sum^\infty_{n=0}(\sum^n_{k=0}\binom nk^2\binom{n+k}k^2)(20n+10-3\sqrt5)\l(\frac{\sqrt5-1}2\r)^{12n+6}$$
$\ell=6_A$$t_n$についてのものであり、佐藤氏が同時に発表していた円周率公式
$$\frac{\sqrt5}\pi=3\sum^\infty_{n=0}(-1)^nt_n(282n+141-50\sqrt5)\l(\frac1{5\sqrt5}\l(\frac{\sqrt5-1}2\r)^{15}\r)^{n+\frac12}\quad((a,b,c)=(11,3,1))$$
$\ell=5$$t_n$の場合となっている。
 しかし Wikipedia で提示されているものを見るに円周率公式を構成できる$(a,b,c)$の取り方はこれだけに尽きないようにも思えるが、どうなのだろうか。このことが気にかかったので一先ずどのような$(a,b,c)$に対してどのように保型形式と保型関数の組$z(\tau),x(\tau)$を探せばいいのか、ということについて次回の記事で見ていこうと思う。

参考文献

[2]
S. Elaydi, An Introduction to Difference Equations, Third edition, Springer, New York, 2005, pp. 387-389
投稿日:20231122
更新日:20231226
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

子葉
子葉
985
226009
主に複素解析、代数学、数論を学んでおります。 私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。 同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中