NNLOJET: a parton-level event generator for jet cross sections at NNLO QCD accuracy

本文介绍了开源 NNLOJET 部分子级事件生成器，该生成器实现了天线减法方法，用于计算电子 - 正电子、轻子 - 强子以及强子 - 强子碰撞中达到 NNLO QCD 精度的喷注截面及相关可观测量。

要点

以下是该论文的通俗解释，通过类比进行翻译。

宏观图景：打造更精准的“水晶球”

想象你是一位物理学家，试图预测两个微小粒子以接近光速相互撞击时会发生什么。你拥有一套规则（标准模型），可以告诉你这些粒子“应该”如何行为。但自然界是混乱的。当粒子碰撞时，它们不仅仅是弹开，而是像从大炮中爆炸出的彩纸一样，喷涌出一场混乱的新粒子阵雨。

为了做出真正准确的预测，你不能只猜测主要的爆炸。你必须计算出因爆炸而产生的那些微小、几乎不可见的涟漪和次级喷流。在粒子物理学中，这些计算被称为NNLO（次次领头阶）。它们极其复杂，就像试图预测飓风中每一粒沙子的确切路径。

本文介绍了一个名为NNLOJET的新开源计算机程序，旨在为“喷注”（粒子喷流）进行这些超精密计算。这就像给科学家提供了一台高清水晶球，能够看清旧模型遗漏的细节。

核心问题：“无限”的混乱

当物理学家试图计算这些碰撞时，他们会遇到一个名为“红外奇点”的数学噩梦。

• 类比：想象你试图计算房间里的人数，但每当有人低语，房间就会变得无限嘈杂。如果你试图累加这些噪音，数字就会趋向无穷大。
• 现实：在粒子物理学中，当粒子移动得非常慢或靠得非常近时，数学计算会爆炸至无穷大。既然现实世界不存在无限能量，这些无穷大仅仅是计算方法中的误差。

为了解决这个问题，作者使用了一种称为天线相减（Antenna Subtraction）的方法。

• 类比：把它想象成降噪耳机。耳机监听恼人的背景噪音（数学上的无穷大），并生成一个“负”声波来完美抵消它。
• 工作原理：该程序分别计算那些混乱的、无穷大的部分，然后将它们减去，只留下干净的、有限的、符合现实世界的答案。

什么是 NNLOJET？

NNLOJET是运行这些计算的“引擎”。在这篇论文之前，许多此类计算就像锁在单一厨师厨房里的秘密食谱。如果你想使用它们，要么得乞求那位厨师，要么得从头开始建造自己的厨房。

NNLOJET 通过开源改变了游戏规则。

• 类比：这就像将完整的食谱书和厨房工具向公众发布。任何拥有电脑的人都可以下载它、烹饪这顿饭，并检查食材。
• “喷注”部分：该程序专门处理“喷注”。在粒子物理学中，当一个夸克或胶子（物质的基本构建块）被射出时，它不会保持孤立。它会瞬间变成一股粒子喷流。我们将这种喷流称为“喷注”。NNLOJET 能够精确预测这些喷流的大小、速度和方向。

如何使用它？（运行卡）

你不需要是编程巫师就能使用 NNLOJET。论文解释说，你通过运行卡（Runcard）来控制程序。

• 类比：把运行卡想象成飞行计划或食谱卡。你不需要知道飞机引擎是如何工作的；你只需要告诉飞行员去哪里。
• 你写什么：你告诉程序：
  • 地点：我们是撞击质子（如在大强子对撞机中）还是电子？
  • 目标：我们要寻找 Z 玻色子、希格斯玻色子，还是仅仅是一股喷注？
  • 规则：喷注应该有多大？多快？
  • 输出：你希望在最后看到什么图表？

论文提供了关于如何编写这份“飞行计划”的详细手册（第 5 和第 6 节），涵盖了从碰撞类型到具体数学设置的所有内容。

工作流程：装配线

计算这些碰撞是如此繁重，以至于一台计算机无法独自完成。那将需要数年时间。

• 类比：想象你需要绘制一幅巨大的壁画。一个人需要花费一生的时间。相反，你雇佣了一支由 100 名画家组成的团队。
• 工具：论文描述了一个工作流脚本（nnlojet-run），它充当工头的角色。
  1. 预热：工头派出一小部分画家去测试墙壁，找出最佳的绘画方式（这是"VEGAS 网格”的自适应过程）。
  2. 生产：一旦计划确定，工头就会派出数百名画家同时开展实际工作。
  3. 收尾：当所有人完成后，工头收集所有部分，检查错误（异常值），并将它们拼接成一幅完美的画面。

它能计算什么？

论文列出了当前食谱（1.0.0 版）中可用的具体“食谱”：

• 电子 - 正电子碰撞：类似于旧的 LEP 对撞机。
• 质子 - 质子碰撞：类似于当前的大强子对撞机（LHC）。
• 特定事件：它可以预测以下粒子的产生：
  • Z 玻色子（光子的重表亲）。
  • W 玻色子（传递弱力的粒子）。
  • 希格斯玻色子（赋予其他粒子质量的粒子）。
  • 光子（光）和喷注。
  • 上述所有组合（例如，希格斯玻色子加一个喷注）。

这为什么重要？

论文总结道，通过使该代码开源且易于使用，它允许科学家将他们的理论预测直接与真实的实验数据进行比较。

• 类比：如果你在建造一辆模型车，你需要知道它是否真的能开。NNLOJET 提供了尽可能准确的“风洞”测试。如果真实数据与 NNLOJET 的预测相匹配，那么我们对宇宙的理解就是正确的。如果它们不匹配，这可能意味着我们发现了一些新东西！

简而言之：这篇论文发布了一个强大、免费且用户友好的工具，帮助物理学家以极高的精度计算粒子碰撞的混乱细节，确保我们的宇宙地图尽可能准确。

技术摘要

NNLOJET 技术摘要：一个用于 NNLO QCD 精度喷注截面的部分子级事件生成器

问题陈述
对标准模型（SM）的精密检验严重依赖于将强子末态（喷注）的实验测量值与理论预测进行比较。这些预测需要计算次次领头阶（NNLO）量子色动力学（QCD）修正，以处理大强耦合常数（$\alpha_s$）以及源于软辐射和共线辐射的复杂红外（IR）奇点结构。虽然次领头阶（NLO）计算已基本实现自动化，但由于需要抵消双虚（VV）、混合实虚（RV）和双实（RR）贡献之间的红外发散，NNLO 计算仍然充满挑战。现有的 NNLO 代码通常是针对特定过程的，未广泛作为开源代码提供，或者需要对每个可观测量进行大量手动调整。因此，亟需一种灵活、开源的部分子级事件生成器，能够在各种对撞机环境（电子 - 正电子、轻子 - 强子以及强子 - 强子）中计算 NNLO 精度的喷注截面。

方法论：天线减法
NNLOJET 代码实现了天线减法（antenna subtraction）方法来处理红外奇点。该方法基于 QCD 实辐射振幅在红外极限下的普适因子化，推导出减法项。

• 核心机制：该方法将未分辨的实辐射构型识别为硬辐射体部分子对。它们之间的辐射由“天线函数”描述。
• 减法项：未积分的减法项由天线函数与较低部分子多重性的约化矩阵元的乘积构建。这些项在天线相空间上积分，生成积分减法项（$d\hat{\sigma}_T$ 和 $d\hat{\sigma}_U$），从而显式地抵消来自虚修正的 $\epsilon$ 极点。
• 实现：代码利用相空间映射，将天线函数的运动学与约化矩阵元因子化。该实现涵盖了末态 - 末态、初态 - 末态以及初态 - 初态辐射体构型。
• 计算框架：NNLOJET 是一个主要用现代 Fortran 编写，并包含 C++ 和 Python 模块的部分子级事件生成器。它使用自适应蒙特卡洛例程 VEGAS 进行相空间积分。代码包含用于通过 HPLOG、TDHPL 和 CHAPLIN 评估特殊函数（调和多对数）的库，并与 LHAPDF 接口以获取部分子分布函数（PDF）。

主要贡献与特性
本文介绍了 NNLOJET 代码的开源版本（v1.0.0），该版本根据 GNU 通用公共许可证 v3.0 发布。主要技术贡献包括：

1. 广泛的过程覆盖：该代码支持以下过程的 NNLO 喷注产生计算：

   • 电子 - 正电子对撞：$e^+e^- \to 2j$ 和 $3j$。
   • 轻子 - 强子对撞：中性流（NC）和带电流（CC）深度非弹性散射（DIS）过程（$ep \to \ell + 2j$，$ep \to \nu + 2j$ 等）。
   • 强子 - 强子对撞：包容性喷注产生（$pp \to j$）、双喷注产生（$pp \to 2j$）、矢量玻色子加喷注（$V+1j$）、光子加喷注（$\gamma+1j$）、双光子（$\gamma\gamma$）以及希格斯玻色子加喷注（$H+1j$）产生。
   • 衍生阶数：该框架隐式地为具有额外一个喷注的过程提供 NLO 修正，为具有额外两个喷注的过程提供 LO 修正（例如，通过具有 2 喷注要求的 $ZJ$过程计算$Z+2j$ 的 NLO 结果）。

2. 灵活的用户界面（Runcard）：代码通过纯 ASCII 格式的"runcard"文件进行控制，允许用户定义：

   • 过程与对撞机：选择对撞机类型（pp, ep, $e^+e^-$）、束流能量以及喷注算法（anti-$k_T$、Cambridge-Aachen、$k_T$、Durham、Jade）。
   • 可观测量与选择：定义唯象截断（fiducial cuts）、对象级选择器（喷注 $p_T$、快度）以及事件级选择器。
   • 直方图：灵活的分箱（线性、对数、自定义边界）以及通过直方图选择器实现的多微分分布。
   • 标度：同时评估多种重整化标度和因子化标度组合以进行不确定性估计。

3. 自动化工作流：发行版包含一个基于 Python 的工作流脚本（nnlojet-run），用于管理作业提交、资源分配和结果合并。该工作流处理：

   • 预热与生产阶段：自适应网格生成，随后进行生产运行。
   • 资源管理：动态调整计算资源以达到目标相对精度。
   • 结果合并：一种稳健的算法（涉及修剪、k 扫描和加权组合），用于合并来自独立作业的局部结果，同时减轻异常值的影响。

4. 验证：该实现包含严格的验证检查：

   • 点测试：验证在奇异极限附近，矩阵元与减法项的比率趋近于 1。
   • 技术截断检查：确保结果独立于技术红外截断。
   • 极点测试：显式抵消 RV 和 VV 贡献中的 $\epsilon$ 极点。

结果与验证
本文详细阐述了各种过程的实现，并引用了具体出版物以提供验证细节。作为一个具体示例，作者展示了LHC 上 $Z + 1j$ 产生过程（$\sqrt{s} = 7$ TeV）的计算：

• 设置：Anti-$k_T$ 喷注算法（$R=0.4$）、NNPDF3.1 NNLO PDFs 以及对轻子和喷注的特定唯象截断。
• 结果：计算得到的总 NNLO 唯象截面为 68.5(2) pb。
• 微分分布：文章展示了领头喷注的横向动量分布，显示了从 LO 到 NLO 再到 NNLO 的收敛性，其中不确定性带代表了 7 点标度变化。

意义与主张
作者声称，NNLOJET 提供了一个普遍适用、开源的框架，用于喷注可观测量的 NNLO QCD 计算，解决了以往代码通常具有专有性质或仅限于特定过程类别的空白。

• 可及性：通过在 GPLv3.0 下分发代码，并提供用户友好的 runcard 语法和工作流脚本，作者旨在使高精度 NNLO 预测对更广泛的唯象学家和实验物理学家群体开放。
• 灵活性：无需重新编译代码即可定义任意唯象截断和可观测量的能力，使得与实验数据的直接比较成为可能。
• 可扩展性：模块化设计允许在未来的版本中添加新过程和特性。

本文结论指出，NNLOJET 是用于精密物理研究、对标准模型进行压力测试以及完善 SM 参数和部分子分布函数确定的成熟工具。该代码可在 HEPForge 获取，作者欢迎用户反馈和错误报告。