The $B^+ \to K^+ \nu \bar \nu$ decay as a QCD axion search: comparing reinterpretation approaches

技术摘要：将 $B^+ \to K^+ \nu\bar{\nu}$ 衰变作为 QCD 轴子搜索的研究：比较不同的重新诠释方法

问题陈述
针对 Belle II 关于 $B^+ \to K^+ \nu\bar{\nu}$ 衰变的两次近期独立的重新诠释分析，在将其重新诠释为寻找两体衰变 $B^+ \to K^+ a$（其中 $a$ 是轻质量不可见粒子，如 QCD 轴子）时，所得出的分支比 $\mathcal{B}(B^+ \to K^+ a)$ 的限制存在约四倍的差异。第一种方法 [1] 利用公开数据进行解析性的运动学变量重建，得出了目前最强的限制。第二种方法 [4] 依赖于已发表的模型无关似然函数 [5]，其结果显著较弱。作者研究了这种差异的起源，并假设其源于运动学变量空间的选择以及用于重新诠释的数据投影的粒度。

方法论
本文对应用于相同 Belle II 数据和模拟的两种不同统计框架进行了对比分析：

1. 细粒度 $q^2_{\text{rec}}$ 方法 (Ref. [1])： 该方法利用重建的二中微子不变质量 ($q^2_{\text{rec}}$) 分布，在 $[-1, 25] \text{ GeV}^2$ 的运动学范围内设有约 21 个分箱。它将背景归一化视为主要的系统不确定度（估计约为 $\sim 1\%$），并将真实的运动学映射到重建变量上。
2. 粗粒度 $q^2_{\text{rec}} \times \eta(\text{BDT}^2)$ 方法 (Ref. [4])： 该方法使用来自 Ref. [5] 的已发表似然函数，该函数将数据投影到由四个变换后的 BDT 输出分箱 ($\eta(\text{BDT}^2)$) 和三个粗略的 $q^2_{\text{rec}}$ 分箱（边界为 $\{-1, 4, 8, 25\} \text{ GeV}^2$）组成的二维空间中。该 BDT 变量是专门针对分离 $B^+ \to K^+ \nu\bar{\nu}$ 信号与背景而训练的。

作者通过五项数值测试来隔离变量粒度、BDT 轴相关性、系统不确定度和标准模型 (SM) 参数处理的影响：

• 测试 1 & 2： 在细粒度方法中改变背景归一化误差（从 1% 到 50%），以测试轴子限制和 $B^+ \to K^+ \nu\bar{\nu}$ 轮廓似然函数的稳健性。
• 测试 3： 直接运行 Ref. [4] 的已发表似然函数，通过改变对 $B^+ \to K^+ \nu\bar{\nu}$ 标度因子 ($\mu_{\nu\bar{\nu}}$) 和形状系统误差的约束，以建立基准。
• 测试 4： 构建一个专门的“混合”分析，将 Ref. [1] 的统计框架应用于 Ref. [4] 的粗粒度 $q^2_{\text{rec}} \times \eta(\text{BDT}^2)$ 分箱。这隔离了仅由运动学轴选择带来的影响。
• 测试 5： 在混合分析中检查 $B^+ \to K^+ \nu\bar{\nu}$ 的轮廓似然函数，以比较形状分离能力。

主要贡献与结果

• 分辨率是关键的品质因子： 灵敏度差异的主要驱动因素是 $q^2_{\text{rec}}$ 的分辨率。对于 QCD 轴子 ($m_a \simeq 0$)，信号在 $q^2=0$ 处是单色的，映射到 $q^2_{\text{rec}} \simeq 0$ 附近的结构化分布。Ref. [4] 似然函数的粗略 3 分箱结构将此信号合并到一个由 $B^+ \to K^+ \nu\bar{\nu}$ 背景主导的单一分箱中，从而稀释了信号。细粒度的 21 分箱结构则能解析出该信号结构，显著提高了灵敏度。
• BDT 轴的不充分性： 变量 $\eta(\text{BDT}^2)$ 是针对 $B^+ \to K^+ \nu\bar{\nu}$ 优化的，而非针对窄带 $B^+ \to K^+ a$ 信号。由于轴子产生单色卡昂，而 $\nu\bar{\nu}$ 产生连续谱，因此 BDT 轴几乎不具备判别能力，甚至可能抑制那些不符合 $\nu\bar{\nu}$ 拓扑结构的信号事件。
• 系统不确定度是保守的： 作者证明，Ref. [1] 方法中忽略次要形状系统误差并不会人为地收紧限制。包含这些系统误差（如 Ref. [4] 中所做）实际上会使限制降低约 20%，因为这允许拟合过程有更多的自由度去适应 $B^+ \to K^+ \nu\bar{\nu}$ 的形状，从而为轴子信号留下的空间更少。因此，Ref. [1] 的界限是保守的。
• 失去“双重探测”特性： 在细粒度 $q^2_{\text{rec}}$ 方法中，轴子限制在统计上与 $B^+ \to K^+ \nu\bar{\nu}$ 测量是解耦的（即“双重探测”特性）。而在粗粒度 $q^2_{\text{rec}} \times \eta(\text{BDT}^2)$ 空间中，取决于 $B^+ \to K^+ \nu\bar{\nu}$ 速率是固定还是允许浮动，轴子限制的变化幅度约为两倍，这表明其稳健性有所丧失。
• 量化灵敏度损失： 混合分析（测试 4）证实，仅运动学轴的变化本身就导致在低轴子质量下灵敏度下降约五倍。混合分析与 Ref. [4] 之间的其余差异归因于 Ref. [4] 中包含了完整的形状系统误差，这进一步降低了限制。

意义与主张
论文得出结论：对于窄带不可见信号的实验数据重新诠释，运动学变量空间的选择是决定性的。由针对特定信号拓扑（如 $B^+ \to K^+ \nu\bar{\nu}$）训练的 BDT 变量所主导的似然函数，在将测量结果重新诠释为具有显著不同形状的信号（如 $B^+ \to K^+ a$）时，其效用有限。

作者主张，实验合作组应在发布基于 BDT 的似然函数的同时，也发布物理变量空间（如细粒度 $q^2_{\text{rec}}$）的似然函数投影。这种做法可以最大限度地提高测量结果的可重新诠释性，以应对广泛的新物理情景，特别是涉及窄带共振或在运动学上与训练数据存在显著差异的信号。本文并非提出新的实验方案，而是呼吁修改数据发布策略，以促进更稳健的理论重新诠释。