Method to study $CP$ violation in $B_s^0\to K_S K^\pm \pi^\mp$ decays

想象一下，宇宙是一个巨大的、复杂的舞池，而被称为 B-介子（B-mesons） 的微小粒子就是其中的舞者。物理学家想要理解这些舞者的舞蹈规则，以观察它们是否符合“标准模型”（当前的物理规则手册），或者是否隐藏着我们尚未发现的秘密新动作（新物理学）。

这篇论文研究的是一种特定的、棘手的舞蹈动作： $B^0_s$ 介子 衰变为三个其他粒子（ $K^0_S$ , $K$ , 和 $\pi$ ）。

以下是作者所做工作的拆解，使用了简单的类比：

1. 谜团：“左手系”与“右手系”舞者

在粒子世界中，有一个概念叫做 CP 破坏（CP violation）。你可以这样理解：如果你拍摄一段粒子衰变的视频，然后倒着播放（或者看它的镜像），它看起来是否完全一样？

标准模型： 通常情况下，是的。舞蹈向前播放和向后播放看起来是一样的。
目标： 作者正在寻找一种在前向版本和后向版本看起来不同的舞蹈。发现这种差异是发现“新物理学”的有力证据。

2. 挑战：两阶段舞蹈程序

通常，物理学家一次只研究一种类型的舞蹈。但这种特定的衰变很特殊，因为它可以通过两种不同的方式（两个末态）发生，而这两个末态互为镜像。

问题： 为了捕捉到“CP 破坏”，你不能只观察一个舞者。你必须同时观察两个舞者，并在时间轴上逐帧对比他们的动作。
类比： 想象一下，试图分辨一名左撇子和一名右撇子网球运动员之间的细微差别。如果你只观察左撇子，你无法判断他是否在做一些奇怪的动作。你必须同时观察球场上的两名选手，对比他们的每一次挥拍、每一步移动，以及他们在场上停留的时间。

3. 方法：“达利兹图（Dalitz Plot）”地图

作者提出了一种分析数据的新方法，称为 达利兹图分析。

地图： 想象一张城市地图，地图上的每一个点都代表粒子飞散出的不同方式。
时间因素： 这不仅仅是一张静态地图；它是一部电影。作者正在创造一种方法，观察舞者如何随着时间在地图上移动。
标记（Tag）： 为了让这项工作奏效，他们需要知道哪个舞者最初是“左手系”版本，哪个最初是“右手系”版本。这被称为“味标记（flavor tagging）”。这就像是在比赛开始时，给一名舞者戴上红帽子，给另一名舞者戴上蓝帽子。

4. 实验：“虚假”舞池

由于他们目前还没有足够的来自 LHCb 实验的真实数据来进行这种复杂的分析，因此他们构建了一个模拟（称为“伪实验”）。

模拟： 他们创建了一个计算机程序，生成了 500 个“虚假”数据集，模拟未来 LHCb 将收集到的真实数据（特别是针对 Run 1、Run 2 和 Run 3）。
测试： 他们将这些虚假数据集输入到他们的新分析方法中，以测试该方法是否能成功找到他们植入代码中的隐藏“CP 破坏”信号。

5. 结果：它奏效了！

该论文声称，他们的新方法是可行的。

成功： 当他们运行“虚假”实验时，该方法成功恢复了隐藏的参数。它能够区分两种不同的舞蹈风格，并能高精度地测量“弱相位差”（CP 破坏的角度）。
精度： 他们发现，利用 LHCb 目前正在收集的数据（Run 1–3），我们可以非常精确地测量这个角度。如果等到未来有更多数据（Run 4–6），精度将会进一步提升。
工具： 他们已经将这种方法集成到了一个名为 Laura++ 的软件包中，其他科学家也可以使用它。

6. 为什么这很重要

新物理学： 如果真实的观测结果（当它到来时）显示出与标准模型预测不同的结果，这意味着有新的物理学隐藏在阴影之中。
蓝图： 这篇论文不仅仅研究一种衰变；它为如何研究任何具有两个镜像结局的复杂粒子衰变提供了一份蓝图（即一份食谱）。

总结

可以将这篇论文看作是一本针对非常困难的侦探游戏的训练手册。作者发明了一种新的方法，用于比较两个镜像对称的粒子舞蹈随时间的变化。他们在计算机模拟中测试了该方法，并证明了其有效性。现在，他们准备好将这种方法应用于来自大型强子对撞机的真实数据，看看宇宙是否还留有任何秘密的、打破规则的舞步等待被发现。

技术摘要：研究 $B^0_s \to K^0_S K^\pm \pi^\mp$ 衰变中 CP 破坏的方法

问题与动机
研究无粲（charmless） $B$ 介子衰变中的 CP 破坏对于测试标准模型（SM）并探测新物理（NP）至关重要。衰变道 $B^0_s \to K^0_S K^\pm \pi^\mp$ 特别值得关注。在标准模型中，当涉及中间 $K^{*0}$ 共振态时，该衰变由单一的 $t$ 圈图（penguin diagram）主导，这意味着不存在直接 CP 破坏，也不存在混合诱导的 CP 破坏（因为弱相位差消失）。然而，高阶 $c$ 圈和 $u$ 圈图的贡献或新物理可能会改变这一预期。相反，当涉及 $K^{*\pm}$ 共振态时，树级图会有贡献，从而可能允许 CP 破坏并对 CKM 角 $\gamma$ 具有敏感性。

为了获取这些衰变中的全套 CP 破坏可观测物理量，需要进行一种带有味标签（flavor-tagged）、随时间变化的 Dalitz 图分析。虽然此类分析已在单自共轭末态（例如 $B^0 \to K^0_S \pi^+ \pi^-$ ）以及四体衰变中得到应用，但从未针对两个电荷共轭末态（ $K^0_S K^+ \pi^-$ 和 $K^0_S K^- \pi^+$ ）进行过同时进行的随时间依赖的分析。目前的 LHCb Run 1 和 Run 2 数据集缺乏进行味标签分析所需的统计量，但 Run 3 数据集（及未来的运行）预计将提供足够的产额。

方法论
作者提出了一种形式化方法和手段，用于对两个末态 $f \equiv K^0_S K^+ \pi^-$ 和 $\bar{f} \equiv K^0_S K^- \pi^+$ 进行同步的、带味标签的、随时间变化的 Dalitz 图分析。

形式化： 最初产生的 $B^0_s$ 和 $\bar{B}^0_s$ 介子衰变为 $f$ 和 $\bar{f}$ 的微分衰变宽度被表示为关于衰变时间 $t$ 和相空间变量 $\Phi_3$ 的函数。这些方程包含了 $B^0_s$ - $\bar{B}^0_s$ 混合参数（ $\Delta m_s, \Delta \Gamma_s, q/p$ ）以及衰变振幅（ $A_f, \bar{A}_f, A_{\bar{f}}, \bar{A}_{\bar{f}}$ ）。CP 破坏表现为振幅大小的差异（直接 CPV）以及混合与衰变之间的干涉（混合诱导 CPV），其特征为有效弱相位差 $\phi^{\text{eff}}_s$ 。
振幅模型： 衰变振幅使用等量模型（isobar formalism）构建，通过叠加中间共振态的贡献来实现。在可行性研究中，采用了仅包含 $K^*(892)$ 共振态（ $K^{*0}, K^{*+}, K^{*-}$ ）的简化模型。复系数（ $c_i$ ）被分解为 CP 守恒部分（ $c$ ）和 CP 破坏部分（ $\Delta c$ ）。
伪实验设置： 该方法在 Laura++ Dalitz 图分析软件包中实现。伪实验是基于 LHCb Run 1–3 数据集（总产额约 130,000，有效标签产额约 6,500）生成的。
拟合策略： 在 500 个伪实验组成的系综上进行非分箱最大似然拟合，以解决拟合歧义，应用了特定约束：
- 参考振幅系数 $c_-$ 被固定为 1。
- $\Delta c_-$ 的虚部被设为 0。
- 对 $A_{\bar{f}}$ 与 $\bar{A}_{\bar{f}}$ 之间的相对相位应用了约束。
- $\phi^{\text{eff}}_s$ 的解被限制在输入值的 $\pm \pi/3$ 范围内，以解决 $\pi$ 歧义。
- 每个伪实验执行 100 次带有随机起始点的拟合，以确保收敛至全局最小值。

核心贡献

新型分析框架： 本文首次概述了对两个电荷共轭末态进行同步、带味标签、随时间变化的 Dalitz 图分析的方法。
实现： 该方法已在 Laura++ 包中实现，为这类测量建立了一个基准框架。
敏感性研究： 作者评估了在各种模型配置下（改变振幅大小、相位，并引入 $K^*$ 振幅中的 CP 破坏），对有效弱相位 $\phi^{\text{eff}}_s$ 和振幅系数中 CP 破坏相位的敏感性。

结果

可行性： 伪实验表明，输入参数（包括 $\phi^{\text{eff}}_s$ 和振幅系数）可以被可靠地确定。
$\phi^{\text{eff}}_s$ 的精度： 对于 LHCb Run 1–3 数据集， $\phi^{\text{eff}}_s$ 可达到约 26 mrad 的统计精度。这与目前使用 $B^0_s \to J/\psi \phi$ 的 LHCb/CMS 测量值（22 mrad）相当，且显著优于 $B^0_s \to \phi \phi$ （75 mrad）。
振幅相位的精度： 对 CP 破坏振幅分量相位（ $\arg(\Delta c_0), \arg(\Delta \bar{c}_0)$ ）的敏感度较低，不确定度约为 200 mrad。这归因于参考带电 $K^*$ 振幅与中性 $K^*$ 振幅之间有限的相空间重叠。
未来预测： 将统计量扩展到预期的 LHCb Run 4（60 fb $^{-1}$ ）和 Run 6（300 fb $^{-1}$ ）数据集， $\phi^{\text{eff}}_s$ 的投影不确定度将分别提高到 17 mrad 和 8.6 mrad，遵循随亮度平方根倒数变化的规律。
模型依赖性： $\phi^{\text{eff}}_s$ 的精度显示出对底层振幅模型的依赖性。在准两体极限下，当强相位差接近 $\pm \pi/2$ 时，敏感度会下降。然而，作者指出，包含额外共振态（如 $K^*_0(1430)$ ）的更现实的模型可能会降低这种依赖性。

意义与主张
本文声称，所提出的方法是可行的，并为利用 LHCb Run 3 数据集测量 $B^0_s \to K^0_S K^\pm \pi^\mp$ 衰变中的 CP 破坏提供了路径。该研究为未来的数据分析建立了必要的基准，这些分析有望在 Run 4 及以后阶段实现精度的进一步提升。作者强调，尽管目前的实验使用了简化模型并忽略了实验效应（背景、分辨率、标签缺陷），但所展示的敏感度表明，可以实现令人感兴趣的精度水平。该方法也被指出可以扩展到具有多个末态的其他多体衰变。对应的 $B^0 \to K^0_S K^\pm \pi^\mp$ 衰变被强调为未来的目标，其中 U-spin 对称性可以联系 $B^0_s$ 和 $B^0$ 系统之间的振幅，从而控制强子学不确定性。

Method to study $CP$ violation in Bs0→KSK±π∓B_s^0\to K_S K^\pm \pi^\mpBs0​→KS​K±π∓ decays