Measurement of the Z $\to$ $\mu^+\mu^-$ angular coefficients in pp collisions at $\sqrt{s}$ = 13 TeV as functions of transverse momentum and rapidity

本文利用 CMS 实验的 140 fb$^{-1}$数据，对 13 TeV 质子 - 质子碰撞中 Drell-Yan $Z \to \mu^+\mu^-$ 过程的八个角极化系数（$A_0$ 至 $A_7$）进行了测量，给出了横向动量和快度的双微分结果，并将其与次次领头阶 QCD 预言进行了比较。

要点

以下是该论文的通俗化解释，辅以生动的类比。

宏观图景：宇宙弹球机

想象一下，位于欧洲核子研究中心（CERN）的大型强子对撞机（LHC）是世界上威力最强大的弹球机。科学家们以接近光速的速度将质子相互撞击。通常，这些碰撞会产生一团混乱的粒子。然而，有时碰撞会生成一种重而不稳定的粒子，称为Z 玻色子，它会立即分裂成两个缪子（电子的重型“表亲”）。

本文关注的是CMS，这是监视这台弹球机的巨型探测器之一。研究团队不仅统计了这种事件发生的次数，更旨在理解缪子是如何飞出的。它们是径直射出？还是旋转飞出？亦或是偏向某个特定方向？

目标：测量碰撞的“自旋”

科学家们测量了八个不同的数值（标记为 $A_0$ 到 $A_7$）。可以将这些数字视为关于 Z 玻色子在爆炸前“姿态”或极化的详细成绩单。

• 类比：想象烟花在天空中爆炸。如果它完美对称地爆炸，火花会呈完美的球形向外飞溅。如果它是倾斜的或在旋转，火花可能会更多地向左、向上射出，或者呈螺旋状飞出。
• 测量：这八个系数（$A_0$ 到 $A_7$）告诉我们爆炸的“形状”究竟是怎样的。它们揭示了 Z 玻色子是否在旋转、摇摆，或者是否被“拉伸”在某个特定方向上。

方法：双重核查

团队分析了 2016 年至 2018 年间记录的140 万亿次碰撞（140 fb⁻¹ 的数据）。他们并非只看整堆数据，而是像切面包一样将其切片，以观察“自旋”是否随质子相互撞击的剧烈程度而变化。

1. 速度（横向动量）：他们观察了那些横向移动缓慢的缪子与那些移动极快的缪子。
2. 角度（快度）：他们观察了那些径直向前飞行的缪子与那些以锐角飞行的缪子。

通过在这些特定切片中测量缪子的角度，他们能够以极高的精度计算出这八个系数。

游戏规则："Lam-Tung"法则

本文讨论了一个著名的物理法则，称为Lam-Tung 关系。

• 类比：想象一条规则说：“如果你把球垂直向上扔，它必须垂直落下。”在粒子物理世界中，在最简单的计算层面上，这两个系数（$A_0$ 和 $A_2$）应该完美地相互抵消（$A_0 - A_2 = 0$）。
• 现实：论文证实，在低速度下，这条规则表现良好；但随着碰撞能量增加（动量更高），规则开始失效。这并非失败，而是一个特征！它表明碰撞中那些“混乱”的部分（例如被踢出的额外粒子）开始变得重要起来。

结果：数据与理论

科学家们将他们的测量结果与现有的最佳计算机模拟（“理论预测”）进行了比较。

• 好消息：对于大多数系数而言，现实世界的数据与计算机模型非常吻合。这意味着我们目前对这些粒子相互作用的理解是坚实的。
• 有趣的张力：在中等速度范围内，特定系数（$A_0$）的数据略高于计算机的预测（偏离约 3 个标准差）。这就像天气预报预测有 50% 的降雨概率，但实际上 80% 的时间都在下雨。这并非灾难，但它暗示计算机模型可能遗漏了微小的细节。
• “幽灵”系数：三个系数（$A_5, A_6, A_7$）理论上应为零或非常接近零。数据显示它们确实极小，与零一致，尽管其中一个（$A_6$）显示出极其微弱、非零的迹象。这就像在安静的房间里听到一声低语；它确实存在，但你需要非常灵敏的耳朵才能听到。

为何这很重要

这篇论文本质上是对物理定律的一次高精度校准检查。

1. 理解“胶水”：这些测量有助于我们理解“部分子动力学”——即质子内部微小构建块（夸克和胶子）在碰撞时的行为。
2. 检验理论通过将 Z 玻色子的“自旋”与复杂的数学（量子色动力学）进行比较，科学家们正在对我们对宇宙的理解进行压力测试。如果数学与数据不匹配，就意味着我们需要发明新的物理学。
3. 基准：这篇论文为未来的实验提供了一个新的、超精确的“标尺”。任何新理论都必须能够解释这八个数字。

总结

简而言之，CMS 团队拍摄了粒子碰撞的巨大快照，测量了由此产生的粒子的确切角度，并计算出了描述事件“自旋”的八个数字。他们发现，虽然我们的现有理论大体正确，但在中等速度范围内存在一些微小而迷人的差异，这让物理学家们时刻保持警惕，确保对宇宙更深层次理解的探索仍在继续。

技术摘要

以下是 CMS 合作组论文《在√s = 13 TeV 的质子 - 质子碰撞中测量 Z → µ+µ− 角系数随横向动量和快度的变化》的详细技术总结。

1. 问题与动机

Drell–Yan (DY) 过程（$pp \to \gamma^*/Z \to \ell^+\ell^-$）是理解质子 - 质子碰撞底层部分子动力学以及电弱矢量玻色子产生机制的基本探针。在玻色子静止系中，末态轻子的角分布编码了关于虚光子或 Z 玻色子极化状态的信息。

该分布由在 Collins–Soper (CS) 框架中定义的八个角极化系数 $A_0$ 至 $A_7$ 进行参数化。

• 物理目标：
  • $A_0$ 和 $A_2$： 参数化纵向极化以及横向态之间的干涉。它们的差值（$A_0 - A_2$）检验Lam–Tung 关系（$A_0 - A_2 = 0$），该关系在 QCD 领头阶（LO）成立，但预期会因高阶 QCD 辐射和部分子内禀横向动量而被破坏。
  • $A_4$： 与前向 - 后向不对称性（$A_{FB}$）相关，用于测量有效弱混合角（$\sin^2 \theta_{eff}^W$）。
  • $A_1, A_3$： 编码纵向态与横向态之间的干涉，并依赖于轴矢量/矢量耦合。
  • $A_5, A_6, A_7$： 定义 T-奇不对称性。这些系数在 LO 和 NLO 预期为零，但由于单圈过程，在次次领头阶（NNLO）可能非零。

尽管此前在 $\sqrt{s}=8$ TeV 和向前区域（LHCb）已有测量，但仍需要在 $\sqrt{s}=13$ TeV 下对全套系数（$A_0$–$A_7$）进行高统计量的全面双微分测量，以约束理论模型并以更高精度检验 NNLO QCD 预测。

2. 方法论

数据与事例选择

• 数据集： CMS 探测器在 2016–2018 年期间于 $\sqrt{s} = 13$ TeV 下采集的质子 - 质子碰撞数据，对应积分亮度为 140 fb$^{-1}$。
• 信号： 不变质量 $81 < m_{\mu\mu} < 101$ GeV（Z 玻色子极点区域）的双轻子事例（$Z \to \mu^+\mu^-$）。
• 运动学 cuts：
  • 双缪子横向动量（$p_T^{\mu\mu}$）：0 至 400 GeV。
  • 双缪子快度（$|y_{\mu\mu}|$）：$< 2.4$。
  • 应用缪子隔离和鉴别标准以确保高纯度。
  • 触发要求：根据数据采集时期不同，采用 $p_T > 24$ 或 $27$ GeV 的单缪子触发。

模拟与修正

• 信号模拟： 使用 POWHEG v2.0 + MiNNLOPS 生成，并接口 PYTHIA 8（用于部分子簇射/强子化）和 PHOTOS++（用于 QED 末态辐射）。使用了 NNPDF3.1 部分子分布函数（PDF）集。
• 背景估计： 对双玻色子（WW, WZ, ZZ）、$t\bar{t}$ 和 $W$+喷注过程进行蒙特卡洛模拟。背景贡献很小（低 $p_T$ 时为 0.03%，高 $p_T$ 时升至约 2.9%）。
• 修正：
  • 尺度因子： 通过“标记 - 探测”（tag-and-probe）方法导出，用于修正触发、鉴别和隔离效率。
  • 堆积重加权： 使模拟的堆积分布与数据匹配。
  • 运动学重加权： 在 $(p_T^{\mu\mu}, y_{\mu\mu})$ 区间内对模拟事例进行重加权，以匹配数据的运动学分布，从而最小化模型依赖性。
  • 预触发修正： 解决 Level-1 触发中束团交叉的对准偏差。

系数提取

提取过程利用模板法来考虑探测器接受度、分辨率以及 $p_T$ 区间间的迁移。

1. 定义： 微分截面展开为三角多项式 $f_i(\theta^*, \phi^*)$ 之和，由系数 $P_i$ 加权（其中 $A_i = P_i/P_8$）。
2. 模板： 从重建的蒙特卡洛事例中构建九个二维模板（$T_i$）。每个模板对应特定的角项 $f_i$，在事件生成中加权以去除极化信息，同时保留探测器效应。
3. 拟合： 在每个运动学区间（8 个 $p_T$ 区间 × 2 个快度区间，共 16 个）的观测角分布（$\cos \theta^*, \phi^*$）上执行二维分箱最大似然拟合（使用 MINUIT）。拟合通过最小化数据与模板线性组合之间的差异来提取系数 $P_i$。

3. 主要贡献

• 13 TeV 下的首套完整测量： 这是首次在 $\sqrt{s}=13$ TeV 的中心快度区域（$|y_{\mu\mu}| < 2.4$）测量全套角系数（$A_0$–$A_7$）。
• 双微分精度： 结果在 8 个 $p_T^{\mu\mu}$ 区间和 2 个 $|y_{\mu\mu}|$ 区间中提供，呈现了角结构的多维视图。
• NNLO 比较： 结果与最先进的 QCD NNLO 理论预测（使用 POWHEG+MiNNLOPS 和 MADGRAPH5_aMC@NLO）进行了比较。
• T-奇不对称性约束： 对中心区域的 T-奇系数（$A_5, A_6, A_7$）提供了迄今为止最严格的约束。

4. 结果

总体一致性

• 测量得到的系数在不确定度范围内通常与 NNLO QCD 预测一致。
• 大数据集（140 fb$^{-1}$）显著降低了统计不确定度，优于之前的 Run 1 测量。

具体观测

• $A_0$ 和 $A_2$：
  • 在中间 $p_T$ 范围（10–35 GeV），$A_0$ 显示数据与 MADGRAPH5_aMC@NLO 预测之间存在约 3 个标准差 的差异，而 POWHEG+MiNNLOPS 吻合良好。
  • 在高 $p_T$（$>50$ GeV）下，MC 模拟略微高估了 $A_2$。
• Lam–Tung 关系（$A_0 - A_2$）：
  • 正如预期因高阶 QCD 效应所导致，该关系被破坏。
  • 在 10–35 GeV 范围内，POWHEG+MiNNLOPS 预测与测量的 $A_0 - A_2$ 值吻合良好。
  • 对于 $p_T > 35$ GeV，两种 MC 预测均低估了测量值。
• $A_1$：
  • 在 $p_T < 35$ GeV 且 $1 < |y_{\mu\mu}| < 2.4$ 区域，POWHEG+MiNNLOPS 显著高估了 $A_1$。
  • MADGRAPH5_aMC@NLO 也在两个快度区域的低 $p_T$ 处高估了 $A_1$。
• T-奇系数（$A_5, A_6, A_7$）：
  • 数值量级为 $10^{-3}$，且在大多数区间内与零一致。
  • 最显著的偏差出现在 $A_6$，区间为 $55 < p_T < 80$ GeV 且 $1.2 < |y_{\mu\mu}| < 2.4$，显示出相对于零的 2.8$\sigma$ 偏差。这与 NNLO 预期一致，但其幅度小于之前 ATLAS 在 8 TeV 下的观测。
• 系统不确定度：
  • 在低 $p_T$（$<55$ GeV）下主导 $A_0$ 和 $A_1$ 的不确定度，在 $p_T < 200$ GeV 下主导 $A_2, A_3, A_4$ 的不确定度。
  • 对于其他系数和更高 $p_T$ 区域，统计不确定度占主导地位。

5. 意义

• 理论验证： 结果为 NNLO QCD 计算提供了严格的基准，证实了高阶修正对于描述 Drell–Yan 产生的角结构至关重要。
• 部分子动力学： 这些测量为部分子分布函数（PDFs）和部分子内禀横向动量的建模提供了新的约束。
• 未来物理： 此处建立的精度和多维分箱为未来的现象学研究以及测试更高阶计算（如 N$^3$LO）提供了关键参考。
• 探测器性能： 展示了 CMS 探测器即使在存在复杂探测器效应和堆积的情况下，也能以高精度测量微妙的角关联的能力。

总之，这篇论文代表了精确电弱物理的一个里程碑，证实了 NNLO QCD 对 Z 玻色子产生机制描述的有效性，同时突出了可能需要进一步完善的理论模型的具体运动学区域。