[论文解读] Extrapolation Algorithms for Infrared Divergent Integrals

Elise de Doncker, J. Fujimoto+6 more|arXiv (Cornell University)|Oct 17, 2011

Numerical methods in inverse problems参考文献 13被引用 5

一句话总结

本文提出一种双外推算法，结合四重精度算术，用于计算量子场论中红外发散一环顶点积分的主导系数。通过对接收正则化参数 ε 和 δ 逐渐减小的积分序列进行非线性外推（采用 ε-算法），该方法在积分区域内部和边界处存在奇异性的情况下，仍能达到接近机器精度的准确性。

ABSTRACT

This paper describes applications of extrapolation for the computation of coefficients in an expansion of infrared divergent integrals. An extrapolation procedure is performed with respect to a parameter introduced by dimensional regularization. While this treats typical IR singularities at the boundaries of the integration domain, special care needs to be taken in cases where the integrand is singular in the interior of the domain as well as on the boundaries. A double extrapolation is devised for a class of massless vertex integrals. Quadruple precision results are presented, demonstrating high accuracy. The computations are supported by the use of general adaptive integration programs from the QUADPACK package, in iterated integrations with highly singular integrand functions.

研究动机与目标

• 解决标准积分方法因不可积奇异性而失效时，计算红外发散圈积分所面临的数值挑战。
• 克服在质量less顶点图中，奇异性同时出现在区域边界和内部时带来的困难。
• 开发一种稳健的数值框架，通过外推技术计算此类积分洛朗展开中系数的数值结果。
• 利用四重精度算术与Quadpack中的自适应积分例程，展示高精度结果。
• 实现在无需事先掌握奇异性结构解析知识的前提下，实现更自动化且可靠的复杂圈积分数值计算。

提出的方法

• 采用维数正则化，引入小参数 ε 来调节红外发散，将积分展开为 ε 的洛朗级数。
• 通过Quadpack中的Dqags/Dqag进行迭代数值积分，计算 ε 值逐渐减小的积分近似值。
• 采用非线性外推技术（使用 ε-算法）加速 ε → 0 时积分近似序列的收敛。
• 对于存在内部奇异性的积分，引入对 δ（分母中的小虚数项 iδ）的第二重外推，以处理不可积奇异性。
• 使用几何序列（如Bulirsch序列）作为 ε 和 δ 的取值序列，以提高外推过程的稳定性和收敛性。
• 由于Quadpack不原生支持复函数，需分别计算复积分的实部与虚部，并对每一部分分别应用外推。

实验结果

研究问题

• RQ1 当奇异性同时出现在积分区域的边界和内部时，如何稳定红外发散顶点积分的数值积分？
• RQ2 外推技术在提高此类积分洛朗展开中系数提取精度方面，能达到何种程度？
• RQ3 双外推（在 ε 和 δ 上）在实现复杂圈积分接近机器精度结果方面表现如何？
• RQ4 不同序列（几何序列与调和序列）对双外推过程的稳定性和收敛性有何影响？
• RQ5 四重精度算术结合自适应积分与外推技术，能否为高精度量子场论计算提供可靠结果？

主要发现

• 双外推方法成功计算了红外发散积分的系数，相对误差低至 5.84×10⁻¹⁶，接近四重精度算术的极限。
• 对于一个外线与两个内线粒子的情况，系数 $ C_{-1} $ 的外推虚部收敛至 1.421026277606×10⁻⁴，误差估计为 5.84×10⁻¹⁶。
• 在八次外推步骤后，$ C_0 $ 的实部计算为 5.24199463865×10⁻⁴，误差估计为 5.64×10⁻¹⁴。
• 虚部收敛速度优于实部，表明奇异性本质导致了不同的收敛行为。
• 使用固定长度（18）的外推序列与 ε-算法可获得稳定结果，但建议未来采用动态终止准则以进一步优化。
• 该方法表明，即使对于具有复杂奇异性的积分，高精度结果仍可实现，支持了圈积分自动化数值计算的可行性。