[论文解读] Numerical Computations Concerning the GRH

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

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

一句话总结

本文提出了两种新颖的算法，用于高效且严格地计算狄利克雷L-函数，实现了对模数 $ 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）。
• 将GRH的计算验证范围扩展至先前结果（此前仅达到 $ q \leq 432 $ 且高度较低）的极限之外。
• 开发并实现两种新型数值算法，能够以高精度高效计算狄利克雷L-函数及其零点。
• 确认所有模数 $ q \leq 2,000,000 $ 的本原特征函数满足 $ L_{\chi}(1/2) \neq 0 $，支持解析数论中的应用（如三元哥德巴赫猜想）。
• 通过区间算术和高精度欧拉-麦克拉林求和，解决函数方程中数值不稳定性与符号模糊性问题。

提出的方法

• 采用改进的布莱斯汀FFT算法计算任意长度的离散傅里叶变换，通过近似函数方程实现L-函数的高效评估。
• 应用泊松求和公式将整数上的求和与傅里叶变换上的求和关联，从而构建基于采样的方法，用于在离散点上评估 $ L_{\chi}(1/2 + it) $。
• 使用归一化后的完整 $ \Lambda_{\chi}(t) $-函数（经伽马因子和指数项调整），确保其为实值并具有对称性，便于检测符号变化。
• 实现两阶段算法：针对 $ q \leq 10,000 $ 的“小 $ q $”算法和针对更高模数的“大 $ q $”算法，两者均可并行化，并针对内存与缓存效率进行了优化。
• 采用双精度和100位MPFI精度的区间算术，解决 $ \Lambda_{\chi}(1/2 + it) $ 中符号变化的模糊性，确保零点的严格检测。
• 应用图灵方法估算区域内的零点数量，并通过实际零点计数进行交叉验证，利用自适应采样和高精度重新计算修正差异。

实验结果

研究问题

• RQ1 能否严格验证所有模数 $ q \leq 400,000 $ 的本原狄利克雷L-函数的GRH？其验证高度在临界带内达到多高？
• RQ2 哪些数值算法能够高效且准确地在高精度下计算狄利克雷L-函数及其零点，特别是针对大 $ q $ 的情况？
• RQ3 如何利用区间算术和严格界值，解决函数方程中出现的数值不稳定性与符号模糊性？
• RQ4 所有模数 $ q \leq 2,000,000 $ 的本原狄利克雷L-函数是否都满足中心值 $ L_{\chi}(1/2) \neq 0 $？该结论能否通过计算验证？
• RQ5 现代计算基础设施与算法进步在多大程度上扩展了GRH验证的范围，超越以往工作？

主要发现

• 广义黎曼猜想对所有模数 $ q \leq 400,000 $ 的本原狄利克雷L-函数成立，验证高度为：当 $ 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) $。
• 计算确认了所有 $ 739,151,526,102 $ 个模数 $ q \leq 2,000,000 $ 的本原特征函数满足 $ L_{\chi}(1/2) \neq 0 $，仅有20个案例需要100位MPFI精度来解决。
• 在所有 $ q \leq 400,000 $ 的情况下，共隔离出约 $ 3.8 \cdot 10^{13} $ 个位于临界线上的零点，未发现任何非临界线上的零点。
• 两种新算法（分别针对小 $ q $ 和大 $ q $）使GRH验证的范围比以往严格计算扩大了六个数量级。
• 区间算术与高精度欧拉-麦克拉林求和成功解决了 $ \Lambda_{\chi}(1/2 + it) $ 中的符号模糊性，确保无零点因数值不稳定性而被遗漏。
• 整个计算过程在英国和法国的高性能计算集群上消耗了约40万核心小时，两种算法在多核系统上均实现了高效的并行化。