[论文解读] Automated one-loop calculations with GoSam

G. Cullen, N. Greiner+6 more|arXiv (Cornell University)|Jan 13, 2012

Particle physics theoretical and experimental studies被引用 54

一句话总结

GoSam 是一个用于高能物理中一阶微扰量子场论振幅自动计算的软件框架，通过费曼图生成积分表达式，并利用 d 维单位性或张量约化方法进行约化。该框架支持量子色动力学（QCD）、电弱理论及新物理（BSM）模型中的精确 NLO 计算，通过 Binoth Les Houches 接口与蒙特卡洛程序（如 Sherpa）完全兼容，在多种多粒子过程的计算中与文献基准结果一致。

ABSTRACT

In this talk, the program package GOSAM is presented, which can be used for the automated calculation of one-loop amplitudes for multi-particle processes. The integrands are generated in terms of Feynman diagrams and can be reduced by d-dimensional integrand-level decomposition, or tensor reduction, or a combination of both. Through various examples we show that GOSAM can produce one-loop amplitudes for both QCD and electroweak theory; model files for theories Beyond the Standard Model can be linked as well.

研究动机与目标

• 自动化生成和数值计算多粒子过程在量子场论中的一阶振幅。
• 提供一个灵活可扩展的框架，通过 FeynRules 和 LanHEP 等模型接口支持 QCD、电弱理论及标准模型之外（BSM）的理论。
• 通过 Binoth Les Houches 协议接口，实现与部分子喷注蒙特卡洛程序的无缝集成，支持 NLO+部分子喷注的匹配计算。
• 通过结合 d 维单位性约化（Samurai）与张量约化（golem95C 或 PJFRY），在重建失败时自动切换，确保数值稳定性和高效性。
• 通过隐式与显式构造方法支持振幅中合理项的计算，保留 $\bar{q}^2$ 与 $\tilde{q}^2$ 分解的精度。

提出的方法

• 程序使用 QGRAF 从费曼图生成一阶振幅，随后通过 FORM 进行符号运算，利用 Spinney 处理旋量代数。
• 采用 't Hooft-Veltman 与 DRED 方案，通过 $q^2 = \hat{q}^2 - \mu^2$ 将环动量分解为 4D 与 $(D-4)$ 维部分。
• 将被积函数分离为动量依赖部分与运动学部分，对非 $\hat{q}$、非 $\mu^2$ 项进行简写以提高效率。
• 使用 haggies 为每张图生成优化的 Fortran95 代码，用于数值计算。
• 约化过程采用 d 维 OPP 方法（通过 Samurai）或张量约化（通过 golem95C 或 PJFRY），在重建失败时自动回退。
• 合理项通过约化过程隐式计算，或利用对 $R_2$ 项有贡献的积分的解析知识显式计算。

实验结果

研究问题

• RQ1 自动化框架能否在高数值稳定性与效率下，为 QCD 与电弱理论中的多粒子过程生成一阶振幅？
• RQ2 d 维单位性约化与张量约化相结合，能否有效处理复杂振幅中的数值不稳定性？
• RQ3 GoSam 框架在使用 FeynRules 与 LanHEP 等标准模型接口工具时，对新物理模型的支持程度如何？
• RQ4 GoSam 在 $W^-$+喷注产生与 $pp \to t\bar{t}$ 这类基准过程中的 NLO 计算中，能否准确重现已知结果？
• RQ5 是否可在不预先分离的情况下可靠计算振幅中的合理项？该方法是否提升性能或精度？

主要发现

• GoSam 成功重现了 25 个基准过程的已知结果，包括 $e^+e^- \to t\bar{t}$、$pp \to W^\pm j$、$pp \to t\bar{t}b\bar{b}$ 与 $\gamma\gamma \to \gamma\gamma\gamma\gamma$，与文献值完全一致。
• 通过 Binoth Les Houches 协议与 SHERPA 的接口，实现了 NLO+部分子喷注的匹配计算，$\sqrt{s} = 7$ TeV 下 $W^-$+喷注产生的横动量与伪快度分布与 MCFM 结果完全一致。
• 混合约化策略——以 Samurai 为主约化，golem95C 或 PJFRY 为备用——确保了鲁棒性与稳定性，尤其在高多重性或低秩振幅中表现优异。
• 在被积函数层面采用张量重构，提升了数值稳定性并可降低计算成本，尤其在外部粒子数较多的振幅中效果显著。
• 框架能正确处理具有复质量的不稳定粒子，支持包含共振中间态（如顶夸克或希格斯玻色子）的计算。
• 程序同时支持隐式与显式合理项计算，后者允许对合理部分进行独立评估，提升了在特定应用中的灵活性。