14
大学数学基礎解説
文献あり

多重S値 (MSV)について、いろいろ

828
0
$$\newcommand{abs}[1]{\left |#1\right |} \newcommand{bm}[1]{\boldsymbol{#1}} \newcommand{C}[0]{\mathbb{C}} \newcommand{Fourier}[2]{\mathcal{F}_{#1}\left [#2\right ]} \newcommand{Hartley}[2]{\mathcal{H}_{#1}\left [#2\right ]} \newcommand{Hilbert}[2]{\mathcal{Hil}_{#1}\left [#2\right ]} \newcommand{inttrans}[3]{\mathcal{#1}_{#2}\left [#3\right ]} \newcommand{invtrans}[3]{\mathcal{#1}^{-1}_{#2}\left [#3\right ]} \newcommand{Laplace}[2]{\mathcal{L}_{#1}\left [#2\right ]} \newcommand{Li}[0]{\operatorname{Li}} \newcommand{Mellin}[2]{\mathcal{M}_{#1}\left [#2\right ]} \newcommand{ord}[0]{\operatorname{ord}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{Res}[1]{\underset{#1}{\operatorname{Res}}} \newcommand{sh}[0]{ш} \newcommand{tLaplace}[2]{\mathcal{B}_{#1}\left [#2\right ]} \newcommand{Weierstrass}[2]{\mathcal{W}_{#1}\left [#2\right ]} $$

先日、Harryさんのツイートを発端として、ある種の三角関数を含んだ重積分(Dirichlet積分の多重化)について色々なことが判明したので、現在わかっていることをまとめます。
インデックス表記については、多重ゼータ値などと統一します。多重S値という名前は、僕が適当につけたものです。

多重S値

$\varepsilon_1,\cdots,\varepsilon_n\in\{0,1\}$に対して、
$$I(\varepsilon_1,\cdots,\varepsilon_n):=\int \cdots \int_{0< t_1<\cdots< t_n}\prod_{i=1}^nA_{\varepsilon_i}(t_i)dt_i$$
とする。ここで、$$A_0(t)=\frac{\cos t}{t},A_1(t)=\frac{\sin t}{t}$$である。このとき、任意のインデックス$\bm k=(k_1,\cdots,k_r)$に対して、
$$ S(\bm k):=I(1,\{0\}^{k_1-1},\cdots,1,\{0\}^{k_r-1}) $$
は収束し、これを多重S値(MSV)と呼ぶ。

$$ \begin {aligned} S(1)&=\int _{0}^{\infty }\frac {\sin x}xdx\\ &=\frac {\pi }2\quad \left (\text{Dirichlet 積分}\right )\\ S(1,2)&=\int _{0< x< y< z}\frac {\sin x}x\frac {\sin y}y\frac {\cos z}zdxdydz\\ &=\int _{0}^{\infty }\frac {\cos z}zdz\int _{0}^{z}\frac {\sin y}zdy\int _{0}^{y}\frac {\sin x}xdx\\ &=\int _{0}^{\infty }\frac {\sin x}xdx\int _{x}^{\infty }\frac {\sin y}ydy\int _{y}^{\infty }\frac {\cos z}zdz \end {aligned} $$

関係式

MSV同士の積について、積分範囲のシャッフルを考えることで以下が成り立ちます。

シャッフル関係式

$$ S(\bm k \sh \bm l)=S(\bm k)S(\bm l) $$

実は、関係式はこれ以外にも存在します。複素積分を使います。
まず複素積分の使い方について具体例から見ていきましょう。
$$ \begin {aligned} f(z)&=\frac {e^{iz}}z\int _{0}^{z}\frac {\sin x}xdx \end {aligned} $$
とします。$f(z)$$\C$全体で正則で、上半平面において、$(\Im(z)^2+1)|f(z)|$は有界になり、$R$に依らない定数$M$以下とできます。
従って、$0$を中心として、$R$から$Ri$へ向かう中心角$\dfrac{\pi}2$の経路を$C$とすると、
$$ \begin {aligned} \left |\int _Cf(z)dz\right | &\leq \int _{0}^{\frac {\pi }2}\left |f(Re^{i\theta })Re^{i\theta }\right |d\theta \\ &\leq \int _{0}^{\frac {\pi }2}\frac {RM}{R^2\sin ^{2}\theta +1}d\theta \\ &\underset{R\to \infty }{\to }0 \end {aligned} $$
となります。コーシーの積分定理を適用すれば
$$ \begin {aligned} \int _{0}^{\infty }f(z)dz&=\int _{0}^{i\infty }f(z)dz\\ &=\int _{0}^{\infty }\frac {e^{-y}}{y}dy\int _{0}^{iy}\frac {\sin x}xdx\\ &=i\int _{0< x< y}\frac {e^{-y}}y\frac {\sinh x}xdxdy \end {aligned} $$
実部,虚部を比較することにより
$$ \begin {aligned} \int _{0< x< y}\frac {\sin x}x\frac {\cos y}ydxdy&=0\\ \int _{0< x< y}\frac {\sin x}x\frac {\sin y}ydxdy&=\int _{0< x< y}\frac {\sinh x}x\frac {e^{-y}}ydxdy \end {aligned} $$
を得ます。なんと$S(2)=0$が求まっちゃいました!!!
ではもう一つの$S(1,1)$についてなんですが、この積分表示を利用せずともシャッフル積から
$$ \begin {aligned} S(1,1)&=\frac {1}2S(1)^2=\frac {\pi ^{2}}8 \end {aligned} $$
と計算することができます。しかし結果に着目すれば$\zeta(2)=\pi^2/6$の有理数倍になっていますしリーマンゼータ値と関係がありそうです。この積分を使うとS(1,1)をリーマンゼータ値と結びつけることができます。
$$ \begin {aligned} \int _{0< x< y}\frac {\sinh x}x\frac {e^{-y}}ydxdy &=\int _{0< s,t}dsdt\int _{0< x< y}e^{-xs}\sinh xe^{-yt}e^{-y}dxdy\\ &=\int _{0< s,t}dsdt\int _{0< y< x<1}x^{s-1}\frac {x^{-1}-x}2y^{t}dxdy\quad \left (e^{-x}\mapsto x,e^{-y}\mapsto y\right )\\ &=\int _{0< s,t}\frac {dsdt}{2(t+1)}\int _{0< x<1}\left (x^{s-2}-x^{s}\right )x^{t+1}dx\\ &=\frac {1}2\int _{0< s,t}\frac {1}{t+1}\left (\frac {1}{s+t}-\frac {1}{s+t+2}\right )dsdt\\ &=\frac {1}2\int _{0< t}\frac {1}{t+1}\ln \frac {t+2}{t}dt\\ &=\frac {1}2\int _{1}^{\infty }\frac {1}t\ln \frac {t+1}{t-1}dt\\ &=\frac {1}2\int _{0}^{1}\frac {1}t\ln \frac {1+t}{1-t}dt\\ &=\int _{0}^{1}\sum _{0\leq n}\frac {t^{2n}}{2n+1}dt\\ &=\sum _{0\leq n}\frac {1}{(2n+1)^{2}}\\ &=\left (1-\frac {1}{2^2}\right )\zeta (2)\\ &=\frac {\pi ^{2}}8 \end {aligned} $$
逆にDirichlet積分から$\zeta(2)$が求まったとも言えますね。
同様にして、複素積分から以下が出てきます。
$$ \begin {aligned} \int _{0< x< y< z}\frac {\sin x}x\frac {e^{iy}}y\frac {e^{iz}}zdxdydz &=i\int _{0< x< y}\frac {\sinh x}x\frac {e^{-y}}y\frac {e^{-z}}zdxdydz \end {aligned} $$
実部と虚部を比較して
$$ \begin {aligned} S(3)-S(1,1,1)&=0\\ S(1,2)+S(2,1)&=\int _{0< x< y< z}\frac {\sinh x}x\frac {e^{-y}}y\frac {e^{-z}}zdxdydz \end {aligned} $$
シャッフル積から
$$ \begin {aligned} S(3)&=S(1,1,1)=\frac {1}6S(1)^2=\frac {\pi ^{3}}{48} \end {aligned} $$
となります。こんな感じで他にもMSV間の関係式
$$ \begin {aligned} S(4)&=S(2,1,1)\\ S(1,2,1)&=-S(1,1,2) \end {aligned} $$
などが出てきます。

特殊値

色々特殊値が求まってきたので整理しましょう。

Weight 1

Dirichlet積分より、
$$ \begin {aligned} S(1)&=\frac {\pi }2 \end {aligned} $$

Weight 2

先に解説した通り、
$$ \begin {aligned} S(2)&=0\\ S(1,1)&=\frac {\pi ^{2}}8 \end {aligned} $$

Weight 3

先に解説した通り、
$$ \begin {aligned} S(3)&=S(1,1,1)=\frac {\pi ^{3}}{48} \end {aligned} $$
です。またシャッフル積より
$$ \begin {aligned} S(1)S(2)&=2S(1,2)+S(2,1) \end {aligned} $$
$S(2)=0$ですから、
$$ \begin {aligned} S(2,1)=-2S(1,2) \end {aligned} $$
が従い、先の解説と同様にして
$$ \begin {aligned} -S(1,2)&=\int _{0< x< y< z}\frac {\sinh x}x\frac {e^{-y}}y\frac {e^{-z}}zdxdydz\\ &=\frac {1}2\int _{0< s,t,u}\frac {1}{(u+1)(t+u+2)}\left (\frac {1}{s+t+u+1}-\frac {1}{s+t+u+3}\right )dsdtdu\\ S(2,1)&=\int _{0< t,u}\frac {1}{(u+1)(t+u+2)}\ln \frac {t+u+3}{t+u+1}dtdu\\ &=\int _{1< u< t}\frac {1}{u(t+1)}\ln \frac {t+2}{t}dudt\\ &=\int _{1}^{\infty }\frac {\ln t}{t+1}\ln \frac {t+2}tdt\\ &=\int _{0}^{1}\frac {\ln \frac {1}t\ln (1+2t)}{t(1+t)}dt \end {aligned} $$
となり、あとは計算をするだけ...
はあぁぁ~~???なんだこれ???
もしかしたら何かしらの対称性とかがあって計算しやすい形に変形できるのかもしれませんがそんなの知ったこっちゃないのでゴリ押します。
不定積分は一応ポリログで書けるので...すが、代入して整理していける量じゃないです。ズルします。その前にちょっと変形してやります。
$$ \begin {aligned} \int _{0}^{1}\frac {\ln \frac {1}t\ln (1+2t)}{t(1+t)}dt &= \int _{0}^{1}\frac {\ln \frac {1}t\left (\ln \left (1+\frac {t}{1+t}\right )+\ln (1+t)\right )}{t(1+t)}dt\\ &=\underbrace{\int _{0}^{1}\frac {\ln \frac {1}t\ln \left (1+\frac {t}{1+t}\right )}{t(1+t)}dt}_{I_1}+\underbrace{\int _{0}^{1}\frac {\ln \frac {1}t\ln (1+t)}{t(1+t)}dt}_{I_2} \end {aligned} $$
まず$I_1$から計算します。
$$ \begin {aligned} I_1&=\int _{0}^{1}\frac {\ln \frac {1}x\ln \left (1+\frac {x}{1+x}\right )}{\frac {x}{1+x}}\frac {dx}{(1+x)^{2}}\\ &=\int _{0}^{1}\frac {\operatorname{Li} _{2} \left (-\frac {x}{1+x}\right )}{x}dx\\ &= \int _{0}^{1}\frac {1}{1+x}\int _{0}^{1}\frac {\ln \frac {1}t}{1+\frac {xt}{1+x}}dtdx\\ &=\int _{0}^{1}\int _{0}^{1}\frac {\ln \frac {1}t}{1+x(1+t)}dxdt\\ &=-\int _{0}^{1}\frac {\ln t\ln (2+t)}{1+t}dt \end {aligned} $$
ここで、こちらの論文(arXiv:2004.06232)$\text{Theorem 1.1}$$(3.3)$を参照すると
$$ \begin {aligned} -I_1+\frac {3}4\zeta (3)&=\frac {5}{24}\zeta (3)\\ I_1&=\frac {13}{24}\zeta (3) \end {aligned} $$
らしいです。次に$I_2$を計算します。
$$ \begin {aligned} I_2&=-\sum _{0< n}\frac {(-1)^{n}}n\sum _{0\leq m}(-1)^m\int _{0}^{1}x^{n+m-1}\ln \frac {1}xdx\\ &=-\sum _{0< n}\sum _{0\leq m}\frac {(-1)^{n+m}}{n(n+m)^{2}}\\ &=-\sum _{0< n\leq m}\frac {(-1)^{m}}{nm^{2}}\\ &=-\zeta^{\star}(1,\bar {2})\\ &=-\zeta(\bar {3})-\zeta (1,\bar {2})\\ &=\frac {3}4\zeta (3)-\frac {1}8\zeta (3)\\ &=\frac {5}8\zeta (3) \end {aligned} $$
最後のAMZVの計算方法は、こちらの記事で紹介されています。
これでようやく準備が整いました。
$$ \begin {aligned} S(2,1)&=I_1+I_2\\ &=\frac {13}{24}\zeta (3)+\frac {5}8\zeta (3)\\ &=\frac {7}{6}\zeta(3) \end {aligned} $$
さらに、
$$ \begin {aligned} S(1,2)&=-\frac {7}{12}\zeta (3) \end {aligned} $$
です。まとめると、
$$ \begin {aligned} S(3)&=S(1,1,1)=\frac {\pi ^{3}}{48}\\ S(2,1)&=-2S(1,2)=\frac {7}6\zeta (3) \end {aligned} $$
となります。

Weight 4

複素積分から出てきた関係式
$$ \begin {aligned} S(4)&=S(2,1,1)\\ S(1,2,1)&=-S(1,1,2) \end {aligned} $$
と、シャッフル積を組み合わせます。
$$ 0=S(2)^2=4S(1,3)+2S(2,2) $$より
$$S(2,2)=-2S(1,3)$$です。この値の多重ゼータ値による表示の予想は現在ありませんが、周期に属することは分かっています。
追記:この記事を書いた直後にAMZVで書けることが発覚しました!(後で解説しています)
$$ \begin {aligned} S(1)S(3)&=2S(1,3)+S(2,2)+S(3,1)=S(3,1)\\ S(1)^4&=24S(1,1,1,1) \end {aligned} $$より
$$ \begin {aligned} S(3,1)&=4S(1,1,1,1)=\frac {\pi ^{4}}{96}=\frac{15}{16}\zeta(4) \end {aligned} $$
です。
$$ \begin {aligned} 0&=S(2)S(1,1)\\ &=3S(1,1,2)+2S(1,2,1)+S(2,1,1)\\ &=-S(1,1,2)+S(2,1,1) \end {aligned} $$
より
$$ S(4)=S(2,1,1)=S(1,2,1)=-S(1,1,2) $$であり、
$$S(1)S(2,1)=2S(1,2,1)+S(2,1,1)=4S(4)$$なので
$$ S(4)=S(2,1,1)=S(1,2,1)=-S(1,1,2)=\frac {7}{48}\pi \zeta (3) $$
です。残った$S(1,3)$$S(2,2)$を求めていきます。
複素積分から
$$ \begin {aligned} \int _{0< w< x< y< z}\frac {\sin w }w\frac {e^{-ix}}x\frac {e^{iy}}y\frac {e^{iz}}zdwdxdydz &=i\int _{0< w< x< y< z}\frac {\sinh w}w\frac {e^{x}}x\frac {e^{-y}}y\frac {e^{-z}}zdwdxdydz \end {aligned} $$
で、この虚部を見ると
$$ \begin {aligned} S(3,1)+S(2,2)-S(1,3)+S(1,1,1,1)&=\int _{0< w< x< y< z}\frac {\sinh w}w\frac {e^{x}}x\frac {e^{-y}}y\frac {e^{-z}}zdwdxdydz \end {aligned} $$
となります。
左辺を整理すると
$$ \begin {aligned} \frac {75}{64}\zeta (4)-3S(1,3) \end {aligned} $$
右辺を計算すると
$$ \begin {aligned} &\int _{0< w< x< y< z}\frac {\sinh w}w\frac {e^{x}}x\frac {e^{-y}}y\frac {e^{-z}}zdwdxdydz\\ &=\frac {1}2\int _{0< r,s,t,u}\frac {1}{(u+1)(t+u+2)(s+t+u+1)}\left (\frac {1}{r+s+t}-\frac {1}{r+s+t+2}\right )drdsdtdu\\ &=\frac {1}2\int _{0< u< t< s< r}\frac {1}{(u+1)(t+2)(s+1)}\left (\frac {1}r-\frac {1}{r+2}\right )drdsdtdu\\ &=\frac {1}2\int _{1< u< t< s}\frac {1}{u(t+1)s}\ln \frac {s+1}{s-1}dsdtdu\\ &=\frac {1}2\int _{0< s< t< u<1}\frac {1}{ut(1+t)s}\ln \frac {1+s}{1-s}dsdtdu\\ &=\frac {1}2\sum _{0< n}\sum _{0\leq m}\frac {1-(-1)^{n}}n(-1)^m\int _{0< s< t< u<1}\frac {s^{n-1}t^{m-1}}{u}dsdtdu\\ &=\frac {1}2\sum _{0< n}\sum _{0\leq m}\frac {(-1)^{m}-(-1)^{n+m}}{n^2(n+m)^2}\\ &=\frac {1}2\left (\zeta^{\star }(\bar {2},\bar {2})-\zeta ^{\star}(2,\bar {2})\right )\\ &=\frac {1}2\left (\zeta (4)+\zeta (\bar {2},\bar {2})-\zeta (\bar {4})-\zeta (2,\bar {2})\right )\\ &=\frac {1}2\left (\zeta (4)-\frac {3}{16}\zeta (4)+\frac {7}8\zeta (4)-\frac {5}{16}\zeta (4)+2\zeta (1,\bar {3})\right )\\ &=\frac {11}{16}\zeta (4)+\zeta (1,\bar {3}) \end {aligned} $$
なので、
$$ \begin {aligned} S(1,3)&=\frac {31}{192}\zeta (4)-\frac {1}3\zeta (1,\bar {3})\\ S(2,2)&=-\frac {31}{96}\zeta (4)+\frac {2}3\zeta (1,\bar {3}) \end {aligned} $$
まとめると、
$$ \begin {aligned} S(3,1)&=4S(1,1,1,1)=\frac{15}{16}\zeta(4)\\ S(4)&=S(2,1,1)=S(1,2,1)=-S(1,1,2)=\frac{7}{48}\pi \zeta(3)\\ S(2,2)&=-2S(1,3)=-\frac {31}{96}\zeta (4)+\frac {2}3\zeta (1,\bar {3}) \end {aligned} $$
以上です。

参考文献

[1]
Tewodros Amdeberhan, Victor H. Moll, Armin Straub and Christophe Vignat, A triple integral analog of a multiple zeta value, International Journal of Number Theory, 2021, pp.223-237
投稿日:2022614
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。

投稿者

便利
便利
262
50347
引き算が苦手です

コメント

他の人のコメント

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