Logarithmically-accurate showers with massive quarks

本文提出了一种包含大质量夸克的 PanScales 末态部分子簇射方案，该方案在保持对质量无关可观测量原有精度的同时实现了次领头对数精度，并通过固定阶检验、全阶重求和比较以及利用 LEP 数据的唯象研究验证了其有效性。

要点

想象宇宙的基本粒子如同一个混乱的高速舞池。在这场舞蹈中，像底夸克和粲夸克这样的重粒子，就像穿着沉重笨重靴子的舞者，而较轻的粒子则如同穿着芭蕾舞鞋的舞者。

几十年来，物理学家一直使用称为“部分子簇射”的计算机模拟来预测这些舞者如何移动和相互作用。然而，大多数此类模拟都将所有人视为穿着芭蕾舞鞋，忽略了重舞者移动方式不同的事实。它们遗漏了舞池的一条关键规则：“死锥”。

死锥：个人空间气泡

当一名重舞者（大质量夸克）旋转时，他们无法像轻舞者那样自由地挥动手臂（辐射能量）。由于他们很沉重，他们在周围形成了一个个人空间气泡，其他人无法靠得太近。这就是“死锥”。如果模拟忽略了这一点，就会预测在重舞者附近辐射出过多的能量，从而导致对舞池形态的错误预测。

新方案：PanScales 簇射

本文作者构建了一种名为PanScales的新颖、更智能的模拟。这相当于升级了舞池的规则手册，将重舞者的靴子纳入其中。

他们并没有仅仅为死锥添加一个简单的“禁止入内”标志，而是重写了舞蹈的物理法则，以确保：

1. 重舞者保持能量：因为他们无法在个人气泡中轻易地辐射能量，所以他们保留了更多原始速度。
2. 轻舞者仍正常起舞：新规则不会扰乱对轻舞者的预测；它们只改变那些重靴子起作用的地方。
3. 数学完美无缺：他们通过将新规则与最精确的数学公式进行比对（就像将食谱与主厨的精确测量进行核对），证明了新规则的有效性。

他们如何测试

为了确保新模拟不仅仅是猜测，他们进行了三种类型的测试：

• “两步”测试（固定阶）：他们模拟了恰好包含两个额外动作的舞蹈，并将其与精确的数学答案进行了比较。即使涉及重舞者，他们的模拟也与数学完美吻合。
• “人群流动”测试（全阶）：他们观察了人群的整体形状（“隆德树”）。他们检查模拟是否正确预测了人群中形成了多少子群。同样，新模拟准确无误，捕捉到了质量如何微妙地改变流动的方式。
• “现实世界”核查（现象学）：他们将模拟与来自 LEP 对撞机（著名的过去粒子加速器）的真实数据进行了比较。他们观察了底夸克如何碎裂成其他粒子。
  • 结果：旧的模拟（忽略质量）预测底夸克会减速并过早碎裂。而新的 PanScales 模拟尊重了“重靴子”，与真实世界数据的吻合度要高得多。

核心结论

本文介绍了一种模拟粒子碰撞的新方法，终于以应有的尊重对待重粒子。它通过正确模拟“死锥”效应，修复了物理模拟中长期存在的盲点。

作者已将此新代码公开，允许其他科学家利用它更好地理解构成我们宇宙的重粒子，确保当我们查看来自大型强子对撞机（LHC）等巨型对撞机的数据时，看到的不是因忽略舞者重量而导致的扭曲画面。

技术摘要

技术摘要：含大质量夸克的对数精确部分子簇射

1. 问题陈述

重夸克（底夸克和粲夸克）在大型强子对撞机（LHC）等对撞机实验中被大量产生，对于理解希格斯玻色子和顶夸克物理至关重要。然而，与无质量情形相比，对包含大质量部分子的喷注中的量子色动力学（QCD）辐射进行建模面临特定挑战。最显著的特征是“死锥”（dead-cone）效应，即在重夸克方向周围大小为 $m_Q/E$ 的锥体内，共线发射受到抑制。

从微扰 QCD 的角度来看，质量尺度的存在在展开式中引入了新的对数增强项（例如 $\ln(m_Q/Q)$）。虽然针对特定可观测量（如碎裂函数）的对数重求和已在次领头对数（NLL）甚至次次领头对数（NNLL）精度下实现，但一个能够系统性地捕捉各类可观测量中的质量效应，同时保持无质量簇射原有精度的通用部分子簇射框架一直缺失。现有的簇射通常仅在领头对数（LL）精度下纳入质量效应，或者在质量相关时无法在全局事件形状上保持 NLL 精度。

2. 方法论

作者构建了一类新的末态部分子簇射，命名为PanScales，该簇射考虑了夸克质量并实现了 NLL 精度。该框架包含两个变体：

• PanLocal：基于偶极子（dipole）的簇射，具有完全局域的动量守恒。
• PanGlobal：类似天线（antenna）的簇射，具有全局横向动量守恒。

2.1 运动学映射与类光矢量

为了在保持 Sudakov 分解结构的同时处理大质量部分子，作者引入了由分支前大质量动量（$\bar{p}_i, \bar{p}_j$）导出的类光参考矢量（$\bar{n}_i, \bar{n}_j$）。这些矢量的定义使得在偶极子静止系中，它们的三维动量与大质量腿的三维动量相匹配，仅在能量分量上有所不同。这一选择确保了：

• 将类光矢量与大质量矢量相关联的系数在无质量极限下为零。
• 大质量粒子的方向得以保留，从而有利于准共线动量分数 $z$ 的有意义定义。
• 对数精度测试的数值稳定性得到改善。

2.2 发射概率

微分发射概率的构建旨在在两个关键极限下重现精确的 QCD 矩阵元：

1. 软极限：簇射必须重现大质量 eikonal 因子，其中包括正比于 $m^2/(p \cdot k)^2$ 的死锥抑制项。
2. 准共线极限：簇射必须重现大质量 DGLAP 分裂函数（$P_{Q \to Qg}$ 和 $P_{g \to Q\bar{Q}}$）。

作者通过修改发射概率因子 $D_i$ 和 $D_j$（及其 PanGlobal 对应项 $\bar{D}_i, \bar{D}_j$）来包含依赖于质量的项，从而抑制死锥区域内的辐射，以此实现这些极限。

2.3 有效耦合与味阈值

有效强耦合 $\alpha_s^{\text{eff}}$ 在变味数方案中定义。作者实施了一种特定的穿越重夸克阈值的方案：

• 由于 CMW 修正（$K_{CMW}$），有效耦合在略偏离物理质量 $m_Q$ 的标度 $\mu^{(n_f)}$ 处保持连续。
• 该偏移量在 $\mathcal{O}(\alpha_s^2)$ 阶推导为 $\ln \mu^{(n_f)} = \ln m_Q + 5/6$，从而确保跑动耦合具有正确的 NLL 行为。

2.4 无质量精度的保持

一个关键要求是，大质量簇射不得降低原始无质量公式的精度。

• PanGlobal 对全局可观测量保持 NNLL 精度。
• PanLocal 保持次次双对数（NNDL）精度。
  这是通过以下方式实现的：
• 使用无质量矩阵元实施次领头阶乘匹配，但在死锥区域内用大质量矩阵元替换簇射矩阵元。
• 使用可变味数（$n_f$）评估 NNLL Sudakov 修正。
• 通过确保双软矩阵元修正在大质量可忽略时退化为无质量结果，并对 $g \to Q\bar{Q}$ 通道在质量阈值以下进行否决（veto），来处理双软矩阵元修正。

2.5 自旋与色关联

• 自旋：Collins-Knowles 算法被改编以包含依赖于质量的螺旋度通道（例如，允许 $Q \to Qg$ 中的螺旋度翻转，以及 $g \to Q\bar{Q}$ 中的相等螺旋度）。
• 色：次领头色效应使用嵌套有序双软（NODS）方法处理，尽管作者指出死锥区域内并未完全计入次领头修正。

3. 验证与结果

作者进行了一系列全面的测试，以验证簇射的对数精度。

3.1 固定阶测试（$\mathcal{O}(\alpha_s^2)$）

• 相空间轮廓：作者验证了第二次发射不会以破坏 NLL 精度的方式改变先前发射的运动学。他们证明，虽然 Dire-v1 簇射在死锥区域未能通过此测试，但 PanScales 簇射（包括 PanLocal 和 PanGlobal）均满足要求，且在深共线区域的偏差被幂次压低（$1/k_t^4$）。
• 矩阵元比较：对于 $e^+e^- \to Q\bar{Q}g$ 和 $e^+e^- \to Q\bar{Q}g_1g_2$ 等过程，簇射权重与通过 MadGraph 生成的精确 QCD 矩阵元之间的微分比较显示，在相空间主体区域及死锥边界附近具有极好的一致性。仅在由幂次修正（$k_t^2/Q^2$）主导的区域或缺乏双软修正（这些修正为 $\mathcal{O}(N_c^{-1})$）的区域观察到偏差。
• 自旋关联：对顺序分裂的 $a_2/a_0$ 自旋关联比值的测试证实，无质量和有质量情形下均与解析预测完美一致。

3.2 全阶对数测试

• Lund 树形状（全局可观测量）：针对 Lund 平面中横向动量之和及最大值等可观测量，将簇射与 NLL 重求和结果进行了对比。结果证实，PanScales 簇射正确预测了死锥抑制和跑动耦合阈值，对于 $\beta_{ps} = 0.5$（PanLocal）和 $\beta_{ps} < 1$（PanGlobal）实现了 NLL 精度。
• 非全局能流：针对流入快度切片（rapidity slice）的能流（一种对非全局对数 NGLs 敏感的可观测量）测试了精度。簇射完美重现了解析重求和结果（包括耦合中的质量阈值），表明大质量 eikonal 因子的完整结构被正确捕捉。
• Lund 子喷注多重性：针对 $k_t$ 截断以上的子喷注数量的测试显示，对于轻夸克和重夸克引发的喷注，均与次双对数（NDL）解析计算结果一致。

3.3 唯象学研究

利用 LEP 数据（Z 峰事件），作者将簇射与实验测量进行了比较：

• $b$ 夸克碎裂函数：大质量簇射正确预测了数据中观察到的更硬的碎裂模式（更高的平均 $x_B$），而无质量簇射无法重现该形状，将平均值推向更低的值。
• 总喷注展宽：与无质量变体相比，大质量簇射在标度 $B_T \approx m_b/Q$ 附近与数据表现出更好的一致性，此处质量效应抑制了辐射。

4. 意义与主张

本文声称提出了首个明确考虑与非零夸克质量相关的 NLL 对数项的末态部分子簇射。

主要贡献包括：

1. 系统框架：一种统一的方法，用于捕捉不同类别可观测量（全局、非全局和多重性）中的质量效应，同时保持底层无质量簇射的高对数精度（NNLL/NNDL）。
2. 验证：通过高达 $\mathcal{O}(\alpha_s^2)$ 的固定阶测试以及与解析重求和的全阶比较进行了严格验证，确认了对死锥抑制、味阈值和自旋关联的正确处理。
3. 唯象学影响：证明了质量修正对于描述 $b$ 夸克碎裂函数至关重要，并且新簇射比无质量替代方案能更好地描述 LEP 数据。
4. 公开可用性：该簇射已实现在 PanScales 代码（v0.4.0）中，并与 Pythia 8 接口以进行强子化。

作者指出，虽然当前的实现达到了 NLL 精度，但它并未在简化的“死锥否决”模型中捕捉非全局对数（讨论见附录 C），这突显了其完整运动学和概率公式的必要性。未来的工作被确定为将这些算法扩展到初态辐射，并实施具有完整质量效应的高阶成分。