[论文解读] New implementation of the sector decomposition on FORM

Takahiro Ueda, J. Fujimoto|arXiv (Cornell University)|Feb 16, 2009

Particle physics theoretical and experimental studies参考文献 2被引用 19

一句话总结

本论文提出了一种在符号计算系统 FORM 上实现的新型扇形分解算法，用于处理多圈费曼积分中的发散问题。通过利用 FORM 高效处理极大规模表达式的能力，该方法通过迭代扇形分裂分离紫外和红外发散，从而实现对复杂圈积分的数值计算。其主要贡献是一个稳健且可扩展的框架，能够对多圈图实现高精度结果，已在基准案例（包括无质量壳上盒图和三盒图）中得到验证。

ABSTRACT

Nowadays the sector decomposition technique, which can isolate divergences from parametric representations of integrals, becomes a quite useful tool for numerical evaluations of the Feynman loop integrals. It is used to verify the analytical results of multi-loop integrals in the Euclidean region, or in some cases practically used in the physical region by combining with other methods handling the threshold. In an intermediate stage of the sector decomposition for the multi-loop integrals, one often has to handle enormously large expressions containing tons of terms. The symbolic manipulation system FORM is originally designed to treat such huge expressions and has a strong advantage for it. In this talk, the implementation of the sector decomposition algorithm on FORM is discussed. A number of concrete examples including cases of multi-loop diagrams are also shown.

研究动机与目标

• 开发一种可扩展的方法，用于数值计算具有复杂发散性的多圈费曼积分。
• 克服标准计算机代数系统在处理扇形分解产生的大规模表达式时的内存限制。
• 将扇形分解与数值外推方法结合，以处理物理阈值奇点。
• 提供从积分输入到 ε 展开中拉伦特系数的可靠、自动化处理流程。

提出的方法

• 扇形分解算法在 FORM 上实现，FORM 是一种专为处理超过主内存规模的表达式而优化的系统，通过基于磁盘的计算实现极小性能损失。
• 该方法通过迭代分割积分区域，将重叠奇点（如紫外和红外奇点）解耦，将参数积分转化为具有因子化奇点的扇形之和。
• 扇形分解后，对 ε 的洛朗展开通过符号计算完成，得到系数 C₂、C₁、C₀ 作为多维积分。
• 使用自适应积分器（如 DCUHRE、DQAGE）对分解后的被积函数进行数值积分，通过相对和绝对容差控制收敛性。
• 对于物理阈值，通过计算有限 iδ 下的积分并外推至 δ→0 来应用数值外推方法，避免路径形变。
• 整个工作流实现自动化：输入积分描述 → FORM 生成被积函数的 Fortran 代码 → 数值积分输出带误差估计的系数。

实验结果

研究问题

• RQ1 能否在 FORM 上高效实现扇形分解，以处理多圈积分中表达式规模的急剧增长？
• RQ2 当处理具有 O(10⁴) 子扇形的扇形分解积分时，FORM 的性能与其它系统相比如何？
• RQ3 数值外推方法在扇形分解后能否有效处理物理阈值奇点？
• RQ4 对于双盒或三盒等多圈图，能够实现何种精度和计算效率？

主要发现

• 该实现成功计算了平面无质量壳上双盒积分至 O(ε⁰) 阶，C₀ 在 134.6 秒内完成，与解析结果的相对误差在 3×10⁻⁴ 以内。
• 对于 s=123、t=−200 的一阶圈壳上盒图，该方法对 C₂、C₁、C₀ 的计算结果与解析值相对误差低于 10⁻¹⁰，验证了数值稳定性。
• 非平面双盒图和三圈传播子已验证至 O(ε²) 阶，证明了该方法在更高圈数和更复杂图中的可扩展性。
• 数值外推方法能有效处理物理阈值，即使在具有复极点的闵可夫斯基区域，结果也与解析值一致。
• 使用 FORM 可处理具有 O(10⁴) 子扇形的表达式，并完成后续的 O(ε⁰) 展开，而不会发生内存溢出，这与标准 CAS 系统形成鲜明对比。
• 非壳上盒图的 C₀ 计算耗时 1105.9 秒，突显了高阶项的计算成本，但实现了低于 10⁻¹⁰ 的精度并具备稳健的误差控制。