[论文解读] Fuchsia and master integrals for splitting functions from differential equations in QCD

O. Gituliar, Vitaly Magerya|arXiv (Cornell University)|Jul 4, 2016

Particle physics theoretical and experimental studies参考文献 18被引用 11

一句话总结

本文提出 Fuchsia，一个罗马·李方法的开源 Python 实现，用于将费曼图主积分的微分方程约化为规范形式。该工具通过福克斯化、归一化和因式分解步骤，实现了自动、系统的约化，成功应用于量子色动力学中 NLO 分裂函数的计算，给出了一个 6×6 系统的显式结果，展示了其在多圈计算中的高精度与可扩展性。

ABSTRACT

We report on the recent progress in reducing differential equations for Feynman master integrals to canonical form with the help of a method proposed by Roman Lee. For the first time, we present Fuchsia --- our open-source implementation of the Lee algorithm written in Python using mathematical routines of a free computer algebra system SageMath. We demonstrate Fuchsia by reducing differential equations for NLO contributions to splitting functions in QCD, which contain both loops and legs integrals.

研究动机与目标

• 为解决在量子场论中，罗马·李方法用于将微分方程约化为规范形式的公开可用实现缺乏的问题。
• 自动化约化多圈和相空间费曼积分中产生的复杂微分方程系统，用于 QCD 计算。
• 提供一个可扩展的开源解决方案，基于 SageMath 进行符号计算，支持未来高能物理的高精度计算。
• 通过系统性地约化主积分系统，实现量子色动力学中下一阶（NLO）和下一下一阶（NNLO）分裂函数的计算。

提出的方法

• Fuchsia 实现了李的三步算法：福克斯化以消除非正则奇点，归一化以使本征值为 ε 的整数倍，以及因式分解以实现与 ε 成比例的规范形式。
• 福克斯化步骤使用改进的 Moser 约化和 P-平衡变换，将具有在 {0,1,∞} 处的 Poincaré 秩 {2,1,0} 的系统转换为 Poincaré 秩全为零的 Fuchsian 系统。
• 通过迭代调整奇异点处的本征值，利用广义本征向量导出的投影算子实现归一化，确保残差之和保持为零。
• 因式分解应用常数变换，将系统表示为规范形式 ∂g/∂x = εS(x)g，其中 S(x) 具有有理系数，ε 为维度正则化参数。
• 该方法支持块三角形和稀疏矩阵结构，提升了大规模系统的性能，并通过用户可控的随机种子探索不同的有效规范形式。
• Fuchsia 基于 SageMath 构建，利用其自由的计算机代数系统进行符号运算，并设计为可扩展，适用于多变量和多尺度系统。

实验结果

研究问题

• RQ1 是否可以使用公开可用的开源工具，系统且自动地计算费曼主积分微分方程的规范形式？
• RQ2 李算法的实现对 NLO QCD 分裂函数中产生的复杂非 Fuchsian 系统约化为规范形式的效率如何？
• RQ3 Fuchsia 工具在具有多个奇异点（包括复数奇异点）的大规模系统中的性能和可扩展性如何？
• RQ4 通过算法选择（如随机种子选择）在多大程度上可使规范形式唯一化或受用户影响？
• RQ5 当前符号计算系统（如 SageMath）在处理复杂奇点和多变量系统方面存在哪些局限性，如何克服？

主要发现

• Fuchsia 顺利地将 QCD 中 NLO 分裂函数的 6×6 微分方程系统约化为规范形式，实现了期望的 ε 成比例结构 ∂g/∂x = εS(x)g。
• 福克斯化和归一化后的归一化系统中，本征值与 ε 成比例，矩阵结构表现出对 x 和 1−x 的有理依赖性。
• 通过常数变换矩阵 T(x,ε) 实现了最终的规范形式，该矩阵被显式计算并证明可生成在 x=0 和 x=1 处具有极点的有理系数 S(x) 矩阵。
• 该工具展示了可扩展性，可在标准笔记本电脑上约一小时内完成包含 20 个奇异点和 3×3 耦合块的 74×74 矩阵的约化。
• 尽管存在挑战，Fuchsia 仍提供了一个稳健的高精度计算框架，预计可扩展至 NNLO 和多尺度系统。
• 在 SageMath 处理复多项式因式分解方面观察到局限性，这目前阻碍了具有复奇异点系统的完全自动化。