Loading [MathJax]/jax/output/CommonHTML/jax.js

【ζ6】フルヴィッツゼータ関数のHermiteの公式(積分表示・Abel-Planaの和公式・ビネの公式)(ゼータ関数の基礎6)

「ゼータ関数の基礎」シリーズ第6回です。前回はゼータ関数の特殊値は関数等式について解説しました:

読まなくても本記事は理解できます。

今日のテーマと概要

エルミートの公式

Hermite's formula for ζ(s)

0<a1 とする。あらゆる sC に対しζ(s,a)=12as+a1ss1+20sin(sarctanya)(a2+y2)s2(e2πy1)dyこれを「エルミートの公式」という。

これとビネの第2公式を比較することによりlims1(ζ(s,a)1s1)=ψ(a)ただし ψ(a) はディガンマ(digamma)関数。

またエルミートの公式によりフルヴィッツゼータ関数の微分係数が求まる。dζ(s,a)ds|s=0=logΓ(a)12log2π

公式の導出について調べるとビネの第2公式から出発するものが出てきますが、ここではWhittaker-Watsonのテキストに習ってアベル・プラナ(Abel-Plana)の和公式からスタートします。どっちみちビネの第2公式も使いますので本質は同じなのかもしれません。

古いですが有名な書物で、どんどん改訂版が出ています。前半は解析学一般、後半は特殊関数という内容で、網羅的に勉強できます。演習問題に解答がないのが昔ながらのものって感じ。2022/11/6現在、最新版は5th Editionで私も所有していますが、廉価な3rdとかでも十分かと。


A Course of Modern Analysis: fifth Edition


A Course of Modern Analysis: Third Edition

予備知識:Abel-Planaの和公式

Abel–Plana summation formula

nk=mf(k)=f(m)+f(n)2+nmf(z)dz+i0f(m+iy)f(miy)f(n+iy)+f(niy)e2πy1dy

これをある関数に適用するところからスタートします。おさえておきましょう。

過去に「ビネの第2公式」を証明するときに、補題としてこれを導出しました:

関数の設定と和公式の適用

複素数 s=σ+it とし σ>1 とします。定数を 0<a1 とし、次の複素関数 f(z) を定義します。f(z)1(z+a)sこの f(z) はAbel-Planaに適用するための条件limyf(x±iy)e2πy=0を明らかに満たしています。

和公式で m=0 , n=N とします。ここで N は大きな自然数です。Nk=01(k+a)s=f(0)+f(N)2+N0dz(z+a)s+i0(a+iy)s(aiy)se2πy1dyi0(a+N+iy)s(a+Niy)se2πy1dy=12as+12(a+N)s+N0dz(z+a)s+i0(a+iy)s(aiy)se2πy1dyi0(a+N+iy)s(a+Niy)se2πy1dy

N の極限をとりましょう。左辺はゼータ関数になります。ζ(s,a)=12as+0dz(z+a)s+i0(a+iy)s(aiy)se2πy1dyilimN0(a+N+iy)s(a+Niy)se2πy1dy右辺第2項は σ>1 に注意して0dz(z+a)s=[(z+a)1ss1]0=a1ss1第3項は分子の2項を極形式にします。(a+iy)s(aiy)s=(a2+y2)s2(eiϕseiϕs)(tanϕ=ya)=2isin(sarctanya)(a2+y2)s2第4項の分子に関しても同様に(a+N+iy)s(a+Niy)s=2isin(sarctanya+N)[(a+N)2+y2]s2これらを(2)に用いると

ζ(s,a)=12as+a1ss1+20sin(sarctanya)(a2+y2)s2(e2πy1)dy2limN0sin(sarctanya+N)[(a+N)2+y2]s2(e2πy1)dy

右辺の最後の項に残った積分は N の極限で 0 となります。それを示しましょう。

極限で消える積分

では0sin(sarctanya+N)[(a+N)2+y2]s2(e2πy1)dyを評価します。分子が厄介です。

分子の評価

s=σ+it , σ>1 に対し|sin(sarctanya+N)|2=sinh2(|t|arctanya+N)+sin2(σarctanya+N)

|sin(a+ib)|2=sin2a+sinh2b を用いました。導出は以下です。|sin(a+ib)|2=|sinacosib+cosasinib|2=|sinacoshb+icosasinhb|2=sin2acosh2b+cos2asinh2b=sin2a(1+sinh2b)+cos2asinh2b=sin2a+sinh2b

不等式その1

y0=(a+N)tanπ2σ とします。するとσarctany0a+N=π2が成立します。

0yy0 ならば 0σarctanya+Nπ2 です。また N を大きくとれば ya+Nπ2 とできます。arctanXX も合わせれば0σarctanya+Nσya+Nπ2[0,π2] での sinx の単調増加性からsin2(σarctanya+N)sin2σya+Nsin2XsinXX ですので結局次のようになります。sin2(σarctanya+N)σya+N

なお y0 はその定義より N の1次式となっています。

不等式その2

t=0 のときはこの節は無視します。t0 として考えましょう。arctanXX の関係および sinhx の単調増加性からsinh2(|t|arctanya+N)sinh2|t|ya+Nここで次のような y10 をとります。sinh2|t|y1a+N=|t|y1a+N数値としては |t|y1a+N0.81 です。

sinhXX

すると 0yy1 なる y についてsinh2|t|ya+N|t|ya+Nsinh2(|t|arctanya+N)|t|ya+Nなお y1 はその定義より N の1次式となっています。

0へ収束

N に依存しない定数 KK=min[tanπ2σ,1|t|]とします。ただし t=0 のときは後者を無視します。すると 0Ky0 かつ Ky1 です。よって 0yK に対して(5)(6)を(4)に同時に適用できます。|sin(sarctanya+N)|2(σ+|t|)ya+N|sin(sarctanya+N)|σ+|t|a+Ny

いよいよ積分がゼロとなることを示します。|0sin(sarctanya+N)[(a+N)2+y2]s2(e2πy1)dy|0|sin(sarctanya+N)|[(a+N)2+y2]σ2(e2πy1)dyK0σ+|t|a+Ny[(a+N)2+y2]σ2(e2πy1)dy+K|sin(sarctanya+N)|[(a+N)2+y2]σ2(e2πy1)dyσ+|t|a+N1(a+N)σK0ydye2πy1+1(a+N)σK|sin(sarctanya+N)|e2πy1dy

ここで|sin(a+bi)|=sin2a+sinh2b1+sinh2b=|coshb|および arctanXπ2 を用いて|0sin(sarctanya+N)[(a+N)2+y2]s2(e2πy1)dy|σ+|t|a+N1(a+N)σ0ydye2πy1+1(a+N)σK|coshπt2|e2πy1dy=σ+|t|a+N1(a+N)σ0ydye2πy1+coshπt2(a+N)σKdye2πy1

第1項は0ydye2πy1=π2(2π)32ζ(32)M1なので N に関係なく有界です。

第2項はKdye2πy1=K(e2πye2πy11)dy=[12πloge2πy1e2πy]K=12πlog11e2πKM2なので N に関係なく有界です。

以上より|0sin(sarctanya+N)[(a+N)2+y2]s2(e2πy1)dy|M1σ+|t|a+N1(a+N)σ+M2coshπt2(a+N)σN0

公式の完成

長くなりましたが(3)および先ほどの結果から

ζ(s,a)=12as+a1ss1+20sin(sarctanya)(a2+y2)s2(e2πy1)dy

なお右辺の積分の収束性については先の積分評価と同様にやると示すことができます。y0=atanπ2σ として・・・のように。

この公式で s=0 としてみましょう。すると ζ(0,a)=12a となります。さらに a=1 はリーマンのゼータ関数に一致しますから ζ(0)=12 です。過去記事の既習事項と一致します。

s=1の極限を調べる

Binetの第2公式を応用する

(7)の積分はあらゆる s で値をとりますから、ζ(s,a)s=1 に1位の極をもつと分かります。そこで極の原因となる 1s1 を左辺に追いやって極限をとってみます。lims1(ζ(s,a)1s1)=12a+lims1a1s1s1+20sin(arctanya)a2+y2(e2πy1)dy=12aloga+20y(a2+y2)(e2πy1)dy

ここで次の関係式を用いました。sin(arctanya)=ya2+y2

Binetの第2公式に関連する等式(※)ddzlogΓ(z)=12z+logz20tdt(z2+y2)(e2πy1)と比較すると

lims1(ζ(s,a)1s1)=ddalogΓ(a)=ψ(a)

を得ます。

a=1 とするとリーマンゼータ関数になりますので

lims1(ζ(s)1s1)=ψ(1)=γ

ζ(s,a)のs=0における微分係数

エルミートの公式(7)を s で微分します。dζ(s,a)ds=loga2asa1slogas1a1s(s1)2+20[arctanyacos(sarctanya)(a2+y2)s2(e2πy1)12log(a2+y2)sin(sarctanya)(a2+y2)s2(e2πy1)]dys0 とするとdζ(s,a)ds|s=0=loga2+alogaa+20arctanyae2πy1dy=(a12)logaa+20arctanyae2πy1dy

Binetの第2公式(前述の記事参照)logΓ(z)=(z12)logzz+12log2π+20arctantze2πt1dtを用いれば次の結論を得ます。

dζ(s,a)ds|s=0=logΓ(a)12log2π

a=1 とするとリーマンゼータ関数となります。

ζ(0)=12log2π

私としてはなかなか計算がヘビーでした。

次はシンプルなお話:

応援のおねがい

Please support me!

まめしば
まめしば

記事を気に入って下さった方、「応援してあげてもいいよ」という方がいらっしゃったら15円から可能なので支援していただければ幸いです。情報発信を継続していくため、サーバー維持費などに充てさせていただきます。

ご支援いただいた方は、こちらで確認できます。

Amazonギフトの場合、
Amazonギフト券- Eメールタイプ – Amazonベーシック
より、金額は空白欄に適当に(15円から)書きこんで下さい。受取人は「mamekebiamazonあっとgmail.com」です(あっとは@に置き換えてください)。贈り主は「匿名」等でOKです。全額がクリエイターに届きます。

OFUSEは登録不要で、100円から寄付できます。金額の90%がクリエイターに届きます。

codocは登録不要で、100円から寄付できます。金額の85%がクリエイターに届きます。

Remaining:
codocの場合

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA