Observation of $ψ(3686) \to n\bar{n}$ and improved measurement of $ψ(3686) \to p \bar{p}$

该研究利用 BESIII 探测器采集的 1.07 亿个 $ψ(3686)$ 事例，首次观测到 $ψ(3686) \to n\bar{n}$ 衰变并以前所未有的精度测量了 $ψ(3686) \to p\bar{p}$ 的分支比及角分布参数。

要点

这篇论文讲述的是物理学家们利用中国北京的高能物理实验装置（BESIII），进行了一次非常精彩的“粒子捉迷藏”游戏。他们主要做了两件事：第一次成功“抓”到了 $\psi(3686)$ 粒子衰变成一对中子和反中子（$n\bar{n}$）的踪迹，并且更精确地测量了它衰变成一对质子和反质子（$p\bar{p}$）的情况。

为了让你更容易理解，我们可以把这篇论文的内容想象成一场**“宇宙级的高能台球赛”**。

1. 舞台与主角：巨大的粒子台球桌

• 主角 $\psi(3686)$：你可以把它想象成一个非常重、能量很高的“台球”。它是由夸克组成的（具体说是粲夸克和反粲夸克）。
• 实验场地：这个台球是在北京正负电子对撞机（BEPCII）上产生的。就像两个台球以极高的速度相撞，产生了一个新的、更重的“母球” $\psi(3686)$。
• 探测器 BESIII：这就像是一个巨大的、360 度无死角的超级高速摄像机，包围着碰撞点。它的作用是记录母球碎裂后飞出来的所有碎片。

2. 任务一：第一次抓到“隐形”的中子对

背景故事：
以前，物理学家已经知道这个“母球”可以分裂成一对质子和反质子（$p\bar{p}$）。质子带电，就像带磁铁的台球，很容易在摄像机里留下轨迹，所以以前很容易抓到。

但是，中子（$n$）和反中子（$\bar{n}$）是不带电的。它们就像隐形的幽灵台球。它们在摄像机里不会留下轨迹，只有当它们撞到探测器最外层的“墙壁”（电磁量能器）时，才会突然爆发出一团能量闪光，然后消失。

挑战：
因为中子太“隐形”了，而且背景噪音（其他杂乱的粒子碰撞）太多，以前没人能确认 $\psi(3686)$ 到底有没有变成过中子对。这就像在嘈杂的夜店里，试图听清一个轻声细语的人说话。

解决方案（AI 侦探）：
为了抓住这些“幽灵”，物理学家们用了一种叫**“多变量分析”（MVA）的高级算法，这就像是一个超级 AI 侦探**。

• 它学习了成千上万个“幽灵”出现时的特征（比如能量爆发的形状、位置等）。
• 它学会了区分哪些是真正的“幽灵”（信号），哪些是普通的“杂音”（背景）。
• 最终，这个 AI 侦探从 1.07 亿次碰撞中，成功识别出了 6000 多次 $\psi(3686) \to n\bar{n}$ 的事件。这是人类第一次正式确认这个衰变过程的存在！

3. 任务二：更精确地测量“带电”的质子对

既然已经能抓“幽灵”了，物理学家回头重新检查了那些容易抓的“带电台球”（质子和反质子）。

• 以前大家测过很多次，但不够精确。
• 这次，利用海量的数据（1.07 亿次碰撞），他们把测量结果做得非常精准，就像用显微镜重新观察了一个熟悉的物体，发现了一些以前忽略的细节。

4. 核心发现：台球飞出的“姿势”

这篇论文不仅数了数有多少个台球，还研究了它们飞出来的角度。

• 比喻：想象母球爆炸后，碎片是像烟花一样向四面八方均匀散开，还是像被挤压过一样，主要沿着某个方向飞出去？
• 公式 $1 + \alpha \cos^2\theta$：物理学家用一个参数 $\alpha$ 来描述这种“飞散姿势”。
  • 如果 $\alpha = 1$，说明碎片像被压扁的橄榄球，主要沿着碰撞轴线飞。
  • 如果 $\alpha = 0$，说明碎片是均匀散开的。
• 结果：
  • 对于质子对（$p\bar{p}$），$\alpha$ 接近 1.03。这意味着它们飞出的姿势非常“标准”，符合某些理论预测（就像两个台球沿着直线弹开）。
  • 对于中子对（$n\bar{n}$），$\alpha$ 是 0.68。这意味着它们的飞散姿势和质子对不一样！

5. 为什么这很重要？（解开谜题）

这就好比我们在研究“为什么两个双胞胎兄弟（质子和中子）在同样的爆炸中，摆出的姿势却不同？”

• 12% 规则：物理界有一个著名的“经验法则”，认为 $\psi(3686)$ 衰变的概率应该是 $J/\psi$（另一个较轻的粒子）衰变概率的 12% 左右。这次测量发现，无论是质子对还是中子对，都符合这个规则。这就像验证了“能量守恒”在微观世界依然有效。
• 未解之谜：虽然概率符合规则，但姿势（$\alpha$ 值）不同却让人困惑。
  • 在 $J/\psi$ 衰变中，质子和中子的姿势是一样的。
  • 但在 $\psi(3686)$ 衰变中，它们不一样了。
  • 这暗示着在 $\psi(3686)$ 这个能量级别，强相互作用（把夸克粘在一起的力）和电磁相互作用（光子交换）之间，可能有一种复杂的“舞蹈”或“相位差”，以前我们没完全看懂。

总结

这篇论文就像是一次高精度的宇宙考古：

1. 新发现：第一次挖出了“中子对”这个隐藏的宝藏。
2. 旧物新测：把“质子对”的测量精度提升到了新高度。
3. 核心启示：发现虽然概率符合预期，但粒子飞出的“姿势”揭示了微观世界强相互作用中尚未被完全理解的复杂机制。

这就好比我们一直以为两个双胞胎兄弟走路姿势完全一样，直到有一天发现，在特定的高速奔跑中，哥哥（质子）是正步走，弟弟（中子）却是侧着走。这个发现将帮助物理学家们重新编写关于“物质如何结合”的教科书。

技术摘要

这是一份关于 BESIII 合作组论文《Observation of $\psi(3686) \to n\bar{n}$ and improved measurement of $\psi(3686) \to p\bar{p}$》（arXiv:1803.02039v2）的详细技术总结。

1. 研究背景与问题 (Problem)

• 物理背景：量子色动力学（QCD）在高能区已被充分验证，但在低能区非微扰效应占主导，理论计算极为复杂。$\psi(3686)$（即 $\psi'$）介子处于微扰与非微扰区域的过渡区，研究其强子和电磁衰变有助于理解该能区的强相互作用结构。
• 核心问题：
  1. 首次观测缺失：尽管 $\psi(3686) \to p\bar{p}$ 的衰变已被测量近四十年，但 $\psi(3686) \to n\bar{n}$ 的衰变从未被观测到。
  2. 对称性检验：测量这两个过程可以检验味 SU(3) 对称性，并确定强相互作用与电磁相互作用振幅之间的相对相位角。
  3. 角分布参数 $\alpha$：矢量共振态衰变到重子 - 反重子对的角分布遵循 $dN/d\cos\theta \propto 1 + \alpha \cos^2\theta$。之前的测量（如 E835 和 BESII）对 $\psi(3686) \to p\bar{p}$ 的 $\alpha$ 值测量精度较低且存在不确定性，且缺乏 $\psi(3686) \to n\bar{n}$ 的 $\alpha$ 值数据。
  4. 12% 规则与 $\rho\pi$ 谜题：微扰 QCD 预测 $\psi(3686)$ 与 $J/\psi$ 衰变到同一末态的分支比之比应约为 12.7%（"12% 规则”）。然而，该规则在某些道（如 $\rho\pi$）中失效。精确测量 $N\bar{N}$ 末态有助于理解这一谜题。

2. 实验数据与方法论 (Methodology)

• 数据来源：

  • 使用 BESIII 探测器在 BEPCII 对撞机上采集的数据。
  • 信号样本：$1.07 \times 10^8$ 个 $\psi(3686)$ 事例。
  • 本底样本：$44 \text{ pb}^{-1}$ 的离共振（off-resonance）数据（质心能量 3.65 GeV），用于确定非共振本底。
  • 蒙特卡洛模拟 (MC)：使用 $1.07 \times 10^8$ 个 inclusive $\psi(3686)$ 模拟事例研究本底，信号 MC 根据角分布 $1+\alpha\cos^2\theta$ 生成。

• $\psi(3686) \to n\bar{n}$ 分析策略：

  • 特征：中子和反中子不带电，在探测器中不留下径迹，主要通过电磁量能器（EMC）中的簇射（shower）被探测。反中子因发生湮灭，在 EMC 中的能量沉积通常高于中子。
  • 事件选择：要求 MDC 中无重建径迹，仅依赖 EMC 信息。定义“第一簇射群”（SG1，能量最高）和“第二簇射群”（SG2）。
  • 多变量分析 (MVA)：由于中子/反中子探测复杂，采用基于 Boosted Decision Tree (BDT) 的多变量分析。
    • 输入变量：共 14 个变量，包括 SG1 和 SG2 的能量、击中数、二阶矩、侧向矩、最大开口角、剩余簇射的数量和能量等。
    • 训练：使用信号 MC 和加权混合的本底数据（离共振数据、inclusive MC、特定衰变道 MC）进行训练。
    • 优化：BDT 切割条件优化后，保留了约 76% 的信号效率，同时拒绝了约 95% 的本底。
  • 效率修正：利用控制样本 $\psi(3686) \to p\bar{n}\pi^-$ 研究数据与 MC 在中子/反中子探测效率上的差异，并作为 $\cos\theta$ 的函数进行修正。

• $\psi(3686) \to p\bar{p}$ 分析策略：

  • 特征：质子和反质子带电，在 MDC 中留下径迹。
  • 事件选择：要求两条带电径迹，动量在 $1.546 < p < 1.628 \text{ GeV}/c$ 范围内，且开口角 $\theta_{open} > 3.1$ rad。
  • 粒子鉴别 (PID)：利用飞行时间探测器（TOF）的似然值区分质子/反质子与 K 介子。
  • 本底估计：主要本底来自连续过程 $e^+e^- \to p\bar{p}$，利用离共振数据外推；非 $p\bar{p}$ 的 $\psi(3686)$ 衰变本底经 MC 检查可忽略。
  • 效率修正：利用 $\psi(3686) \to p\bar{p}$ 样本本身，通过单径迹重建和反冲侧径迹重建的方法，确定追踪和 PID 效率，并作为 $\cos\theta$ 的函数进行修正。

• 拟合与参数提取：

  • 对 $\theta_{open}$ 分布进行拟合以提取信号事例数。
  • 对 $\cos\theta$ 分布进行拟合，函数形式为 $N_{sig}(1+\alpha\cos^2\theta)\epsilon(\theta) + N_{bg}f_{bg}$，提取 $\alpha$ 参数。

3. 主要结果 (Key Results)

基于 $1.07 \times 10^8$ 个 $\psi(3686)$ 事例，研究得出以下结果：

A. 分支比 (Branching Fractions)

• $\psi(3686) \to n\bar{n}$：首次观测到该衰变。
  • $B(\psi(3686) \to n\bar{n}) = (3.06 \pm 0.06 \text{ (stat)} \pm 0.14 \text{ (syst)}) \times 10^{-4}$
• $\psi(3686) \to p\bar{p}$：测量精度显著提高。
  • $B(\psi(3686) \to p\bar{p}) = (3.05 \pm 0.02 \text{ (stat)} \pm 0.12 \text{ (syst)}) \times 10^{-4}$

B. 角分布参数 $\alpha$

假设角分布服从 $1 + \alpha \cos^2\theta$：

• $\psi(3686) \to n\bar{n}$：
  • $\alpha_{n\bar{n}} = 0.68 \pm 0.12 \text{ (stat)} \pm 0.11 \text{ (syst)}$
• $\psi(3686) \to p\bar{p}$：
  • $\alpha_{p\bar{p}} = 1.03 \pm 0.06 \text{ (stat)} \pm 0.03 \text{ (syst)}$
  • 该结果接近 1.0，且精度优于以往所有测量（如 E835, BESII, CLEO 等）。

C. 其他发现

• 奇次项检验：检查 $2\gamma$ 交换过程导致的 $\cos\theta$ 奇次项（$\beta \cos\theta$），结果显示 $\beta$ 与零一致，未观测到显著贡献。
• 12% 规则检验：
  • $B(\psi(3686) \to p\bar{p}) / B(J/\psi \to p\bar{p}) = (14.4 \pm 0.6)\%$
  • $B(\psi(3686) \to n\bar{n}) / B(J/\psi \to n\bar{n}) = (14.8 \pm 1.2)\%$
  • 两个比值均与微扰 QCD 预测的 12.7% 在误差范围内一致。

4. 系统误差分析 (Systematic Uncertainties)

论文详细评估了多种系统误差来源，包括：

• $\psi(3686)$ 事例数：0.7%。
• 探测效率修正：中子/反中子效率修正的不确定性较大（对 $\alpha_{n\bar{n}}$ 影响约 12.8%），质子/反质子效率修正影响较小。
• 本底估计：连续谱本底和 inclusive MC 本底的不确定性。
• 分辨率与重建：动量分辨率、$\theta_{open}$ 分辨率、簇射变量（如剩余簇射能量）的模拟偏差。
• 触发效率：针对中性事例和带电事例的触发效率评估。
• 物理模型：$\alpha$ 值变化对效率计算的影响。

5. 意义与结论 (Significance)

1. 填补空白：首次观测到 $\psi(3686) \to n\bar{n}$ 衰变，填补了重子对产生谱系中的重要空白。
2. 精度提升：将 $\psi(3686) \to p\bar{p}$ 的分支比和 $\alpha$ 参数的测量精度提升至新水平，特别是 $\alpha_{p\bar{p}} \approx 1.03$ 的结果支持了 QCD 螺旋度守恒规则（$\alpha=1$）在 $\psi(3686)$ 能区的适用性，尽管之前的测量倾向于 $\alpha < 1$。
3. 物理机制启示：
   • 在 $J/\psi$ 衰变中，$p\bar{p}$ 和 $n\bar{n}$ 的分支比和 $\alpha$ 值非常接近，暗示强相互作用主导且相对相位接近 90°。
   • 在 $\psi(3686)$ 衰变中，虽然分支比依然接近，但 $\alpha_{n\bar{n}}$ (0.68) 与 $\alpha_{p\bar{p}}$ (1.03) 存在显著差异。这表明 $\psi(3686) \to N\bar{N}$ 的衰变机制可能比 $J/\psi$ 更复杂，可能涉及长程过程或连续谱成分的贡献，使得简单的相位角提取变得困难。
4. 理论挑战：结果与"12% 规则”一致，但 $\alpha$ 值的差异为理解强子化过程中的非微扰效应和可能的“强子过剩”（hadronic excess）提供了新的实验约束，有助于解决 $\rho\pi$ 谜题相关的理论争议。

综上所述，该工作利用 BESIII 的大统计量数据，通过先进的多变量分析技术，实现了对 $\psi(3686)$ 衰变到核子对的首次精确测量，为理解强相互作用在过渡能区的性质提供了关键实验数据。