[论文解读] Numerical integration of massive two-loop Mellin-Barnes integrals in Minkowskian regions

Ievgen Dubovyk, J. Gluza+2 more|arXiv (Cornell University)|Jan 1, 2016

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

一句话总结

本论文提出 MBnumerics.m，一种新颖的数值工具，用于在闵可夫斯基运动学区域中高精度地计算大量两圈Mellin-Barnes积分。通过结合轮廓变形、变量旋转和优化积分网格，该方法稳定了振荡性被积函数，实现了对阈值效应和多尺度效应的精确计算，并通过Z→bb顶点修正的解析结果验证了其准确性。

ABSTRACT

Mellin-Barnes (MB) techniques applied to integrals emerging in particle physics perturbative calculations are summarized. New versions of AMBRE packages which construct planar and nonplanar MB representations are shortly discussed. The numerical package MBnumerics.m is presented for the first time which is able to calculate with a high precision multidimensional MB integrals in Minkowskian regions. Examples are given for massive vertex integrals which include threshold effects and several scale parameters.

研究动机与目标

• 解决在闵可夫斯基区域中评估Mellin-Barnes积分时缺乏稳健数值工具的问题，特别是针对多尺度和阈值敏感振幅。
• 克服由两圈费曼图中高维MB积分产生的振荡被积函数和收敛缓慢所导致的数值不稳定性。
• 开发一个完全自动化的高精度数值框架，能够处理复杂物理运动学，包括质量传播子和多尺度。
• 实现对电弱两圈修正中解析结果的可靠交叉验证，特别是Z→bb衰变等过程。
• 提供一个实用且公开可用的工具（MBnumerics.m），其性能优于现有方法，尤其在大多数数值包失效的闵可夫斯基区域。

提出的方法

• 在复平面上应用轮廓变形，通过同时旋转所有MB积分变量来稳定被积函数，避免极点并改善收敛性。
• 通过变量旋转 z → z₀ + (i + θ)t 来抑制振荡行为，其中θ的选择使得指数的实部为负，以实现最佳收敛。
• 采用Cuhre算法（来自CUBA库）进行确定性积分，以在最优轮廓上实现高精度评估，避免蒙特卡洛噪声。
• 将该方法与AMBRE生成的MB表示相结合，用于平面和非平面图，实现从图到数值结果的完全自动化。
• 应用正切映射（tᵢ = 1/tan(−πTᵢ)）以改善奇点和被积函数尾部附近的采样。
• 在欧几里得和闵可夫斯基区域中，将结果与解析解及其他数值工具（如FIESTA、SecDec）进行验证。

实验结果

研究问题

• RQ1 尽管存在振荡行为和收敛缓慢，是否能在闵可夫斯基区域中高精度地数值计算Mellin-Barnes积分？
• RQ2 轮廓变形和变量旋转在稳定多尺度和阈值效应的高维MB积分方面有多有效？
• RQ3 MBnumerics.m是否能在闵可夫斯基区域中优于现有数值工具（如MB(Vegas)、FIESTA、SecDec），特别是在复杂物理振幅中？
• RQ4 该方法在处理具有质量内部线和阈值效应的非平面图方面的能力如何？
• RQ5 轮廓变形与确定性积分（Cuhre）的结合是否足以在多尺度两圈计算中实现8位数字的精度？

主要发现

• MBnumerics.m 实现了对闵可夫斯基区域中四维质量两圈Mellin-Barnes积分的高精度数值评估，结果与解析值匹配至8位数字精度。
• 对于Z→bb顶点修正，MBnumerics.m 以10⁻⁹量级的相对误差重现了解析结果，而MB(Vegas)和SecDec则因大误差或NaN值而失败。
• 该方法成功处理了阈值效应和多尺度运动学，如对具有质量内部线的3D非平面顶点积分常数部分的精确评估所示。
• 通过θ = 0.7的轮廓变形稳定了被积函数渐近行为，将一个收敛缓慢的积分（KC₁）转变为收敛迅速的积分（KC₂），从而实现了可靠的数值积分。
• 在闵可夫斯基区域中，MBnumerics.m 得到复数结果为 −0.778599608324769 − 4.123512600516016i，与解析值 −0.778599608979684 − 4.123512593396311i 匹配至8位数字。
• 即使对于具有一个质量交叉线的非平面图，该方法也表现出鲁棒性，经由PlanarityTest.m工具包验证，并在图6中可视化。