[論文レビュー] Effective method of computing Li's coefficients and their properties

K. Maślanka|ArXiv.org|Feb 11, 2004

Analytic Number Theory Research参考文献 5被引用数 12

ひとこと要約

この論文は、リーマン予想（RH）の中心的役割を果たすリウヴィル係数 $\lambda_n$ を効果的に数値計算するための手法を提示する。$\lambda_n$ を厳密に増加するトレンドと微小な振動に分解することで、3,300個を超える係数を計算し、安定したトレンドに従属する微小な振動が存在することが明らかになった。これは、有限サンプリングによる決定的検証が不可能であるにもかかわらず、リーマン予想の数値的裏付けを提供する。

ABSTRACT

In this paper we present an effective method for computing certain real coefficients $λ_{n}$ which appear in a criterion for the Riemann hypothesis proved by Xian-Jin Li. With the use of this method a sequence of over three-thousand $λ_{n}$'s has been calculated. This sequence reveals a peculiar and unexpected behavior: it can be split into a strictly growing extit{trend} and some tiny extit{oscillations} superimposed on this trend.

研究の動機と目的

• リーマン予想（RH）のためのリウヴィル係数 $\lambda_n$ を実用的かつ高精度に計算するための手法を開発すること。
• 既存の結果を超えて $n \approx 3300$ まで $\lambda_n$ の数値的挙動を調査すること。
• $\lambda_n$ を支配的トレンドと残差振動に分解し、それぞれの寄与を明確に分析可能にする。
• 今後の研究のための公開データセットとして、$\lambda_n$ および $\eta_n$ の係数を提供すること。

提案手法

• 方法は、$\lambda_n$ をトレンド成分 $\tilde{\lambda}_n$ と振動成分 $\overset{-}{\lambda}_n$ に分解し、$\xi(s)$ の対数形式の $n$ 階微分から導出される。
• トレンド $\overset{-}{\lambda}_n$ は、$s=1$ における $\ln\left(\pi^{-s/2}\Gamma(1+s/2)\right)$ の $n$ 階微分から計算され、解析的に取り扱えるため厳密に増加する。
• 振動成分 $\overset{\sim}{\lambda}_n$ は、$s=1$ における $\ln\left((s-1)\zeta(s)\right)$ の $n$ 階微分から得られ、トレンドからの逸脱を捉える。
• Stieltjes定数 $\gamma_n$ は、Mathematica の任意精度算術を用いて800桁の高精度で計算され、$\eta_n$ の正確な評価を可能にする。
• $\eta_n$ は、$s=1$ の周りの $\zeta(s)$ のローラン級数から導かれた再帰関係式を用いて再帰的に計算され、$\lambda_n$ の効率的計算を可能にする。
• 数値的妥当性は、Mathematica と Maple の両方を用いて検証され、高精度算術における誤差蓄積にもかかわらず、正確性を保証するための相互検証が実施された。

実験結果

リサーチクエスチョン

• RQ1 既存の文献を超えて、大きな $n$ についてリウヴィル係数 $\lambda_n$ を効率的かつ高精度に計算できるか？
• RQ2 数列 $\lambda_n$ は、滑らかなトレンドと微小な振動に構造的に分解可能か？
• RQ3 $\lambda_n$ の振動は有界か、それともトレンドより遅く成長するか？ これはリーマン予想に何を意味するか？
• RQ4 高精度な $\lambda_n$ の数値的計算は、リーマン予想の支持または反証にさらなる証拠を提供できるか？
• RQ5 トレンド成分および振動成分の $\lambda_n$ の漸近的挙動はいかなるものか？

主な発見

• 既存の結果を著しく超える範囲まで、3,300個を超える $\lambda_n$ の値が高精度で計算された。
• 係数 $\lambda_n$ は、$a(1 + n\ln n) + cn$ でよく近似できる厳密に増加するトレンド $\overset{-}{\lambda}_n$ を示し、$a \approx 0.5$ および $c \approx -1.130330701$ となる。これは正確な値 $c = \frac{1}{2}(\gamma - 1 - \ln 2\pi)$ と一致する。
• 振動成分 $\overset{\sim}{\lambda}_n$ の大きさは極めて小さく、トレンドより遅く成長するため、有界または支配的でない場合、リーマン予想が支持される。
• トレンド成分 $\overset{-}{\lambda}_n$ は厳密に増加であることが証明されており、その漸近展開にはベルヌーイ数が関与するが、形式的には発散する級数である。
• 誤差蓄積のため、$n$ が増加するにつれて $\lambda_n$ および $\eta_n$ の有効数字の桁数が減少し、$n \approx 3300$ を超えると精度が制限される。
• 数値的データは、振動がトレンドに比べて無視できるほど小さく、$\lambda_n \geq 0$ の違反が観測されなかったことから、リーマン予想を強く支持する。