A subthreshold pole in electron-positron annihilation into the $D_s^+D_s^-$ final state

该研究利用 BESIII 合作组的精确数据，在$e^+ e^-\rightarrow D_s^+ D_s^-$过程中发现了质量为$(3896^{+13}_{-22})$ MeV 的亚阈极点，并有力支持将其识别为在$D\bar D$过程中表现为共振态的$G(3900)$粒子。

要点

这篇论文讲述了一个关于微观粒子世界的“捉迷藏”故事。为了让你更容易理解，我们可以把高能物理实验想象成在一个巨大的粒子游乐场里进行的一场精密的“听音辨位”游戏。

1. 背景：游乐场里的“听音”游戏

想象一下，科学家们在**北京正负电子对撞机（BEPCII）**这个巨大的游乐场里，让两束粒子（正电子和负电子）像碰碰车一样高速相撞。

• 碰撞结果：当它们撞在一起时，会瞬间产生能量，并“变”出新的粒子。在这个实验中，科学家特别关注它们变成了什么——主要是$D_s^+ D_s^-$ 这种成对的粒子（你可以把它们想象成游乐场里的一对对“双胞胎”）。
• 任务：科学家通过测量不同能量下产生这些“双胞胎”的数量（也就是“截面”），试图画出一种**“反应地图”。这张地图应该显示出，在什么能量下，产生粒子的概率会突然变大（就像听到一声巨响），这通常意味着那里有一个“共振态”**（Resonance），也就是一个不稳定的新粒子在“跳舞”。

2. 遇到的问题：地图上的“幽灵”

科学家（Peter Lichard 和 Josef Jurán）拿着最新的、非常精确的数据（来自 BESIII 合作组），试图用现有的理论模型（就像用标准的乐谱）来解释这张地图。

• 常规操作：他们先假设地图上的高峰是由几个已知的“跳舞粒子”（共振态）造成的。
• 发现异常：无论他们怎么调整已知的粒子参数，地图的起步阶段（也就是能量刚刚超过产生粒子的最低门槛时）总是对不上号。数据在这里的表现，就像有一个看不见的幽灵在推波助澜，让反应概率比理论预测的要高。

3. 核心发现：亚阈值的“隐形 pole"

为了解释这个“幽灵”，科学家提出了一个大胆的想法：这里可能藏着一个**“亚阈值极点”（Subthreshold Pole, SP）**。

• 什么是“亚阈值”？
  想象你要推一辆很重的车（产生一对粒子），你需要达到一定的力气（能量阈值）才能推得动。

  • 普通粒子：力气够了，车动了，你看到了它。
  • 亚阈值极点：这个“幽灵”需要的力气比推车的最低门槛还要小。也就是说，在能量还没达到能产生这对粒子的程度时，这个“幽灵”就已经存在了。
  • 为什么看不见？ 因为它太“稳定”了（或者说它处于一种特殊的束缚状态），在能量不够产生粒子时，它不会像普通共振态那样“爆炸”出来被直接看到。它就像藏在门槛下面的一个隐形弹簧，虽然你看不到它，但它会弹起，影响上面推车的过程。

• 证据确凿：
  科学家把数学模型改了一下，在这个“门槛下面”加了一个这样的“隐形弹簧”（极点）。奇迹发生了！

  • 原本对不上的数据，现在完美吻合了。
  • 统计数据显示，这个发现不是运气，其可信度高达 7.4 个标准差（在科学界，这相当于“铁证如山”，比发现希格斯玻色子时的证据还要强）。
  • 这个“幽灵”的质量大约是 3896 MeV，它位于产生 $D_s^+ D_s^-$ 粒子所需能量的下方（大约低 42 MeV）。

4. 这个“幽灵”是谁？

既然找到了这个“隐形弹簧”，它到底是谁呢？

• 嫌疑犯 A：它可能是 $D_s^+ D_s^-$ 分子的一个束缚态（就像两个粒子手拉手，还没完全分开，但能量又低到不能算作自由粒子）。
• 嫌疑犯 B（最可能的）：它极有可能是著名的 G(3900) 粒子。
  • 在其他的实验中（比如产生普通 D 介子对时），G(3900) 表现得像一个普通的、会“跳舞”的共振态（你能听到它的声音）。
  • 但在产生 $D_s^+ D_s^-$ 这个特定过程中，它却“隐身”了，变成了门槛下的“亚阈值极点”。
  • 比喻：就像同一个人，在热闹的舞会上（普通 D 介子过程）会大声唱歌（共振），但在安静的图书馆里（$D_s^+ D_s^-$ 过程）却只能压低声音，甚至躲在桌子底下（亚阈值极点），虽然你看不到他，但他依然影响着周围的气氛。

5. 总结与意义

这篇论文告诉我们：

1. 世界比看起来更复杂：在粒子物理中，有些东西虽然能量不够“现身”，但它们依然真实存在，并且对物理过程有巨大的影响。
2. 新视角：通过这种“听音辨位”的高级数学分析，我们发现了 G(3900) 粒子在不同环境下的不同“性格”。
3. 未来：这暗示 G(3900) 可能是一个真正的强子态（由夸克组成的粒子），并且它与 $D_s^+ D_s^-$ 系统有着非常强的联系。

一句话总结：
科学家利用最新的精密数据，发现了一个藏在能量门槛之下的“隐形粒子”（亚阈值极点），它完美解释了实验数据的异常，并很可能就是那个著名的 G(3900) 粒子在特定情况下的“隐身”形态。这就像在门槛下发现了一个看不见的弹簧，虽然看不见，但它确实把门推开了。

技术摘要

这是一份关于论文《A subthreshold pole in electron-positron annihilation into $D^+_s D^-_s$ final state》（电子 - 正电子湮灭产生 $D^+_s D^-_s$ 末态中的亚阈值极点）的详细技术总结。

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

• 核心问题：在电子 - 正电子湮灭（$e^+e^-$）产生强子的过程中，除了通常观测到的共振态（表现为复平面第二黎曼叶上的共轭极点）外，是否存在位于第一黎曼叶（物理叶）实轴上、但在反应阈值以下的“亚阈值极点”（Subthreshold Poles, SP）？
• 物理意义：亚阈值极点对应的是对特定末态稳定的束缚态。虽然它们无法直接在实验中作为共振峰被观测到（因为能量低于阈值），但根据解析性（causality and unitarity），它们会显著影响物理区域的散射振幅和激发函数形状。
• 具体目标：利用 BESIII 合作组发布的最新高精度数据，研究 $e^+e^- \to D^+_s D^-_s$ 过程，寻找并确认是否存在亚阈值极点，并探讨其物理本质（如是否与 $G(3900)$ 态有关）。

2. 方法论 (Methodology)

• 数据来源：使用了 BESIII 合作组在 BEPCII 对撞机上测量的 $e^+e^- \to D^+_s D^-_s$ 的 Born 截面数据。数据覆盖了从阈值（3.938 GeV）到 4.951 GeV 的 138 个能量点，包含统计误差和系统误差。
• 理论模型：
  • 基于矢量介子主导模型 (VMD) 的修正公式来描述 $e^+e^- \to D\bar{D}$ 过程。
  • 截面公式 $\sigma(s)$ 包含多个中间矢量介子（共振态）的贡献，形式为：
    $$\sigma(s)= \frac{\pi \alpha^2}{3s} \left(1-\frac{4m^2_{D_s}}{s}\right)^{3/2} \left| \sum_{k=1}^n \frac{\sqrt{Q_k} e^{i\delta_k}}{s - M^2_k + iM_k\Gamma_k} \right|^2$$
  • 关键创新：在拟合中引入亚阈值极点 (SP)。将 SP 视为宽度 $\Gamma=0$ 且位于实轴 $s < 4m^2_{D_s}$ 处的极点。
• 拟合策略：
  • 采用标准 $\chi^2$ 最小化准则拟合实验数据。
  • 逐步增加共振态数量（从单共振到多共振），并引入亚阈值极点进行对比。
  • 比较不同模型（纯共振模型 vs. 共振 + 亚阈值极点模型）的拟合优度（$\chi^2/DOF$ 和 $p$ 值）及统计显著性（通过似然比变化 $\delta(-2 \ln L)$ 计算）。

3. 关键贡献与结果 (Key Contributions & Results)

A. 亚阈值极点的发现

• 初步发现：在尝试用四个共振态拟合数据时，第四个“共振”的参数显示其质量 $M \approx 3891$ MeV，宽度 $\Gamma \approx 0$ MeV。这表明该结构并非普通共振，而是一个位于 $D^+_s D^-_s$ 阈值以下的亚阈值极点。
• 最佳拟合模型 (SP+4R)：
  • 引入一个亚阈值极点（SP）和四个共振态（4R）的模型（SP+4R fit）获得了最佳拟合效果。
  • 拟合优度：$\chi^2/DOF = 124.5/120$，$p$-value 提升至 37%（相比之下，纯 3 共振模型 $p$-value 仅为 $5\times 10^{-8}$）。
  • 统计显著性：亚阈值极点的存在具有 7.4$\sigma$ 的统计显著性（在 SP+3R 模型中）；第四个共振态的显著性为 5.5$\sigma$。

B. 极点参数

• 亚阈值极点质量：
  • 在 SP+3R 拟合中：$M = 3890 \pm 13$ MeV。
  • 在 SP+4R 拟合中：$M = 3896 \pm 12$ MeV。
  • 经过更详细的误差分析（$\chi^2_{min}+1$ 准则），最终确定质量为 $M = (3896^{+13}_{-22})$ MeV。
• 结合能：该极点对应 $D^+_s D^-_s$ 系统的束缚态，结合能 $B = 2m_{D_s} - M \approx (42^{+22}_{-13})$ MeV。
• 量子数：由于出现在 $e^+e^-$ 湮灭中，其量子数为 $I^G(J^{PC}) = 0^-(1^{--})$。

C. 与其他态的对比与排除

• 与 PDG 数据的对比：该研究使用的 BESIII 最新数据与 PDG 2024 基于旧数据（>18 年）得出的 $\psi(4040)$ 参数存在显著差异。新拟合得到的 $\psi(4040)$ 质量 ($4026.4 \pm 1.6$ MeV) 与之前 BES 和 BESIII 的联合分析结果高度一致。
• 排除分子态解释：该极点的量子数 ($1^{--}$) 排除了它是 $D^+_s D^-_s$ 分子基态（应为 $0^{++}$）的可能性，也不符合激发态的预期。
• 与 $G(3900)$ 的关联：
  • 该亚阈值极点的质量与 BESIII 在 $e^+e^- \to D\bar{D}$ 过程中发现的 $G(3900)$ 共振态质量 ($3872.5 \pm 14.2 \pm 3.0$ MeV 或全局耦合道分析中的 $3898.4 \pm 0.9$ MeV) 高度吻合。

4. 结论与意义 (Significance)

1. 证实亚阈值极点的存在：该研究提供了强有力的证据，证明在 $e^+e^- \to D^+_s D^-_s$ 过程的振幅中存在一个显著的亚阈值极点。这一发现强调了在拟合强子产生截面时，考虑阈值以下极点的重要性，仅靠共振态无法完美描述数据。
2. 重新诠释 $G(3900)$：研究提出，$G(3900)$ 可能是一个真实的强子态，具有强烈的 $D^+_s D^-_s$ 耦合。
   • 在 $e^+e^- \to D\bar{D}$ 过程中，它表现为一个共振态。
   • 在 $e^+e^- \to D^+_s D^-_s$ 过程中，由于能量低于阈值，它表现为一个亚阈值极点。
3. 方法论启示：展示了利用高精度实验数据结合解析性约束（VMD 模型扩展），可以有效探测那些在实验上不可直接观测（低于阈值）但对物理过程有重要影响的量子态。
4. 物理图景：暗示 $G(3900)$ 可能是一个具有轨道角动量 $L=1$ 的 $D^+_s D^-_s$ 激发态，或者是某种具有强耦合特性的奇特强子态，而非简单的分子基态。

总结：这篇论文通过重新分析 BESIII 的高精度数据，利用改进的 VMD 模型，首次明确指出了 $D^+_s D^-_s$ 产生过程中的亚阈值极点，并将其与 $G(3900)$ 态联系起来，为理解重夸克偶素及奇特强子态的性质提供了新的视角。