[論文レビュー] Numerical Computations Concerning the GRH

David J. Platt|arXiv (Cornell University)|May 14, 2013

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

ひとこと要約

本稿では、ディリクレL関数の効率的かつ厳密な計算のための2つの新規アルゴリズムを提示する。これにより、モジュラス $ q \leq 400,000 $ のすべての原始的ディリクレL関数について、指定された高さまで一般化リーマン予想（GRH）の検証が可能になった。この計算は、高速並列計算クラスタと、FFTおよび区間演算を含む高度な数値技法を用いて実施され、290億件を超えるL関数についてGRHが確認された。また、$ q \leq 2,000,000 $ のすべての原始的特徴に対して中心点での非零性が立証された。

ABSTRACT

We describe two new algorithms for the efficient and rigorous computation of Dirichlet L-functions and their use to verify the Generalised Riemann Hypothesis for all such L-functions associated with primitive characters of modulus q<=400,000. For even q, we check to height t_0=max(1e8/q,7.5e7/q+200) and for odd q to height t_0=max(1e8/q,3.75e7/q+200).

研究の動機と目的

• モジュラス $ q \leq 400,000 $ のすべての原始的ディリクレL関数について、一般化リーマン予想（GRH）を厳密に検証すること。
• 従来の結果（$ q \leq 432 $ かつ低い高さまで）を超えて、GRHの計算的検証を拡張すること。
• ディリクレL関数およびその零点を高精度で効率的に計算できる2つの新規数値アルゴリズムの開発と実装。
• すべての原始的特徴について $ q \leq 2,000,000 $ で $ L_{\chi}(1/2) \neq 0 $ であることを確認し、三項ガウスト問題などの解析的整数論の応用を支援すること。
• 区間演算と高精度のオイラー＝マクローリン和公式を用いて、関数方程式における数値的不安定性と符号の曖昧さを解消すること。

提案手法

• 任意の長さの離散フーリエ変換を効率的に計算できるように改造したブルーステインFFTアルゴリズムを採用し、近似関数方程式を用いたL関数の評価を可能にした。
• ポissonの和公式を適用して整数上の和をフーリエ変換上の和に関連づけ、$ L_{\chi}(1/2 + it) $ を離散点で評価するサンプリングベースの手法を構築した。
• ガンマ因子と指数項で正規化された完成形 $ \Lambda_{\chi}(t) $-関数を用いて、実数値性と対称性を確保し、符号の変化検出を容易にした。
• 2段階のアルゴリズムを実装した：$ q \leq 10,000 $ のための「小規模 $ q $」アルゴリズムと、より大きなモジュラスのための「大規模 $ q $」アルゴリズムで、両者とも並列化可能で、メモリおよびキャッシュ効率に最適化された設計であった。
• 区間演算を用い、倍精度および100ビットMPFI精度を組み合わせて、$ \Lambda_{\chi}(1/2 + it) $ の曖昧な符号の変化を解消し、厳密に零点を検出することを保証した。
• 領域内の零点数を推定するためのターリングの方法を適用し、実際の零点数と照合することで、不一致を是正した。この際、適応的サンプリングと高精度再計算を用いた。

実験結果

リサーチクエスチョン

• RQ1 モジュラス $ q \leq 400,000 $ のすべての原始的ディリクレL関数について、一般化リーマン予想（GRH）を厳密に検証できるか。また、臨界帯内のどの高さまで検証可能か。
• RQ2 大規模な $ q $ に対して、高精度で効率的かつ正確にディリクレL関数およびその零点を計算できる数値的アルゴリズムは何か。
• RQ3 区間演算を用いて厳密な境界を保証することで、関数方程式における数値的不安定性と符号の曖昧さをどのように解消できるか。
• RQ4 すべての原始的ディリクレL関数について $ q \leq 2,000,000 $ で中心値 $ L_{\chi}(1/2) \neq 0 $ であるか。計算的に確認可能か。
• RQ5 現代の計算インfraストラクチャとアルゴリズムの進歩により、従来の研究をはるかに超えてGRH検証の範囲を拡大できるか。

主な発見

• モジュラス $ q \leq 400,000 $ のすべての原始的ディリクレL関数についてGRHが成立し、偶数 $ q $ については $ t_0 = \max\left(\frac{10^8}{q}, \frac{7.5 \cdot 10^7}{q} + 200\right) $、奇数 $ q $ については $ t_0 = \max\left(\frac{10^8}{q}, \frac{3.75 \cdot 10^7}{q} + 200\right) $ の高さまで検証された。
• 計算により、$ q \leq 2,000,000 $ のすべての原始的特徴（合計739,151,526,102個）について $ L_{\chi}(1/2) $ が非零であることが確認された。そのうち20件のみが100ビットMPFI精度を要した。
• すべての $ q \leq 400,000 $ に対して、臨界線上に約 $ 3.8 \cdot 10^{13} $ 個の零点が特定され、非臨界線上には零点は検出されなかった。
• 小規模および大規模 $ q $ に最適化された2つの新規アルゴリズムにより、GRHの検証範囲が従来の厳密計算より6桁も広がった。
• 区間演算と高精度のオイラー＝マクローリン和公式により、$ \Lambda_{\chi}(1/2 + it) $ の符号の曖昧さが厳密に解消され、数値的不安定性による零点の漏れは発生しなかった。
• 全計算は、英国およびフランスのHPCクラスタを用いて約40万コアアワーの計算資源を消費した。両アルゴリズムともマルチコアシステム上で効率的に並列化された。