[论文解读] MBnumerics: Numerical integration of Mellin-Barnes integrals in physical regions

Johann Usovitsch, Ievgen Dubovyk+1 more|arXiv (Cornell University)|Jan 1, 2018

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

一句话总结

本论文提出 MBnumerics，一种用于在闵可夫斯基运动学中评估梅林-巴恩斯积分的数值方法——这对Z玻色子共振处的电弱二圈修正至关重要。该方法引入了 contour shift（轮廓偏移）和余切映射，以稳定多项式发散的被积函数，从而实现标准方法失效时费曼积分的精确计算。

ABSTRACT

We introduce techniques to treat numerically Mellin-Barnes integrals in physical regions, which arise in the need of the computation of Feynman integrals for the electroweak two-loop corrections to the pseudo observables at the Z-boson resonance.

研究动机与目标

• 解决在闵可夫斯基运动学中由于渐近行为为多项式而非指数衰减而引起的梅林-巴恩斯积分的数值不稳定性。
• 克服在高能区域中因收敛性问题和振荡行为导致的标准数值积分困难。
• 在存在大质量传播子的情况下，实现对Z玻色子共振处物理可观测量的二圈费曼积分的精确计算。
• 通过梅林-巴恩斯表示降低积分维度，为多尺度积分提供一种鲁棒的替代 sector decomposition（扇区分解）的方法。
• 建立系统化的轮廓偏移与留数修正框架，以提高数值收敛性和精度。

提出的方法

• 在复平面上应用轮廓偏移（zi = xi + iti + ni），以改善梅林-巴恩斯被积函数的渐近行为，减少多项式增长。
• 使用余切映射（t → cot(πt)）变换积分变量，确保边界为有限值，并避免对数映射中常见的奇点。
• 通过对对数 Gamma 函数进行变换：∏Γi → exp(∑logΓi)，提升大 |zi| 时的数值稳定性。
• 在轮廓偏移穿过极点时实施留数修正，增加 (n−1) 重积分，其数值计算更为简便。
• 与梅林-巴恩斯主公式结合：1/(a+b)^ν = ∫ dz/(2πi) a^z b^{-z-ν} Γ(−z)Γ(ν+z)/Γ(ν)，该公式在 |arg a − arg b| < π 时成立。
• 使用费曼参数表示和西曼尼克多项式（U(x), F(x)）在数值计算前推导梅林-巴恩斯形式。

实验结果

研究问题

• RQ1 如何在闵可夫斯基运动学中对具有多项式发散渐近行为的梅林-巴恩斯积分进行稳定化处理，以实现数值计算？
• RQ2 轮廓偏移和变量变换在多尺度二圈积分中能在多大程度上改善收敛性并减少数值误差？
• RQ3 在物理区域积分中，梅林-巴恩斯方法是否能在积分维度和计算效率方面优于扇区分解？
• RQ4 轮廓偏移后留数修正对最终结果的精度和稳定性有何影响？
• RQ5 映射选择（余切 vs. 对数）如何影响数值稳定性和边界行为？

主要发现

• 通过将轮廓偏移量设为 n1 = −2，原始积分值从 0.3923828588857 + 0.7456388536613i 降低至 −0.00974965823202，表明其量级显著减小。
• 在修正被偏移轮廓所包围的三个极点的留数后，修正结果与原始积分高度一致：0.402132517117807 + 0.745638853661318i。
• 对于非平凡的三尺度积分（图4），MBnumerics 在 9000 万次迭代下与扇区分解（SD）在 ε⁻²、ε⁻¹ 和有限阶项上均达到 10⁻⁸ 的一致精度。
• 该积分的梅林-巴恩斯表示最多为三维，而扇区分解需五维，证实了维度上的关键优势。
• 该方法成功计算了包含大质量传播子（MZ, MW, mt）的物理区域积分，验证了其在真实电弱二圈计算中的适用性。
• 余切映射保持了有限的边界极限，避免了数值不稳定性，而对数映射则在边界处引入无穷大。