Search for the charmonium weak decay $ψ(2S)\to D_s^-π^+ + c.c.$ and $ψ(2S)\to D_s^-ρ^+ + c.c.$

BESIII 合作组利用 $(2712.4\pm14.3)\times 10^6$ 个 $\psi(2S)$ 事例，首次对 $\psi(2S)\to D_s^-\pi^+ + c.c.$ 和 $\psi(2S)\to D_s^-\rho^+ + c.c.$ 弱衰变过程进行了搜索，未发现信号超出，并给出了相应的分支比上限。

要点

这篇论文就像是一份来自粒子物理侦探社的“寻人启事”和“安全报告”。

简单来说，BESIII 合作组（一支由全球科学家组成的超级侦探团队）利用巨大的“粒子对撞机”（BEPCII）和超级相机（BESIII 探测器），试图捕捉一种极其罕见、甚至可能从未被人类亲眼见过的“粒子变身”现象。

以下是用大白话和生动比喻对这篇论文的解读：

1. 侦探的任务：寻找“不可能”的变身

想象一下，$\psi(2S)$ 是一个性格非常稳定的“粒子大叔”。

• 他的常规操作：通常，他要么保持原样，要么通过“强力”（强相互作用）把自己拆成两个更小的粒子（比如 $D\bar{D}$ 对）。这就像是他把家里的家具（能量）分给两个孩子，大家都能吃饱。
• 物理定律的限制：但是，这个“大叔”的体重（质量）有点尴尬，刚好比“拆成两个大孩子”所需的能量少了一点点。所以，按照常规玩法，他根本没法把两个大孩子都生出来。
• 这次的任务：科学家们想看看，这位大叔会不会偷偷用一种极慢、极微弱的“弱力”（弱相互作用），把自己的一部分能量“变”成一个带 charm 夸克的粒子（$D_s^-$）和一个普通粒子（$\pi^+$ 或 $\rho^+$）。
  • 这就像是一个胖子，明明不够钱买两张机票（能量不够），却试图通过某种极其隐秘的魔法（弱力），把自己变成一张机票加一个行李箱。
  • 为什么这很重要？ 在标准模型（物理界的“官方剧本”）里，这种变身发生的概率极低，低到几乎不可能被现在的设备抓到（就像中彩票头奖的概率）。但如果真的抓到了，而且频率比剧本写的要高，那就说明**“剧本”错了**，或者宇宙里藏着**“新物理”**（比如超对称粒子、额外维度等）。

2. 侦探的装备：超级显微镜

为了抓这个“幽灵”，BESIII 团队动用了他们的看家本领：

• 海量数据：他们收集了约 27 亿 个 $\psi(2S)$ 粒子事件。这就像是在大海里捞针，而且他们把整个大海都捞了一遍。
• 精密筛选：
  • 他们不直接看“大叔”怎么变，而是看变出来的“孩子”（$D_s^-$）去了哪里。
  • $D_s^-$ 很调皮，它会迅速衰变成其他粒子（比如 $\phi$ 介子、电子 $e^-$ 和中微子 $\nu$）。
  • 中微子是个“隐形人”，探测器看不见它。所以侦探们玩了一个“反向推理”的游戏：既然看不见中微子，那我就算算“大叔”原本的能量和动量，减去所有看得见的“孩子”（电子、介子等），剩下的“失踪人口”就是中微子。如果算出来的“失踪人口”符合中微子的特征，那就有戏了！

3. 调查过程：盲测与排雷

为了不让科学家自己的“先入为主”影响判断，他们采用了**“盲测”**（Blind Analysis）：

• 第一步：先不看真实数据，只在电脑模拟（蒙特卡洛模拟）里把筛选规则定好。
• 第二步：只挑出 10% 的真实数据来“试跑”一下，看看规则有没有漏洞。
• 第三步：最后，把剩下的 90% 数据（真正的“宝藏”）打开，看看有没有抓到“鱼”。

在这个过程中，他们还要排除各种“捣乱分子”（背景噪音）：

• 有些粒子衰变看起来很像目标，但其实只是普通的强相互作用产物。
• 有些是宇宙射线或者机器本身的杂讯。
• 他们像过滤网一样，一层层把这些噪音过滤掉，只留下最纯净的候选者。

4. 调查结果：虽然没有抓到“大鱼”，但排除了“嫌疑”

经过严密的筛选和计算，结果出来了：

• 现状：在 27 亿个样本里，没有发现明显的“信号过剩”。也就是说，没有看到比背景噪音多出来的“异常变身”事件。
• 结论：虽然没抓到“新物理”的实锤，但这并不是失败。
  • 科学家们给出了一个**“上限”（Upper Limit）：如果这种变身真的存在，它的概率最多**只能是多少多少（$\psi(2S) \to D_s^- \pi^+$ 的概率小于 $1.4 \times 10^{-6}$，$\psi(2S) \to D_s^- \rho^+$ 的概率小于 $7.0 \times 10^{-6}$）。
  • 这个上限虽然比标准模型预测的（$10^{-10}$ 级别）要高，说明现在的设备还不够灵敏，还没法完全排除标准模型的预测，但也没有发现那些“新物理模型”预言的、能增强几百倍的信号。

5. 总结：下一步怎么办？

这就好比侦探在案发现场转了一圈，虽然没找到凶手，但排除了很多嫌疑人，并且确认了案发地点的安保级别（灵敏度）还不够高，抓不住那种“隐形”的凶手。

• 好消息：标准模型目前依然坚挺，没有被推翻。
• 坏消息：现有的数据量还不够大，灵敏度还不够高，无法探测到那些极其微小的“新物理”效应。
• 未来计划：科学家们表示，“我们需要更多的数据！” 只要继续积累更多的粒子碰撞数据，或者升级更灵敏的探测器，未来就有希望真正揭开这个“粒子大叔”是否拥有隐藏魔法的秘密。

一句话总结：
BESIII 团队用 27 亿个粒子做了一次超级大海捞针，虽然这次没捞到传说中的“新物理”金条，但他们成功画出了一张更精确的“藏宝图”，告诉全世界：如果金条真的存在，它一定藏在这个范围之外，我们需要挖得更深、更仔细才行。

技术摘要

以下是基于 BESIII 合作组论文《Search for the charmonium weak decay $\psi(2S) \to D_s^- \pi^+ + \text{c.c.}$ and $\psi(2S) \to D_s^- \rho^+ + \text{c.c.}$》的详细技术总结：

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

• 物理背景：$\psi(2S)$ 共振态的质量略低于 $D\bar{D}$ 对产生的阈值，因此其强衰变到 $D\bar{D}$ 在运动学上是禁戒的。然而，在标准模型（SM）中，$\psi(2S)$ 通过弱相互作用衰变到单个粲介子（如 $D$ 或 $D_s$）是允许的。
• 科学挑战：
  • 极低的分支比：标准模型预测 $\psi(2S)$ 弱衰变的分支比（BF）极低（约为 $10^{-10}$ 量级或更低），远低于当前实验的探测能力。
  • 新物理探针：尽管 SM 预测极低，但多种新物理模型（如 Top-color 模型、最小超对称标准模型 MSSM、双希格斯二重态模型等）预测这些弱衰变的分支比可能比 SM 预期高出 2-3 个数量级。
  • 研究空白：此前实验尚未观测到粲偶素的弱衰变。BESIII 此前已对 $J/\psi$ 的弱衰变进行了搜索，但针对 $\psi(2S)$ 的特定模式（$D_s^- \pi^+$ 和 $D_s^- \rho^+$）尚未进行系统性搜索。
• 核心目标：利用 BESIII 探测器收集的大样本数据，首次搜索 $\psi(2S) \to D_s^- \pi^+ + \text{c.c.}$ 和 $\psi(2S) \to D_s^- \rho^+ + \text{c.c.}$ 这两个弱衰变道，以检验标准模型并寻找新物理迹象。

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

• 数据来源：
  • 利用 BESIII 探测器在质心能量 $\sqrt{s} = 3.686$ GeV 处收集的 $(2712.4 \pm 14.3) \times 10^6$ 个 $\psi(2S)$ 事例。
  • 使用基于 GEANT4 的蒙特卡洛（MC）模拟样本来确定选择标准、计算探测效率并估计背景。
• 衰变道重建策略：
  • 信号特征：两个衰变道最终都产生四个带电粒子：$K^+, K^-, \pi^+, e^-$（对于 $\rho^+$ 模式，$\pi^0$ 衰变为 $\gamma\gamma$）。
  • $D_s^-$ 重建：由于 $D_s^-$ 包含中微子，无法完全重建。采用半轻子衰变道 $D_s^- \to \phi e^- \bar{\nu}_e$（其中 $\phi \to K^+ K^-$），利用反冲质量法（Recoil Method）来识别 $D_s^-$ 候选者。
    • 计算 $D_s^-$ 的能量和动量：$E_{D_s^-} = E_{\psi(2S)} - E_{\pi^+/\rho^+}$，$\vec{p}_{D_s^-} = \vec{p}_{\psi(2S)} - \vec{p}_{\pi^+/\rho^+}$。
    • 定义反冲质量 $M_{D_s^-}$ 并设定信号窗口。
  • 粒子鉴别 (PID)：结合飞行时间（TOF）、电离能损（dE/dx）和电磁量能器（EMC）能量沉积（$E/p$）来区分 $K, \pi, e$。
  • 光子重建：针对 $\rho^+$ 模式中的 $\pi^0 \to \gamma\gamma$，对光子候选者施加能量和角度约束。
• 背景抑制：
  • 缺失动量与能量：利用中微子特性，要求缺失动量 $|\vec{p}_{miss}| > 0.02$ GeV/c 和缺失能量变量 $U_{miss}$ 接近零，以抑制不含中微子的强衰变背景。
  • 额外光子抑制：针对 $\psi(2S) \to D_s^- \rho^+$ 模式，额外要求除 $\pi^0$ 衰变光子外的剩余光子总能量 $E_{\gamma}^{rest} < 0.12$ GeV，以抑制多光子背景。
• 分析方法：
  • 采用逐步盲分析（Stepwise Blind Analysis）策略，先基于 MC 确定选择标准，再用 10% 的随机数据验证，最后分析全样本，以避免分析者偏差。
  • 背景估计分为两部分：来自 $\psi(2S)$ 衰变的背景（利用包容性 MC 估计）和连续谱背景（利用邻近能量点的数据估计）。
  • 由于未观测到显著信号，采用**似然函数（Profile Likelihood）**方法计算 90% 置信水平（C.L.）下的分支比上限。

3. 关键贡献 (Key Contributions)

• 首次搜索：这是人类历史上首次对 $\psi(2S) \to D_s^- \pi^+$ 和 $\psi(2S) \to D_s^- \rho^+$ 这两个特定弱衰变道进行实验搜索。
• 高精度上限设定：基于目前世界上最大的 $\psi(2S)$ 数据样本，给出了这两个衰变道最严格的分支比上限。
• 系统误差评估：详细评估了包括 MC 模型、跟踪效率、PID 效率、中间分支比、光子探测效率等在内的多种系统误差来源，确保了结果的可靠性。
• 盲分析实践：在稀有衰变搜索中严格实施了盲分析流程，保证了统计结果的客观性。

4. 研究结果 (Results)

• 观测事例数：
  • 在 $\psi(2S) \to D_s^- \pi^+$ 信号区域内观测到 9 个事例。
  • 在 $\psi(2S) \to D_s^- \rho^+$ 信号区域内观测到 19 个事例。
• 背景估计：
  • $\psi(2S) \to D_s^- \pi^+$ 的估计背景为 $5.6 \pm 2.3$ 个事例。
  • $\psi(2S) \to D_s^- \rho^+$ 的估计背景为 $12.4 \pm 3.2$ 个事例。
  • 观测值与预期背景一致，未发现显著的信号超出。
• 分支比上限 (90% C.L.)：
  • $\mathcal{B}(\psi(2S) \to D_s^- \pi^+ + \text{c.c.}) < 1.4 \times 10^{-6}$
  • $\mathcal{B}(\psi(2S) \to D_s^- \rho^+ + \text{c.c.}) < 7.0 \times 10^{-6}$
• 系统误差：总系统误差分别为 10.1% ($\pi^+$ 模式) 和 14.2% ($\rho^+$ 模式)。

5. 意义与结论 (Significance)

• 与标准模型对比：实验给出的上限（$10^{-6}$ 量级）虽然比标准模型预测值（$10^{-10}$ 量级）高出约 4 个数量级，但与标准模型预测一致，即目前的实验数据尚未发现新物理存在的证据。
• 新物理探测能力：结果表明，当前的 $\psi(2S)$ 数据样本灵敏度还不足以探测到那些预测分支比增强 2-3 个数量级的新物理模型效应。
• 未来展望：
  • 该研究为未来更高统计量的实验（如超级 $\tau$-charm 工厂）设定了基准。
  • 结论明确指出，为了探测潜在的弱衰变信号或新物理效应，需要积累更多的数据。
  • 该工作填补了粲偶素弱衰变实验研究的空白，为理解强子弱衰变机制提供了重要的实验约束。

总结：BESIII 合作组利用海量 $\psi(2S)$ 数据，通过精密的盲分析策略，首次对 $\psi(2S)$ 的弱衰变进行了严格搜索。虽然未观测到信号，但给出了目前最严格的分支比上限，证实了标准模型的低分支比预测，并强调了未来更高亮度实验的必要性。