[论文解读] Precise Numerical Evaluation of the Scalar One-Loop Integrals with the Infrared Divergence

Fukuko Yuasa, Elise de Doncker+4 more|ArXiv.org|Sep 6, 2007

Particle physics theoretical and experimental studies参考文献 14被引用 10

一句话总结

本文提出了一种新颖的数值方法，用于计算具有红外发散的标量一 Loop 积分，结合自适应数值积分、在调节参数 $\epsilon$ 上的外推法以及多精度算术。即使在光子质量小至 $10^{-150}$ GeV 时，该方法仍能与解析结果保持十位有效数字的一致性，证明了在高精度量子场论计算中，精度控制对于获得稳定可靠结果至关重要。

ABSTRACT

We present a new approach for obtaining very precise integration results for infrared vertex and box diagrams, where the integration is carried out directly without performing any analytic integration of Feynman parameters. Using an appropriate numerical integration routine with an extrapolation method, together with a multi-precision library, we have obtained integration results which agree with the analytic results to 10 digits even for such a very small photon mass as $10^{-150}$ GeV in the infrared vertex diagram.

研究动机与目标

• 开发一种针对具有红外奇点的标量一 Loop 顶点和盒子图的鲁棒数值方法。
• 解决当虚构光子质量 $\lambda$ 极其小时（如 $10^{-150}$ GeV）仍保持数值精度的挑战。
• 证明在高精度量子场论计算中，精度控制对于获得稳定且准确结果至关重要。
• 通过在广泛 $\lambda$ 值范围内将数值结果与解析解进行比较，验证该方法的有效性。

提出的方法

• 该方法使用数值积分程序（DQAGE）结合自适应求积法，直接计算圈积分，无需进行解析参数积分。
• 在传播子分母中引入调节项 $i\epsilon$ 以避免发散，其中 $\epsilon$ 被视为有限正参数。
• 计算一系列随 $\epsilon_l$ 递减的积分 $I(\epsilon_l)$，随后外推至 $\epsilon \to 0$ 以恢复物理极限。
• 采用多精度算术（双精度、四重精度，最高至 32 精度）控制数值误差，尤其在 $\lambda$ 极小时更为关键。
• 该方法将精度控制与外推法相结合，以稳定结果并检测计算过程中丢失的位数。
• 该方法被应用于一 Loop 顶点图和盒子图，参数根据每张图及 $\lambda$ 值进行调整。

实验结果

研究问题

• RQ1 是否能够通过完全的数值方法，使一 Loop 顶点图和盒子图在具有红外奇点时，与解析结果达到十位有效数字的一致性？
• RQ2 当光子质量 $\lambda$ 降低至 $10^{-30}$ GeV 以下时，所需数值精度如何变化？
• RQ3 精度控制在稳定具有红外发散的圈积分数值结果中起到何种作用？
• RQ4 当标准双精度或四重精度失效时，$\epsilon$ 外推能否可靠地恢复 $\epsilon \to 0$ 的极限？
• RQ5 为保持极高精度而计算极小 $\lambda$ 值时，相关的计算成本如何？

主要发现

• 对于 $\lambda = 10^{-150}$ GeV 的一 Loop 顶点图，实部的数值结果与解析值一致至 11 位有效数字。
• 在 $\lambda = 10^{-150}$ GeV 时，使用 32 精度算术，顶点图的虚部也实现了 10 位有效数字的一致性。
• 对于 $\lambda = 10^{-25}$ GeV 的一 Loop 盒图，数值结果与解析结果一致至 7 位；而在 $\lambda = 10^{-30}$ GeV 时，若缺乏足够的精度控制，一致性下降至 4 位。
• 该方法即使在 $\lambda = 10^{-150}$ GeV 时，仍能与解析结果保持十位有效数字的一致性，证实了在极端情况下高精度算术的必要性。
• 在 Xeon 3.06 GHz 处理器上，顶点图的 $\lambda = 10^{-150}$ GeV 计算耗时达 6.8 天，凸显了高精度计算的计算成本。
• 本研究证实，精度控制并非可选项，而是确保具有红外发散的圈积分结果可靠性的关键，尤其当 $\lambda < 10^{-30}$ GeV 时。