[论文解读] New developments in PJFry

J. Fleischer, T. Riemann+1 more|DESY (CERN, DESY, Fermilab, IHEP, and SLAC)|Oct 15, 2012

Black Holes and Theoretical Physics参考文献 29被引用 6

一句话总结

本论文提出 PJFry，一个用于量子场论中数值稳定的一圈张量积分计算的开源 C++ 软件包。它利用维度递推关系与级数展开，避免了逆格拉姆行列式的问题，从而能够稳健地计算最高秩为 6 的五点与六点积分，在高能物理振幅及自动化振幅计算中具有应用价值。

ABSTRACT

We report on recent progress in numerical evaluation of one loop tensor integrals. A public C++ package PJFry implementing algorithms from arXiv:1009.4436 and its extension to hexagons up to rank 6 are presented.

研究动机与目标

• 解决多外腿过程中由于逆格拉姆行列式导致的一圈张量积分数值不稳定性问题。
• 开发一个稳健、公开可用的工具，用于四点函数以外的高阶张量积分。
• 将张量积分约化方法扩展至更高点函数（最高至六点），支持任意质量与任意秩。
• 通过渐近展开在格拉姆行列式较小时保持数值稳定性，实现张量形式因子的精确计算。
• 提供高性能、可重用的框架，用于自动化振幅计算与圈振幅库中。

提出的方法

• 该软件包利用维度递推关系，将张量积分表示为分母幂次偏移的高维标量积分。
• 使用递推关系式 (6) 和 (7)，将张量形式因子约化为高维空间中的标量积分，避免了格拉姆行列式的直接求逆。
• 在格拉姆行列式较小的区域，应用由递推关系式 (10) 导出的级数展开，以保持数值稳定性。
• 通过组合系数与幂次提升/降低算符，将张量形式因子映射为广义标量积分。
• 采用缓存机制，重用低阶点张量积分以加速计算并减少冗余计算。
• 该实现支持多种接口（C、FORTRAN、C++、Mathematica、GoSam），并通过 qd 库支持扩展精度，适用于复质量情况。

实验结果

研究问题

• RQ1 如何在一元格拉姆行列式较小或趋近于零时，稳定地计算一圈张量积分？
• RQ2 对于五点与六点函数，最高效且数值鲁棒的高阶张量积分约化方法是什么？
• RQ3 能否将维度递推与级数展开结合，以在格拉姆行列式较小时保持精度，同时避免引入逆行列式？
• RQ4 如何设计一个通用、高性能的库，以支持任意质量配置下最高秩为 6 的张量积分？
• RQ5 在张量积分约化中，缓存与重用中间结果的最优策略是什么，以提升性能？

主要发现

• PJFry 能够在格拉姆行列式较小的区域中，对五点与六点函数的张量积分实现最高秩为 6 的完全数值稳定计算。
• 该软件包通过基于递推的约化与渐近级数展开避免了逆格拉姆行列式，确保了在所有相空间区域中的鲁棒性。
• 通过缓存机制重用低阶点张量积分，实现高性能计算，且无需额外开销。
• 在包含六个外线粒子与一条内部质量线的测试相空间点上，标量六点积分 $F_0$ 在 $\epsilon=0$ 时的计算结果为 $-30.3579 - 205.8213i$，且 $\epsilon=0,1,2$ 时结果一致。
• 向量六点积分分量通过同一框架准确计算，其结果在复平面上表现出预期行为。
• Mathematica 接口支持 LoopTools 中的标准约定，可直接用于符号振幅计算。