Study of the Magnetic Dipole Transition of $J/ψ\toγη_c$ via $η_c\to p\bar{p}$

基于 BESIII 探测器收集的约 100 亿 $J/\psi$ 事例，该研究首次对 $J/\psi\to\gamma p\bar{p}$ 过程进行了振幅分析，精确测定了相关分支比，并由此计算出 $J/\psi\to\gamma\eta_c$ 和 $\eta_c\to\gamma\gamma$ 的分支比，其结果与最新格点量子色动力学计算一致。

要点

这篇论文就像是一次高能物理界的“精密侦探破案”，由中国的 BESIII 实验团队（北京谱仪 III 合作组）完成。他们利用巨大的粒子加速器，解开了一个困扰物理学家多年的谜题：关于一种叫"J/ψ"的粒子如何变成另一种叫"ηc"的粒子，并在这个过程中发射出一个光子。

为了让你轻松理解，我们可以把整个实验过程想象成一场**“超级粒子魔术秀”**。

1. 舞台与演员：巨大的粒子工厂

想象一下，北京有一个巨大的环形跑道（BEPCII 加速器），里面有两股“粒子流”——正电子和负电子，它们像两列高速对撞的火车，在特定的地点猛烈相撞。

• J/ψ粒子：当它们相撞时，会产生一种非常重、非常不稳定的“明星粒子”，叫 J/ψ。你可以把它想象成一个**“超重且充满活力的杂技演员”**。
• 实验数据：这次魔术秀非常盛大，他们收集了超过 100 亿次 这样的碰撞事件。这就像是在一场盛大的马戏团里，记录了 100 亿次杂技演员的表演瞬间。

2. 魔术过程：变身与发光

在这个实验中，科学家们想观察 J/ψ 演员的一个特殊动作：

• 动作描述：J/ψ 演员突然“瘦身”并变身，从一种状态（J/ψ）变成了另一种状态（ηc），同时吐出一个光子（γ，也就是一束光）。
• ηc 的结局：变身后的 ηc 演员非常不稳定，它几乎瞬间就会“爆炸”成一对质子和反质子（p 和 p̄）。
• 侦探任务：科学家们的任务就是在这 100 亿次表演中，精准地找出那些“吐光变身”并“炸成质子对”的特定瞬间。

3. 难点：在噪音中找信号

这就好比在一个嘈杂的万人体育馆里，你想找出那一个特定的魔术师变出的特定鸽子。

• 背景噪音：除了我们要找的 J/ψ 变身，还有很多其他杂乱的粒子碰撞（背景噪音），它们也会产生质子和光子，看起来很像，但其实是“冒牌货”。
• 传统方法 vs. 新方法：
  • 以前，科学家像用**“单筒望远镜”**看，只能大概估算有多少个这样的魔术，容易看走眼，或者被背景噪音干扰。
  • 这次，BESIII 团队用上了**“全景 3D 摄像机”**（这就是论文中提到的“振幅分析”）。他们不仅看结果，还详细分析了每一个粒子飞出的角度、能量和相互关系。这就像不仅数鸽子，还分析了魔术师挥手的角度、帽子的倾斜度，从而能极其精准地把“真魔术”和“假把式”区分开。

4. 破案成果：精确测量

通过这种高精度的“全景分析”，他们得到了两个惊人的结果：

1. 找到了“真魔术”的数量：他们极其精确地算出了 J/ψ 变成 ηc 并放出光子的概率（分支比）。这个精度比以前的记录提高了10 倍！就像以前只能猜“大概有 100 次”，现在能精确到“就是 100.2 次”。
2. 解开了“理论 vs. 现实”的矛盾：
   • 谜题：过去几十年，物理学家用超级计算机（格点量子色动力学，LQCD）模拟这个魔术，算出来的结果总是比实验测出来的大一倍。就像理论预测魔术师应该变出 200 只鸽子，但实验只看到了 100 只。这一直是个大谜团。
   • 新发现：这次 BESIII 用新方法测出来的结果，完美地落在了理论预测的范围内！
   • 意义：这就像侦探终于找到了关键证据，证明了之前的理论模型（超级计算机模拟）其实是对的，而以前的实验测量可能因为方法不够精细（没排除噪音）而看错了。

5. 总结：为什么这很重要？

这篇论文不仅仅是在数数粒子，它是在验证我们理解宇宙基本力量的规则。

• 强相互作用：J/ψ 和 ηc 是由夸克组成的，它们之间的变身过程涉及“强相互作用”（把原子核粘在一起的力）。
• 量子色动力学 (QCD)：这是描述强相互作用的理论。这次实验结果与最新的理论计算一致，意味着我们对宇宙中最基础、最复杂的“胶水”力量的理解又向前迈进了一大步。

一句话总结：
BESIII 团队利用 100 亿次粒子碰撞数据，用一种全新的“全景透视”技术，精准捕捉了粒子变身的瞬间，不仅把测量精度提高了 10 倍，还成功消除了理论与实验之间长达几十年的巨大分歧，证明了我们的宇宙理论模型是靠谱的。

技术摘要

这是一份关于 BESIII 合作组最新研究成果《通过 $\eta_c \to p\bar{p}$ 研究 $J/\psi \to \gamma\eta_c$ 磁偶极跃迁》的详细技术总结。该论文旨在解决重夸克偶素系统中磁偶极（M1）跃迁宽度理论与实验长期存在的显著差异问题。

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

• 核心物理问题：$J/\psi \to \gamma\eta_c$ 是重夸克偶素系统中两个最低能态（$1^3S_1$ 和 $1^1S_0$）之间的磁偶极（M1）跃迁。这是检验量子色动力学（QCD）在微扰和非微扰区域理论的理想实验室。
• 理论与实验的矛盾：
  • 在非相对论极限下，理论预测的跃迁宽度比实验测量值大 2 到 3 倍。
  • 尽管有多种理论尝试（如色散求和规则、QCD 求和规则、相对论夸克模型、有效场论、光锥求和规则等）来解释这一差异，但格点 QCD (LQCD) 的最新计算结果系统性地比粒子数据组 (PDG) 的平均实验值大约大两倍。
  • 此前 BESIII 在 $J/\psi \to \gamma\eta_c$ 中测量 $\eta_c \to \gamma\gamma$ 分支比时，发现若使用 PDG 的 $J/\psi \to \gamma\eta_c$ 分支比，计算出的 $\eta_c \to \gamma\gamma$ 分支比与 LQCD 预测及 PDG 全局拟合值存在超过 $3\sigma$ 的偏差。
• 现有测量的局限性：之前的实验测量（包括 CLEO-c, KEDR, Crystal Ball 及 BESIII 早期工作）多采用一维拟合，往往忽略了非共振（NR）振幅与 $\eta_c$ 共振态之间的干涉效应，或者未包含所有可能的非共振贡献（如 $J^P C \neq 0^{-+}$ 的态），这可能导致测量结果的偏差。

2. 研究方法 (Methodology)

• 数据来源：利用 BESIII 探测器在 BEPCII 对撞机上采集的 $(10.087 \pm 0.044) \times 10^9$ 个 $J/\psi$ 事例。
• 分析通道：首次对 $J/\psi \to \gamma p\bar{p}$ 过程进行振幅分析 (Amplitude Analysis)，选取 $p\bar{p}$ 不变质量 $M_{p\bar{p}}$ 在 $\eta_c$ 质量区域 $[2.70, 3.05] \text{ GeV}/c^2$ 范围内的事例。
• 事件选择与背景抑制：
  • 要求两个带电径迹（质子/反质子）电荷代数和为零，并满足动量、飞行时间 (TOF) 和比电离能 ($dE/dx$) 的粒子鉴别 (PID) 要求。
  • 光子候选者需满足能量阈值及电磁量能器 (EMC) 的时间窗要求，并排除质子/反质子附近的簇射以抑制次级光子背景。
  • 实施 4C 运动学拟合（能量动量守恒），选择 $\chi^2$ 最小的组合。
  • 主要背景来源：$J/\psi \to p\bar{p}\gamma_{FSR}$ (FSR 光子) 和 $J/\psi \to p\bar{p}\pi^0$。通过控制样本和蒙特卡洛 (MC) 模拟精确估计背景贡献（分别约占 22.8% 和 4.2%）。
• 振幅模型构建：
  • 采用协变张量振幅形式，包含 $J/\psi$ 极化矢量、光子极化矢量及 $p\bar{p}$ 自旋波函数。
  • 共振态描述：$\eta_c$ 线形采用相对论性 Breit-Wigner 函数，质量 $M_{\eta_c}$ 和宽度 $\Gamma_{\eta_c}$ 为自由参数。
  • 阻尼因子：引入 KEDR 形式的阻尼因子（优于 CLEO-c 形式）以抑制 $\eta_c$ 长尾发散。
  • 非共振 (NR) 贡献：考虑了所有可能的角动量和自旋耦合的 NR 波（共 13 种潜在贡献）。通过似然比检验，保留了统计显著性大于 $3\sigma$ 的 5 个分量（包括 $0^{-+}, 1^{++}, 2^{++}, 2^{--}$ 等波）。
  • 干涉处理：在拟合中显式地处理了 $\eta_c$ 与 NR 振幅之间的干涉，特别是 $\eta_c$ 与 $0^{-+}$ NR 贡献之间的相位角。
• 拟合策略：使用最大似然法拟合，构建包含数据项和背景项的对数似然函数。通过 Toy MC 样本评估拟合偏差。

3. 主要贡献与结果 (Key Contributions & Results)

A. 物理量测量结果

1. $\eta_c$ 共振参数：

   • 质量：$M_{\eta_c} = (2984.55 \pm 0.09_{\text{stat}} \pm 0.77_{\text{syst}}) \text{ MeV}/c^2$
   • 宽度：$\Gamma_{\eta_c} = (29.74 \pm 0.17_{\text{stat}} \pm 0.33_{\text{syst}}) \text{ MeV}$
   • 结果与 PDG 世界平均值高度一致，且精度显著提高。

2. 分支比乘积测量：

   • 首次通过振幅分析精确测定：
     $$\mathcal{B}(J/\psi \to \gamma\eta_c) \times \mathcal{B}(\eta_c \to p\bar{p}) = (2.11 \pm 0.02_{\text{stat}} \pm 0.07_{\text{syst}}) \times 10^{-5}$$
   • 相比之前的测量，精度提高了约一个数量级。

3. 独立分支比提取：
   结合 PDG 平均的 $\mathcal{B}(\eta_c \to p\bar{p}) \times \mathcal{B}(\eta_c \to \gamma\gamma)$ 和 BESIII 最新测得的 $\mathcal{B}(J/\psi \to \gamma\eta_c) \times \mathcal{B}(\eta_c \to \gamma\gamma)$，利用时间反演不变性，解耦得到：

   • $J/\psi \to \gamma\eta_c$ 分支比：
     $$\mathcal{B}(J/\psi \to \gamma\eta_c) = (2.29 \pm 0.01_{\text{stat}} \pm 0.04_{\text{syst}} \pm 0.18_{\text{opbf}})\%$$
     (其中 opbf 代表来自其他分支比输入的不确定性)
   • $\eta_c \to \gamma\gamma$ 分支比：
     $$\mathcal{B}(\eta_c \to \gamma\gamma) = (2.28 \pm 0.01_{\text{stat}} \pm 0.04_{\text{syst}} \pm 0.18_{\text{opbf}}) \times 10^{-4}$$
   • $\eta_c \to p\bar{p}$ 分支比：
     $$\mathcal{B}(\eta_c \to p\bar{p}) = (0.92 \pm 0.01 \pm 0.02 \pm 0.07) \times 10^{-3}$$

B. 关键发现

• 解决长期矛盾：新测得的 $\mathcal{B}(J/\psi \to \gamma\eta_c)$ 和 $\mathcal{B}(\eta_c \to \gamma\gamma)$ 与 PDG 全局拟合值存在约 $3\sigma$ 的偏差，但与最新的格点 QCD (LQCD) 计算结果（如 HPQCD 和 CLQCD 合作组）高度吻合。
• 质量分裂：基于 $J/\psi$ 和 $\eta_c$ 的质量，测得 $c\bar{c}$ 系统的自旋 - 轨道分裂为 $(112.35 \pm 0.77) \text{ MeV}/c^2$，与理论计算一致。

4. 科学意义 (Significance)

1. 验证 LQCD 理论：该结果强有力地支持了格点 QCD 对重夸克偶素 M1 跃迁宽度的预测，表明之前的实验与理论之间的巨大差异可能源于实验分析中对非共振干涉效应的处理不当。
2. 方法论突破：首次在该过程中应用全振幅分析，成功分离了共振态与非共振背景，并精确量化了干涉效应。这为未来研究类似的强子衰变过程提供了重要的分析范式。
3. 解决“谜题”：这一发现有望解决重夸克偶素物理中长达数十年的 M1 跃迁宽度理论与实验不符的谜题，确认了非相对论势模型和 LQCD 在描述强相互作用非微扰区域的有效性。
4. 精度提升：将相关分支比的测量精度提升至 10% 以内，为未来更精确的 QCD 检验和寻找新物理提供了坚实的基础数据。

总结：BESIII 合作组通过大规模数据样本和先进的振幅分析技术，重新测定了 $J/\psi \to \gamma\eta_c$ 及其后续衰变的分支比。新结果不仅修正了 PDG 的推荐值，更重要的是消除了实验数据与前沿格点 QCD 理论之间的显著矛盾，标志着我们对强相互作用非微扰性质的理解迈出了重要一步。