[论文解读] Multi-gluon one-loop amplitudes numerically

Achilleas Lazopoulos|ArXiv.org|Dec 16, 2008

Particle physics theoretical and experimental studies参考文献 22被引用 24

一句话总结

本文提出了一种C++实现，用于在NLO QCD中数值计算多胶子振幅的一圈单位性切割算法，采用五维中带有胶子和标量的树幅值色序振幅作为构建模块。该方法能够准确评估虚振幅的有理部分和可切割构造部分，通过四重精度回退确保数值稳定性，对于最多22个外部胶子的结果与已发表结果完全一致。

ABSTRACT

A c++ implementation of the D_s-dimensional unitarity cut algorithm for the numerical evaluation of the virtual contribution to NLO QCD amplitudes is presented. The current version includes an arbitrary number of external gluons with gluonic propagators in the loop. The building blocks are tree level color-ordered amplitudes with gluons and with gluons and two scalars in five dimensions. Numerical stability issues are addressed and agreement has been reached with the results published in the literature.

研究动机与目标

• 开发一种通用的、数值稳定的C++工具，用于计算NLO QCD中任意外部胶子数的一圈虚振幅。
• 实现$D_s$维单位性切割方法（EGKM算法），以高效且准确地评估振幅的有理部分和可切割构造部分。
• 通过集成四重精度算术作为回退机制，确保高多重性过程中的数值鲁棒性。
• 在相同框架内支持未来扩展至费米子和规范玻色子，支持LHC中具有物理意义的过程。

提出的方法

• 该算法在被积函数层面使用OPP约化方法，将一圈振幅分解为标量盒、三角形、泡和 tadpole 积分。
• 虚振幅表示为$A^{D_s} = A^0 + A^1 \cdot D_s$，其中$D_s$为内部粒子的维度，通过在两个$D_s$值上进行数值计算，提取$\epsilon$依赖部分和有理部分。
• 使用Berends-Giele递推关系计算五维中带有胶子和两个标量的树幅值色序振幅。
• 该方法在$D_s$维中使用复单位性切割来重构振幅，通过使用QD库的四重精度版本代码保持数值稳定性。
• 该算法应用于具有任意外部手性及相空间点的N胶子振幅，性能通过在$10^5$个相空间点上的平均CPU时间进行衡量。
• 通过与特定相空间点的已发表结果进行对比，确认在数值精度范围内一致。

实验结果

研究问题

• RQ1 如何在NLO QCD中对多胶子过程的完整一圈虚振幅（包括有理部分）进行数值计算？
• RQ2 在具有大量外部粒子的高多重性振幅中，哪些数值策略能有效保持稳定性？
• RQ3 在通用C++框架中，$D_s$维单位性切割方法在多大程度上可以实现高效且稳健的实现？
• RQ4 该算法的性能如何随外部胶子数量变化而变化？对于具有物理意义的过程，需要哪些计算资源？

主要发现

• C++实现成功计算了最多22个外部胶子的一圈完整振幅，典型CPU时间从4个胶子的3.4毫秒到22个胶子的1.85秒不等，运行于Intel Xeon X5450处理器上。
• 对于$N=6$个胶子，当开关精度设为$10^{-4}$时，约5%的相空间点需要四重精度算术，表明性能开销可控。
• $N=6$时双极点和单极点的相对误差呈对数分布，随着$N$增加，误差峰值向更大值移动，但有限部分在大多数情况下保持稳定。
• 该代码与文献中已发表结果完全一致，包括在$D_s = 5$和$6$时的验证，确认了实现的正确性。
• 该方法正确恢复了振幅$\epsilon$展开中的$1/\epsilon^2$和$1/\epsilon$极点，以及有理有限部分。
• 该框架可扩展至费米子和规范玻色子，如先前工作所示，现已建立实现完整物理过程的路径。