Search for new physics in $B \to K \pi \pi \gamma$ with Belle II data

本文利用 Belle II 数据讨论了$B \to K_{res} \gamma \to K \pi^+ \pi^- \gamma$衰变的振幅分析，这对于区分 CP 本征态与非 CP 本征态、进而通过时间依赖的 CP 不对称性测量来约束新物理至关重要。

要点

以下是用简单语言和创造性类比对该论文的解读。

宏观图景：在机器中猎捕幽灵

将物理学的标准模型想象成一本非常严格、编写完善的宇宙运行规则手册。几十年来，科学家们一直在核查这本手册，它通常完美契合。然而，他们怀疑可能有“幽灵”（新的、未知的粒子）隐藏在阴影中，以微妙的方式破坏规则。

本文讲述的是日本 Belle II 实验室（由 Belle II 合作组运营）进行的一项特定实验，旨在捕捉这些幽灵。他们正在观察一种特定类型的粒子衰变：一种被称为B 介子的重粒子分解为一个K 介子、两个π介子和一个光子（光）。

谜团：光的“手性”

在标准模型中，当 B 介子衰变为光子时，光子几乎总是“左旋”的（就像左旋螺丝）。如果科学家发现显著数量的“右旋”光子，那将是新物理的铁证。

为了测量这一点，他们观察CP 不对称性。这就像粒子与其镜像双胞胎（反粒子）之间的一场舞蹈。

• 如果舞蹈完全对称，规则就是标准的。
• 如果舞蹈失衡，就有某种新事物在推动舞者。

然而，存在一个问题。衰变的最终结果（K 介子 + 两个π介子）可以通过许多不同的“路径”或“路线”达到。其中一些路径是"CP 本征态”（完全对称的舞蹈），而另一些则是“非 CP 本征态”（混乱、不对称的舞蹈）。

类比：想象试图在拥挤的房间里听一段特定的小提琴独奏（信号）。但房间里充满了人们交谈、唱歌和鼓掌的声音（背景噪声和不同的衰变路径）。如果你只是听整个房间的声音，小提琴独奏就会被淹没。你需要将独奏与噪声分离开来，才能知道小提琴实际上有多响。

解决方案：“振幅分析”

论文解释说，为了发现新物理，他们必须执行振幅分析。这就像一位超级音响工程师，能够隔离管弦乐队中的每一件乐器，以确切了解它们是如何协同演奏的。

1. 管弦乐队：衰变并非沿一条直线发生。B 介子先转变为一个“共振态”（一种临时的重粒子），然后分解。存在许多可能的共振态（如 $K_1$、$K^*$ 等），每种都有不同的自旋和属性。
2. 干涉：这些不同的路径并非依次发生，而是同时发生并彼此“干涉”，就像池塘中的波浪相互碰撞。有时它们会增强信号；有时它们会抵消信号。
3. 目标：科学家们构建了一个复杂的数学模型（“衰变模型”），描述了所有可能的路径及其干涉方式。他们利用该模型计算出一个"稀释因子"。
   • 类比：如果混乱的舞蹈（非 CP 本征态）占人群的 90%，它们就会“稀释”对称舞蹈的信号。稀释因子能确切告诉他们信号被冲淡了多少，以便他们对此进行修正。

他们是如何做到的（实验室工作）

• 数据：他们使用了来自 SuperKEKB 对撞机的数据，该对撞机将电子和正电子相互撞击，产生数十亿个 B 介子。
• 过滤器：他们使用了一种称为sPlot的统计技巧，将真实的 B 介子衰变与背景噪声（看起来相似但并非如此随机碰撞）分离开来。
• 模拟：用于模拟这些事件的标准计算机程序不够好，因为它们无法理解不同路径之间复杂的“干涉”。因此，团队使用了一个名为AmpGen的新工具，来创建如果他们的模型正确，这些粒子应如何行为的逼真模拟。

目前的结果

本文展示的是初步工作。

• 他们已成功构建了数学模型，描述了 B 介子衰变为 K 介子和两个π介子的所有可能方式。
• 他们已在模拟数据上测试了该模型，并表明它能成功“拟合”数据，这意味着它能找出每条不同路径的强度和相位。
• 下一步：既然“引擎”已构建完成，他们需要对其进行微调（测试其稳健性），然后将其应用于 Belle II 收集的真实数据。

这为何重要

一旦他们将此模型应用于真实数据，他们就能计算出真实的 CP 不对称性，而不会受到混乱衰变路径造成的“稀释”。这将使他们能够精确测量光子的“左旋”与“右旋”。

如果结果偏离标准模型的预测，那将不仅仅是一个小误差；它将证明有一个新的重粒子潜伏在量子回路中，正在改变宇宙的法则。

简而言之：本文是关于构建一个复杂的数学滤波器，以在复杂的粒子衰变中将“信号”与“噪声”分离，从而使科学家最终能够看清宇宙是否正在打破其自身的规则。

技术摘要

基于所提供的预印本，以下是 Sahil Saha 代表 Belle II 合作组对论文“利用 Belle II 数据搜寻 $B \to K\pi\pi\gamma$ 中的新物理”的详细技术总结。

1. 问题陈述

本文旨在通过测量辐射衰变 $B^0 \to K_S^0 \pi^+ \pi^- \gamma$ 中的时间依赖 CP 不对称性，来搜寻超出标准模型（BSM）的新物理。

• 理论背景： 在标准模型（SM）中，味改变中性流（FCNC）跃迁 $b \to s\gamma$ 在树图阶被禁止，仅通过企鹅图发生。它受到 Glashow–Iliopoulos–Maiani (GIM) 机制的高度抑制。电弱圈图中的新重粒子可能会改变威尔逊系数（$C_7$ 和 $C'_7$），从而导致 CP 不对称性的偏差。
• 挑战： 测得的 CP 不对称性（$S_{B^0 \to K_S^0 \pi^+ \pi^- \gamma}$）是一个“有效”参数，由于末态中存在非 CP 本征态而被稀释。末态 $K_S^0 \pi^+ \pi^-$ 是各种中间共振态（奇异共振态，$K_{res}$）的叠加。其中一些共振态是 CP 本征态（例如 $\rho^0$），而另一些则不是。
• 目标： 为了在不假设特定场成分的情况下约束新物理（NP），必须从测得的有效不对称性中提取“纯”CP 不对称性（$S_{CP}$）。这需要计算一个稀释因子（$D$），该因子取决于不同衰变振幅之间的干涉。目前，若无对衰变动力学的完整振幅分析，无法确定该因子。

2. 方法论

本文概述了一种针对带电同位旋伙伴衰变 $B^+ \to K^+ \pi^+ \pi^- \gamma$ 的模型依赖型振幅分析形式体系的开发，该衰变在统计上更为有利，且假设不存在直接 CP 破坏。

A. 衰变模型构建

该分析将衰变建模为涉及三个步骤的相干叠加：

1. 产生： $B^+ \to K_{res}^+ \gamma$。
2. 共振态衰变： $K_{res}^+ \to K^+ \pi^+ \pi^-$，通过 various 中间态进行。
3. 末态： $K^+ \pi^+ \pi^-$ 态。

• 中间态： 模型包含了一个全面的奇异共振态（$K_{res}$）列表，其自旋 - 宇称 $J^P = 1^\pm, 2^\pm$。包括：
  • $K_1(1270)^+$, $K_1(1400)^+$
  • $K^*(1410)^+$, $K^*(1680)^+$
  • $K_2^*(1430)^+$, $K_2(1770)^+$
• 衰变链： 每个共振态衰变为 $K^* \pi$、$\rho K$、$\omega K$ 或 $K_2^* \pi$ 等组合。
• 振幅形式体系： 总概率定义为 $P_{total} = |\sum c_i A_i|^2$，其中 $c_i$ 为复权重，$A_i$ 为振幅。
  • 振幅 $A_i$ 被分解为自旋相关项（$S$）、线形（$T$）和 Blatt-Weisskopf 势垒因子（$BL$）。
  • 运动学： 分析利用五个可观测量：三个不变质量（$m_{K\pi\pi}$、$m_{K\pi}$ 等）和两个角度（$\cos\theta$、$\phi$），定义在 $K_{res}$ 静止系中。
  • 干涉： 模型使用相干和明确考虑了所有衰变模式之间的量子干涉。

B. 数据与模拟

• 数据来源： 来自 SuperKEKB 对撞机（$\sqrt{s} = 10.58$ GeV，$\Upsilon(4S)$ 共振态）$e^+e^-$ 碰撞收集的 Belle II 数据。
• 背景扣除： 使用 sPlot 技术提取拟合变量的背景扣除分布。区分变量为 $\Delta E$（能量差）。
• 蒙特卡洛生成： 标准生成器 EvtGen 不足，因为它缺乏 $K_{res}$ 模式之间正确的角结构和干涉。因此，合作组使用 AmpGen 生成 $B^+ \to K^+ \pi^+ \pi^- \gamma$ 衰变，随后通过全局 Belle II 模拟和重建软件进行处理。

C. 拟合策略

• 拟合确定每个中间衰变路径的复权重（幅值和相位）。
• 共振态极点质量和宽度固定为粒子数据组（PDG）的平均值。
• 拟合输出提供完整的振幅结构，用于通过对相空间积分（公式 1.5）计算稀释因子 $D$。

3. 主要贡献

• 首个形式体系开发： 这项工作首次提出了专门针对 Belle II 上 $B^+ \to K^+ \pi^+ \pi^- \gamma$ 衰变的振幅分析形式体系，涵盖了所有可能的中间奇异共振态。
• 工具集成： 成功将 AmpGen 生成器与 Belle II 模拟链集成，克服了 EvtGen 在角分布和干涉方面的局限性。
• 稀释因子计算： 该方法建立了计算稀释因子 $D$ 的必要框架，这是将测得的有效 CP 不对称性转化为对威尔逊系数（$C_7, C'_7$）约束的关键缺失环节。
• 综合模型： 模型包含了广泛的共振态谱（$K_1, K^*, K_2$）及其各种衰变链，确保了对 $K\pi\pi$ 末态的完整描述。

4. 结果

• 初步状态： 本文报告了初步进展以及在模拟数据上对拟合策略的验证。
• 模拟验证： 图 3.2 展示了使用完整衰变路径组合对生成器级模拟数据的成功拟合。拟合成功恢复了输入的振幅和相位，证实了模型解析复杂干涉模式的能力。
• 数据应用： 虽然形式体系已准备就绪，但对迄今为止注册的 Belle II 实际数据集的分析尚待进行。当前工作侧重于定义无偏拟合策略并测试其稳健性，然后再将其应用于真实数据。
• 当前测量背景： 本文引用了之前利用 Belle 和 Belle II 数据对有效 CP 参数 $S_{B^0 \to K_S^0 \pi^+ \pi^- \gamma} = -0.29 \pm 0.11 \pm 0.05$ 的测量，该结果有待稀释因子修正后才能被解释为新物理约束。

5. 意义

• 新物理灵敏度： 辐射 $B$ 衰变中 CP 不对称性的测量是探测新物理最敏感的探针之一。通过实现稀释因子的计算，这项工作使 Belle II 合作组能够将当前的实验限制转化为对 $b \to s\gamma$ 跃迁手征结构（特别是比率 $C'_7/C_7$）的严格约束。
• 模型无关性： 该方法允许在不指定新重粒子确切性质的情况下约束新物理，而是依赖于有效场论（EFT）框架。
• 未来前景： 该分析是 Belle II 下一代精密测量的先决条件。一旦应用于完整数据集，它将显著降低 CP 不对称性参数的不确定性，可能揭示出偏离标准模型的迹象，从而指向新粒子或新相互作用。