8
現代数学解説
文献あり

アペリーの定理とWZ method

951
0

はじめに

 この記事では WZ method の応用としてApéryの数列
an=k=0n(nk)2(n+kk)2cn,kbn=k=0n(nk)2(n+kk)2(cn,k=m=1n1m3+m=1k(1)m12m3(nm)(n+mm))
に関する謎について考察していきます。

Apéryの定理

  WZ methodの記事 ではApéryの定理

Apéryの定理

 ζ(3)は無理数である。

の証明に使われた公式
ζ(3)=52n=1(1)n1n3(2nn)
がWZ methodによって導出できることを紹介しました。
 しかしApéryの定理の証明にはこれ以外にもWZ的な現象が深く関わっています。

Apéryの手法

 まずApéryの手法について簡単に紹介しておきましょう。
 発想としてはζ(3)に急速に収束する有理数列pn/qnを構成することで次の補題を適用することにあります。

 実数βと整数列pn,qnに対し
βpnqn(nZ0),limn(qnβpn)=0
が成り立つとき、βは無理数である。

 ある整数p,qを用いてβ=p/qと表せるとすると、十分大きいあるnに対し
0<|qnppnq|<1
が成り立つことになるが、これはqnppnqが整数であることに矛盾。よってβは無理数でなければならない。

 まあこのような有理数列pn/qnを簡単に構成できれば苦労はしません。それが未だにできていないから我々は数々の定数の無理性を確かめられずにいるわけです。
 しかしApéryは冒頭で紹介した奇妙な数列
an=k=0n(nk)2(n+kk)2cn,kbn=k=0n(nk)2(n+kk)2(cn,k=m=1n1m3+m=1k(1)m12m3(nm)(n+mm))
を持ち出すことによりβ=ζ(3)に対しそれを成し遂げました。以下でこの数列によってζ(3)の無理性を示せることを見ていきましょう。

 un=an,bnは同じ漸化式
(n+1)3un+1(2n+1)(17n2+17n+5)un+n3un1=0
を満たす。ただし初項は
a0=0,a1=6,b0=1,b1=5
である。

 ちなみにbnの初期値についてはb1=5の代わりにb1=0によって特徴付けることもあります。

 後で説明する。

anbn1an1bn=6n3
が成り立つ。

 後で説明する。

limnanbn=ζ(3)
が成り立つ。

|m=1k(1)m12m3(nm)(n+mm)|m=112m3(nm)(n+mm)=m=112m3n2=ζ(3)2n2
と評価できるので
|ζ(3)anbn|(ζ(3)m=1n1m3)+ζ(3)2n20(n)
を得る。

0<ζ(3)anbn<6ζ(3)bn2
が成り立つ。

An=ζ(3)anbn
とおいたとき補題4,5から
An1An=6n3bn1bn,limNAN=0
が成り立つこと注意すると
An=limN(AnAN)=limNk=n+1N(Ak1Ak)=k=n+16k3bk1bk
と表せるので、bnの単調増加性に注意すると
0<An6bn2k=n+11k3<6ζ(3)bn2
を得る。

 ここでbnは常に整数であるのに対しanは整数とは限らないことに注意し、Dnanが整数となるような数列Dnを取ると
0<Dnbnζ(3)Dnan<6ζ(3)Dnbn
と評価できるので補題2を適用するためには
limnDnbn=0
が成り立つことを示さなければなりません。が、このあたりの議論についてはこの記事の本題から反れるので詳細については折り畳みにて解説しておきます。

dn=lcm(1,2,3,,n)
とおくと
2dn3(n+kk)cn,kZ
が成り立つ。特に2dn3anは整数となる。

証明

 各mに対し
dn3(n+kk)m3(nm)(n+mm)Z
が成り立つこと、特に任意の素数pに対し
ordp(dn3(n+kk)m3(nm)(n+mm))0
となることを示せばよい。
 いまLegendreの公式(やKummerの定理)から任意の素数pに対し
ordp(nm)logpnordpm=ordpdnordpm
と評価できるので
(n+mm)(n+kk)=(n+m)!k!(n+k)!m!=(km)(n+kkm)
に注意すると
ordp(m3(nm)(n+mm)(n+kk))ordp(m3(nm)(km))3ordpm+ordpdn+ordpdk2ordpm=3ordpm+ordpdn+ordpdk3ordpdn
を得る。

 ζ(3)は無理数である。

証明

 上での議論と補題7から
limndn3bn=0
が成り立つことを示せばよい。
 いま素数定理から十分大きい任意のnに対して
dn=pnplogpnpnn=nπ(n)<n(1+ε)nlogn=e(1+ε)n(<3n)
のように評価でき、またbnの満たす漸化式の特性方程式
x234x+1=0
を解くと この記事 にて紹介したPerronの定理から
bn(17+122)n(>33n)
と評価できるので、十分大きい任意のnに対し
dn3bn<27n33n0(n)
が成り立つ。よって主張を得る。

ζ(2)の無理性

 また興味深いことに同様の数列
an=k=0n(nk)2(n+kk)cn,kbn=k=0n(nk)2(n+kk)(cn,k=m=1n(1)m1m2+m=1k(1)m+n12m2(nm)(n+mm))
を用いることで
η(2)=m=1(1)m1m2=12ζ(2)(=π212)
の無理性についても同じ議論で示すことができます。
 やることは上と同じなので対応する補題の内容だけ紹介しておきます。

 un=an,bnは同じ漸化式
(n+1)2un+1(11n2+11n+3)unn2un1=0
を満たす。ただし初項は
a0=0,a1=5/2,b0=1,b1=3
である。

anbn1an1bn=(1)n152n2
が成り立つ。

limnanbn=η(2)
が成り立つ。

0<|η(2)anbn|<5ζ(2)2bn2
が成り立つ。

 2dn2anは整数となる。

 ζ(2)は無理数である。

Apéry数列の謎

 見ての通りApéryの手法には多くの謎が秘められています。例えば簡単に思い浮かぶもので言えば

  • なぜan,bnan,bnは同じ漸化式を満たすのか。
  • なぜanbn1an1bn=6/n2のような関係が成り立つのか。
  • なぜ
    (n+1)3un+1(2n+1)(17n2+17n+5)un+n3un1=0
    のような漸化式のある有理数列による二解の比(Apéry limits)がζ(3)ζ(2)に収束するのか。
  • なぜ
    (n+1)3un+1(2n+1)(17n2+17n+5)un+n3un1=0
    のような漸化式の有理数列による解の分母の増大度が階乗n!ではなく指数関数dne(1+ε)nで抑えられるのか。
  • そもそもan,bnという数列はどこから持ってきたのか。

といった疑問が挙げられます。
 これらの現象には一体どのような一般論が眠っているのか、あるいは限られた場合にしか成り立たない特殊な現象なのかが非常に気になるところではありますが、残念ながらまだその全てが解明されているわけではなさそうです。
 しかし上二つの現象についてはある程度確かなことがわかっているので今回の記事ではそのことについて簡単に考察していこうと思います。

カゾラーティアン

 まず簡単に解決できる問題である
anbn1an1bn=6n2
という等式について考察していきましょう。
 素朴な方法としてはan,bnが同じ漸化式
(n+1)3an+1(2n+1)(17n2+17n+5)an+n3an1=0(n+1)3bn+1(2n+1)(17n2+17n+5)bn+n3bn1=0
を満たすことからこの両辺にそれぞれbn,anを掛けて差を取ることで
(n+1)3wn+1n3wn=0(wn=anbn1an1bn)
がわかるので
wn+1=n3(n+1)3wn=w1(n+1)3=6(n+1)3
と示せます。
 正直この問題についてはこれで十分なのですが、少し味気ないのでもう少し一般的な主張を紹介しておきましょう(かなり余談なので読み飛ばしてもらっても構いません)。

Casoratian

 d+1項間線形漸化式
k=0dpk(n)u(n+k)=0
の解u1(n),u2(n),,ud(n)に対しCasoratianw(n)
w(n)=det(u1(n)u2(n)ud(n)u1(n+1)u2(n+1)ud(n+1)u1(n+d1)u2(n+d1)ud(n+d1))
によって定める。

 なおd+1項間漸化式は、簡単のためpd(n)=1とし
u(n)=(u(n)u(n+1)u(n+d1))T
とおくと
u(n+1)=(010001p0(n)p1(n)pd1(n))u(n)
という一項間漸化式に帰着できるのでもう少し一般に次のような行列式をCasoratianと言うことにします。

 d次元の線形漸化式
u(n+1)=P(n)u(n)
の解u1(n),u2(n),,ud(n)に対しCasoratianw(n)
w(n)=det(u1(n)u2(n)ud(n))
と定める。

 このとき以下が成り立ちます。

 Casoratianの一般項は
w(n)=det(P(n1)P(n2)P(0))w(0)
と表せる。特に
u1(n),u2(n),,ud(n)が線形独立常にw(n)0u1(n),u2(n),,ud(n)が線形従属常にw(n)=0
が成り立つ。

w(n+1)=det(P(n))w(n)
を示せばよい。そのことは
w(n+1)=det(u1(n+1)u2(n+1)ud(n+1))=det(P(n)u1(n)P(n)u2(n)P(n)ud(n))=det(P(n))det(u1(n)u2(n)ud(n))
とわかる。

 ちなみにCasoratianは線形微分方程式におけるWronskianの類似物であり、Wronskianの関係式
W(t)=exp(0ttrA(s)ds)W(0)
と比較すると面白いですね。
 またこれをd+1項間漸化式に関する話に引き戻すと以下が成り立ちます。

 d+1項間線形漸化式
k=0dpk(n)u(n+k)=0
に関するCasoratianの一般項は
w(n)=(1)dnw(0)k=0n1p0(k)pd(k)
と表せる。

 例えばζ(2)の無理性証明に出てきた数列an,bn
(n+1)2un+1(11n2+11n+3)unn2un1=0
という漸化式を満たしていたためwn=anbn1an1bn
wn+1=n2(n+1)2wn=(1)nw1(n+1)2=(1)n52(n+1)2
と求まっていたわけです。

漸化式とポテンシャル

 次に今回の記事の本題である、なぜ
an=k=0n(nk)2(n+kk)2cn,kan=k=0n(nk)2(n+kk)cn,kbn=k=0n(nk)2(n+kk)2bn=k=0n(nk)2(n+kk)(cn,k=m=1n1m3+m=1k(1)m12m3(nm)(n+mm))(cn,k=m=1n(1)m1m2+m=1k(1)m+n12m2(nm)(n+mm))
によって定まる数列an,bnan,bnが同じ漸化式を満たすのか、ということについて考察していきましょう。

ポテンシャル

 ここで WZ methodの記事 において
ζ(3)=52n=1(1)n1n3(2nn)
の証明に出てきたWZ-pair
H(n,k)=(1)k1(nk)(n+kk)F(n,k)=12(k+1)3H(n,k+1)G(n,k)=nk(k+1)2(n+1)2H(n,k+1)
を思い出しましょう。これを用いると上のcn,k
cn,k=m=1nG(m1,0)+m=1kF(n,m1)
と表せます(cn,kζ(2)に関するWZ-pairを用いて同様に表せます)。
 またグリッド上の経路
C:(0,0)(1,0)(n,0)(n,1)(n,2)(n,k)
を考えるとこれはWZ form
ω=F(n,k)δk+G(n,k)δn
の経路和分
cn,k=Cω
として表せます。
 このようにcn,kはいわゆる"渦なし"差分形式の(0,0)から(n,k)までの経路和分となっており、これはベクトル解析で言うところのポテンシャルに類似するものであることがわかります。

 数列c(n,k)がWZ form
ω=F(n,k)δk+G(n,k)δn
ポテンシャルであるとは、その偏差分が
c(n,k+1)c(n,k)=F(n,k)c(n+1,k)c(n,k)=G(n,k)
を満たすことを言う。

a(n),b(n)の満たす漸化式

 さてcn,kにポテンシャルという名前が付いたことにより次のような一般的な問いを立てることができます。

 properな超幾何数列b(n,k)とWZ form
ω=F(n,k)δk+G(n,k)δn
のポテンシャルc(n,k)に対して
a(n)=k=0nb(n,k)c(n,k)b(n)=k=0nb(n,k)
とおくと、これらの満たす漸化式はどのような関係にあるだろうか。

 そしてこの問題はZeilbergerによって次のような解答が与えられました。

 上の問題設定においてある多項式R,Sおよび
B(n,k)b(n,k),D(n,k)b(n,k)F(n,k)
n,kについての有理関数とするような超幾何数列B,Dであって
R(n,N)b(n,k)=B(n,k+1)B(n,k)S(n,N)R(n,N)(b(n,k)c(n,k))=A(n,k+1)A(n,k)(A(n,k)=S(n,N)(B(n,k)c(n,k))+D(n,k))
を満たすようなものが存在する。ただしNnについての前進作用素Nf(n,k)=f(n+1,k)とした。

  Zeilberger's method により
R(n,N)b(n,k)=B(n,k+1)B(n,k)
を満たすようなR,Bの存在はわかる。
 このとき
R(n,N)=i=0Iri(n)Ni
とおくと
=R(n,N)(b(n,k)c(n,k))(R(n,N)b(n,k))c(n,k)=i=0Iri(n)b(n+i)(c(n+i,k)c(n,k))=i=0Iri(n)b(n+i)j=0i1(c(n+j+1,k)c(n+j,k))=i=0Iri(n)b(n+i)j=0i1G(n+j,k)=(R(n,N)b(n,k))c(n,k)=(B(n,k+1)B(n,k))c(n,k)=B(n,k+1)c(n,k+1)B(n,k)c(n,k)B(n,k+1)(c(n,k+1)c(n,k))=B(n,k+1)c(n,k+1)B(n,k)c(n,k)B(n,k+1)F(n,k)
が成り立つので
H(n,k)=i=0Iri(n)b(n+i)j=0i1G(n+j,k)B(n,k+1)F(n,k)
とおくと
R(n,N)(b(n,k)c(n,k))=B(n,k+1)c(n,k+1)B(n,k)c(n,k)+H(n,k)
と表せる。
 このとき
H(n,k)b(n,k)F(n,k)=i=0Iri(n)b(n+i)b(n)j=0i1G(n+j,k)F(n,k)B(n,k+1)b(n,k)
n,kについての有理関数となることに注意してZeilberger's methodにより
S(n,N)H(n,k)=D(n,k+1)D(n,k)
なるS,Dを取ることで、主張のようなR,S,B,Dが得られる。

 そして適当な条件下でこの等式
R(n,N)b(n,k)=B(n,k+1)B(n,k)S(n,N)R(n,N)(b(n,k)c(n,k))=A(n,k+1)A(n,k)
をそれぞれkについて足し合わせることで
R(n,N)b(n)=0,S(n,N)R(n,N)a(n)=0
というa(n),b(n)についての漸化式が得られることとなります。
 特にこのことからa(n),b(n)は同じ漸化式
S(n,N)R(n,N)a(n)=S(n,N)R(n,N)b(n)=0
を満たすこともわかります。

残された謎

 これによりApéryの数列
an=k=0n(nk)2(n+kk)2cn,kan=k=0n(nk)2(n+kk)cn,kbn=k=0n(nk)2(n+kk)2bn=k=0n(nk)2(n+kk)(cn,k=m=1n1m3+m=1k(1)m12m3(nm)(n+mm))(cn,k=m=1n(1)m1m2+m=1k(1)m+n12m2(nm)(n+mm))
が同じ漸化式を満たすことはある程度必然であることがわかりました。
 しかし依然として不思議なのがan,bnが"丁度"同じ漸化式を満たす、つまり上の定理におけるS(n,N)としてS(n,N)=1が取れる、あるいは同じことですが上の証明におけるH(n,k)kについてGosper総和可能となる、ということにあります。与えられたポテンシャルc(n,k)に対して都合のいいb(n,k)が存在するのか、あるいは与えられたb(n,k)に対して都合のいいポテンシャルc(n,k)が存在するのか、はたまた上の例が特殊なだけで一般論なんて存在しないのか。やはりまだまだ疑問は尽きませんね。

おわりに

 以上がApéryの定理の証明とその手法の背景に関する考察でした。
 今回の記事では個人的に疑問であった「なぜan,bnのような奇妙な数列が同じ漸化式を満たすのか」という問題に一つの解答が与えられて満足していますが(ちなみにCasoratianの項は思い付きで書いたおまけです)、上でも言及したようにApéryの手法にはまだまだ興味深い謎が数多く残されています。皆さんも興味があればApéryの手法について考察してみたり、他の定数の無理性証明への応用を考えてみてはいかがでしょうか。
 とりあえず今回はこんなところで。では。

おまけ:log2の無理性

 上でも登場したZeilbergerの文献にApéryの手法によるlog2の無理性証明が載っていたので簡単に紹介しておこうと思います。

H(n,k)=(1)k2nn!k!(n+k)!=(1)n2k(n+kk)F(n,k)=1k+1H(n,k+1)G(n,k)=12(k+1)H(n,k+1)
とおくとF,GはWZ-pairをなす。

an=k=0n(nk)(n+kk)cn,kbn=k=0n(nk)(n+kk)(cn,k=m=1k(1)m1m+m=1n(1)k2mm(k+mm))
とおくとun=an,bnは同じ漸化式
(n+1)un+13(2n+1)un+nun1=0
を満たす。ただし初項は
a0=0,a1=2,b0=1,b1=3
である。

 あとの議論は
log2=n=1(1)n1n
に注意すればわかる(と思う)ので省略します。
 ちなみにanは実は
an=k=0n(nk)(n+kk)m=1k(1)m1m
と表せる、つまり以下の事実が成り立つことが示せます。

k=0n(nk)(n+kk)m=1n(1)k2mm(k+mm)=0
が成り立つ。

証明

 多項式列Pn(x)
Pn(x)=1n!dndxn(x(1x))n=k=0n(1)k(nk)(n+kk)xk
によって定めると上の式は
=k=0n(nk)(n+kk)m=1n(1)k2mm(k+mm)=m=1n12mk=0n(1)k(nk)(n+kk)k!(m1)!(k+m)!=m=1n12mk=0n(1)k(nk)(n+kk)01xk(1x)m1dx=m=1n12m01(1x)m1Pn(x)dx
と表せる。
 また部分積分を繰り返すことで任意の関数fに対し
01Pn(x)f(x)dx=(1)nn!01(x(1x))ndndxnf(x)dx
が成り立つことと、1mnに対し
dndxn(1x)m1=0
となることに注意すると主張を得る。

参考文献

[1]
R. Apéry, Irrationalité de ζ(2) et ζ(3), Astérisque, 1979, 11 - 13
[2]
A. van der Poorten, A proof that Euler missed..., The Mathematical Intelligencer, 1979, 195 - 203
投稿日:2024322
更新日:202483
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

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

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. Apéryの定理
  3. Apéryの手法
  4. ζ(2)の無理性
  5. Apéry数列の謎
  6. カゾラーティアン
  7. 漸化式とポテンシャル
  8. ポテンシャル
  9. a(n),b(n)の満たす漸化式
  10. 残された謎
  11. おわりに
  12. おまけ:log2の無理性
  13. 参考文献