Amplitude Analysis and Branching Fraction Measurement of $D^+ \to \pi^+\pi^0\pi^0$

利用 BESIII 探测器采集的 20.3 fb$^{-1}$ $e^+e^-$ 对撞数据，本文首次对 $D^+ \to \pi^+\pi^0\pi^0$ 衰变进行了振幅分析，确认 $D^+ \to \rho(770)^+\pi^0$ 成分占主导地位，并报告了总分支比、中间拟合分数及 CP 不对称性的精确测量结果。

要点

想象一下，将亚原子世界比作一个高风险的舞池，粒子在其中碰撞、旋转，有时还会碎裂成更小的碎片。本文是一份来自 BESIII 实验的详细报告，该实验位于中国，相当于一个巨大的“相机”（探测器），它观测了数百万次这样的微小舞步，旨在理解一种特定的碎裂过程：一种名为$D^+$介子的粒子分裂成三个π介子（一种类似于质子“较轻表亲”的粒子）。

以下是他们发现的故事，无需复杂的数学公式即可理解。

1. 背景：一本庞大的相册

科学家们不仅仅拍了一张快照；他们汇编了一本庞大的相册。他们收集了来自203 亿次电子 - 正电子碰撞的数据（想象一下将两个微小的磁铁以接近光速的速度相互撞击）。如此庞大的数据量使他们能够观察到在较小样本中不可见的罕见事件。

他们的目标是研究衰变过程 $D^+ \to \pi^+ \pi^0 \pi^0$。

• $D^+$介子：开始舞蹈的舞者。
• $\pi^+$和两个$\pi^0$：舞者碎裂成的三个部分。

2. 谜团：碎裂是如何发生的？

当一个粒子碎裂成三块时，很少会一次性完成。通常，这是一个两步过程。想象一下，一位家长（$D^+$）将一个玩具拆分成三部分。

• 情景 A：家长先将玩具拆成一大块和一小块，然后大块再次碎裂。
• 情景 B：家长先将玩具拆成两块中等大小的部分，然后其中一块再次碎裂。

在物理学中，这些“块”被称为中间共振态。科学家们想知道：$D^+$介子走了哪条路径？

3. 主要发现：舞台上的“明星”

利用一种称为振幅分析的技术（这就像使用超级计算机，根据舞者的最终位置逆向推导舞蹈步骤），他们发现了一条特定的路径是绝对的赢家。

• 赢家：$D^+$介子几乎总是先转变成一个$\rho(770)^+$粒子和一个$\pi^0$。随后，$\rho(770)^+$迅速衰变成剩余的$\pi^+$和$\pi^0$。
• 类比：想象一位魔术师从帽子里变出一只兔子，但这只兔子实际上是一顶里面藏着更小兔子的帽子。这顶“大帽子”（$\rho$）是魔术发生的最常见方式。
• 结果：这条特定路径约占所有碎裂事件的63.5%。科学家们测量了这种发生的频率（即“分支比”），发现大约每1000个$D^+$介子中就有3个发生这种情况。

4. 配角

虽然$\rho(770)^+$是主角，但还有其他较少见的碎裂方式：

• $\rho$粒子的较重版本（$\rho(1450)$）。
• 一种名为$f_2(1270)$的不同粒子。
• 一个"S 波”态（一种模糊的、非共振的粒子云）。
• “干涉”效应：有时，这些不同的路径同时发生并相互干扰，就像两列声波相互抵消或增强噪音一样。科学家们测量了这些“干涉份额”，以了解不同路径是如何混合的。

5. “镜像”测试：寻找差异（CP 破坏）

物理学中最大的问题之一是：宇宙对待物质和反物质是否完全相同？

• $D^+$是物质。它的孪生兄弟$D^-$是反物质。
• 如果物理定律完全对称，$D^+$和$D^-$应该以完全相同的方式、相同的速率碎裂。
• 如果它们以不同的方式碎裂，则称为CP 破坏（这暗示宇宙略微偏爱物质而非反物质）。

结果：科学家们比较了$D^+$和$D^-$的“舞步”。他们发现没有显著差异。在误差范围内，两者的速率完全相同。

• 类比：这就像观察一位左撇子舞者和一位右撇子舞者表演完全相同的舞蹈。他们的手部动作略有不同，但整体的速度和风格是一样的。在这里没有发现任何“新物理”（如隐藏的力量）。

6. 这为何重要？

• 检验规则：理论物理学家建立了模型（如“极点模型”或“因子化”）来预测这些碎裂发生的频率。BESIII 的结果就像是这些模型的期末考试。
• 得分：主导路径（$\rho(770)^+$）与某些预测相符，但与另一些略有出入。这有助于科学家完善关于“强相互作用力”（将粒子粘合在一起的胶水）的理论，而该力的计算 notoriously 困难。
• 精度：通过测量这些事件的确切频率（约每1000次总衰变中有4.84次），他们为未来的实验提供了坚实的基准。

总结

BESIII 合作组利用庞大的粒子碰撞数据集，对$D^+$介子如何碎裂成三个π介子进行了详细的“法医分析”。他们发现，这种碎裂主要由涉及$\rho(770)^+$粒子的特定中间步骤主导。他们还证实，物质和反物质在此过程中表现完全一致，未发现任何可能解释为何我们的宇宙由物质构成的神秘"CP 破坏”的证据。这项工作提供了精确的数据，帮助物理学家调整他们对亚原子世界的理论。

技术摘要

技术摘要：$D^+ \to \pi^+\pi^0\pi^0$ 的振幅分析与分支比测量

问题与动机
粲介子的非轻子衰变为研究非微扰量子色动力学（QCD）效应和检验弱相互作用机制提供了关键平台。三体衰变 $D^+ \to \pi^+\pi^0\pi^0$ 在理论上预期主要由准二体道 $D^+ \to \rho(770)^+\pi^0$ 主导。尽管极点模型、拓扑图方法（TDA）和因子化辅助拓扑振幅（FAT）方法等理论框架为该主导模式的分支比（BF）提供了预测，但仍需精确的实验约束来完善这些模型。此外，虽然在中性 $D^0$ 衰变中已观测到直接 CP 破坏，但在 $D^+$ 衰变中尚未确立明显的 CP 不对称性。在单 Cabibbo 抑制的衰变 $D^+ \to \pi^+\pi^0\pi^0$ 中搜寻 CP 破坏对于探测潜在的新物理具有重要意义，特别是多体衰变中的 CP 不对称性可能源于中间共振态之间的长程干涉效应，即使全局不对称性很小，也可能揭示出局部的不对称性。

方法论
本分析利用 BESIII 探测器在质心能量 $\sqrt{s} = 3.773$ GeV 处收集的 $e^+e^-$ 对撞数据，对应积分亮度为 20.3 fb$^{-1}$。研究采用“标记”（tag）技术，其中 $D^+D^-$ 对通过 $\psi(3770)$ 衰变产生。一个 $D^-$ 介子在特定的强子衰变模式中被重建（单标记，ST），而信号 $D^+$ 则在 $D^+ \to \pi^+\pi^0\pi^0$ 模式中被重建（双标记，DT）。

• 事例选择：利用 BESIII 探测器组件（MDC、TOF、EMC）选择带电径迹和光子。候选事例利用束流约束质量（$M_{BC}$）和能量差（$\Delta E$）变量进行筛选。应用特定的否决条件以抑制背景，包括 $K_S^0$ 否决以排除 $D^+ \to K_S^0\pi^+$，以及运动学约束以排除来自 $D^0\bar{D}^0$ 对的错误分割背景。
• 信号提取：对 $M_{BC}^{tag}$ 与 $M_{BC}^{sig}$ 分布进行无分箱二维最大似然拟合，以提取信号产额，将其与组合背景及错误重建事例分离。
• 振幅分析：信号样本使用等标量模型（isobar model）进行无分箱振幅分析。总振幅是中间过程的相干叠加，包括 $D^+ \to \rho(770)^+\pi^0$、$D^+ \to \rho(1450)^+\pi^0$、$D^+ \to f_2(1270)\pi^+$ 以及 $\pi^0\pi^0$ S 波（通过 K 矩阵参数化描述）。分析考虑了两个全同 $\pi^0$ 介子之间的玻色对称性。探测效率和背景形状利用蒙特卡洛（MC）模拟进行建模，并针对径迹重建、粒子鉴别（PID）和 $\pi^0$ 重建中的数据-MC 差异进行了修正。
• CP 不对称性：对 $D^+$ 和 $D^-$ 样本进行联合拟合，允许等标量振幅的复系数和 K 矩阵 S 波的产生参数在电荷之间独立变化。
• 分支比：分支比通过双标记与单标记产额的比率计算得出，并针对探测效率和子衰变分支比（特别是 $\pi^0 \to \gamma\gamma$）进行了修正。

主要贡献与结果
本工作首次呈现了 $D^+ \to \pi^+\pi^0\pi^0$ 衰变的振幅分析。

1. 分支比测量：
   测得的总分支比为：
   $$\mathcal{B}(D^+ \to \pi^+\pi^0\pi^0) = (4.84 \pm 0.05_{\text{stat}} \pm 0.05_{\text{syst}}) \times 10^{-3}$$
   总分支比的 CP 不对称性测得为：
   $$A_{CP} = (-1.4 \pm 1.0_{\text{stat}} \pm 0.6_{\text{syst}})\%$$
   在积分产额中未观测到 CP 破坏的证据。

2. 振幅分析与拟合份额：
   分析证实该衰变主要由 $D^+ \to \rho(770)^+\pi^0$ 成分主导。测得的拟合份额（FF）和中间分支比如下：

   • $D^+ \to \rho(770)^+\pi^0$：FF = $63.5 \pm 2.0 \pm 1.2\%$；$\mathcal{B} = (3.08 \pm 0.10_{\text{stat}} \pm 0.07_{\text{syst}}) \times 10^{-3}$。
   • $D^+ \to \rho(1450)^+\pi^0$：FF = $5.2 \pm 0.8 \pm 0.7\%$。
   • $D^+ \to f_2(1270)\pi^+$：FF = $4.5 \pm 0.3 \pm 0.2\%$。
   • $D^+ \to (\pi^0\pi^0)_{S\text{-wave}}\pi^+$：FF = $11.6 \pm 0.9 \pm 0.5\%$。
     由于相消干涉，拟合份额之和小于 100%，而净干涉是相长的，其中 $f_2(1270)\pi^+ \times \rho(770)^+\pi^0$ 和 $\rho(1450)^+\pi^+ \times (\pi^0\pi^0)_{S\text{-wave}}\pi^+$ 贡献显著。

3. 中间态的 CP 不对称性：
   测量了特定中间振幅的 CP 不对称性。结果在误差范围内与零一致：

   • $A_{CP}(\rho(770)^+\pi^0) = +3.7 \pm 2.9 \pm 1.5\%$
   • $A_{CP}(\rho(1450)^+\pi^0) = -17.5 \pm 9.8 \pm 7.0\%$
   • $A_{CP}(f_2(1270)\pi^+) = -2.2 \pm 6.6 \pm 3.6\%$
   • $A_{CP}((\pi^0\pi^0)_{S\text{-wave}}\pi^+) = +10.2 \pm 6.6 \pm 3.8\%$

意义
本文声称，这些测量结果为粲介子衰变的理论模型提供了严格的约束。具体而言，将测得的主导道 $D^+ \to \rho(770)^+\pi^0$ 的分支比与极点模型、TDA 和 FAT 方法的预测进行了比较。作者指出，虽然极点模型的预测与其结果一致，但 FAT 的预测与测量值相差约 2.3 个标准差。此外，对中间态（包括复杂的 S 波结构）的拟合份额和相位的精确测定，为理解非轻子粲衰变中的拓扑振幅组成提供了新数据。尽管对 Dalitz 图中的局部效应具有敏感性，但该道中未观测到 CP 破坏，这为单 Cabibbo 抑制的 $D^+$ 衰变中潜在的新物理贡献设定了限制。