Measurement of muon (anti-)neutrino charged-current quasielastic-like cross section using off-axis NuMI beam at ICARUS

本文利用离轴 NuMI 束流数据，首次报告了 ICARUS 探测器的中微子截面测量结果，给出了不同运动学变量下类电荷电流准弹性事例的通量平均微分截面，以检验中微子事例生成器对复杂核效应的描述能力。

要点

想象一下，在伊利诺伊州深处地下，放置着一台巨大的高科技水下相机。这台名为ICARUS的相机内装有 760 吨液态氩（冷冻的氖气）。它的工作是拍摄来自太空以及附近粒子加速器不断轰击地球的幽灵般粒子——中微子的“照片”。

本文报告了这台特定相机首次成功详细测量中微子与氩相互作用的情况。以下是他们所做工作及发现的具体分解，并辅以简单的类比。

设置：一场与幽灵的台球赛

中微子就像看不见的幽灵。它们极少与任何事物发生碰撞。当它们确实击中某物时，就像幽灵般的台球击中了真实的台球。

• 源头：科学家使用了从费米实验室（Fermilab，一个巨型粒子加速器）发射的中微子束。由于相机位于略微偏离中心的位置（并非正对束流中心），击中相机的中微子具有特定的、较低能量的“速度”。
• 目标：目标就是相机内的液态氩。
• 目标：他们希望研究一种特定的碰撞，称为“准弹性”（Quasi-Elastic）。想象一下，一个中微子击中一个质子（原子的基本组成部分）并将其撞出，同时中微子转变为一个μ子（电子的重表亲）。这里的关键规则是：不允许产生π介子。如果碰撞产生了π介子（另一种类型的粒子），那就是另一场游戏了。他们只想要干净的“击出”碰撞。

挑战：“原子核迷雾”

本文解释说，研究这些碰撞之所以困难，是因为氩原子核不仅仅是一个单独的质子；它是一个拥挤的质子和中子房间。

• 类比：想象试图在黑暗拥挤的房间里看到一颗台球击中另一颗台球。房间里的其他球可能会撞击移动的球，改变其方向，或者在它甚至离开房间之前将其吸收。
• 问题：科学家们有不同的“规则手册”（计算机模型）来预测这个拥挤房间的行为。有些模型认为球之间会频繁反弹；另一些模型则认为它们会粘在一起。这种不确定性是未来试图测量宇宙秘密的实验面临的最大难题。

他们做了什么：“相册”

研究人员收集了 2.5 × 10²⁰个质子撞击靶标的数据（海量数据）。然后，他们使用计算机程序筛选数百万个事件，以寻找特定的“干净”碰撞，其中：

1. 有一个μ子飞出。
2. 有一个质子飞出。
3. 没有其他东西（没有π介子，没有额外碎片）飞出。

他们测量了这些碰撞的四个具体方面，就像在撞击后测量台球一样：

1. μ子的角度：μ子朝哪个方向飞？
2. μ子与质子之间的角度：它们彼此飞离了多远？
3. 两个“不平衡”测量值：动量是否完美平衡，还是来自拥挤房间（原子核）的“踢击”使情况发生了偏差？

结果：规则手册匹配吗？

一旦获得测量数据，他们便将其与各种计算机模型（“规则手册”）的预测进行了比较。

• 裁决：他们收集的数据符合预测。这些模型并没有错；只是目前还不够精确，无法判断哪一个是现实的最佳描述。
• 局限性：论文指出，他们的“不确定度预算”（测量中的误差范围）目前太宽。这就像试图用模糊的相机分辨两种非常相似的蓝色色调。他们能看到蓝色，但还无法确切地说出具体是哪种色调。
• 主要罪魁祸首：最大的误差来源并非中微子本身，而是探测器。相机的灵敏度及其记录粒子“照片”的方式引入了最大的不确定性。

结论

这篇论文是一个里程碑，因为这是这台特定相机（ICARUS）首次测量这些特定的中微子与氩的相互作用。

• 为何重要：未来的实验（如 DUNE）将使用类似的探测器和靶标。为了理解宇宙，他们需要确切知道中微子撞击氩时的行为。
• 要点：科学家们提供了一组新的“地面真实”数据。虽然目前的模型通过了测试，但数据还不够精确，无法在不同理论中选出优胜者。要做到这一点，他们需要更多的数据，以及对相机工作原理更清晰的理解。

简而言之：他们建造了一台高科技相机，拍摄了中微子撞击的百万张照片，并确认了我们目前关于这些粒子行为的大致地图是正确的，但我们需要更精确的地图来导航未来。

技术摘要

以下是论文《利用 ICARUS 在离轴 NuMI 束流下测量缪子（反）中微子带电流准弹性类截面》的详细技术总结。

1. 问题陈述

当前及未来长基线中微子振荡实验（如 DUNE）的主要目标是精确测量振荡参数，包括质量顺序和 CP 破坏。然而，这些测量目前受到系统不确定性的限制，而这些不确定性主要由中微子 - 原子核相互作用的建模主导。

• 挑战： 中微子与原子核（在此背景下特指氩核）的相互作用非常复杂。核效应——例如核子的初始态动量分布，以及粒子在原子核内发生再散射或被吸收的末态相互作用（FSI）——会扭曲末态的运动学特征和粒子成分。
• 差距： 许多事例生成器依赖于简化的因子化方案，将顶点相互作用与强子传播分开处理。迫切需要针对氩靶进行高精度的截面测量，特别是针对对这些核效应敏感的运动学变量，以验证和调节理论模型。

2. 方法论

实验设置

• 探测器： ICARUS T600 探测器，一个 760 吨的液态氩时间投影室（LArTPC），位于费米实验室束流原点 600 米处。
• 束流： 运行在中微子模式下的 NuMI（主注入器中微子） 束流。探测器位于离轴位置（100 毫弧度），导致中微子能谱向较低能量偏移（峰值在 1–2 GeV 左右），与在轴束流相比有所不同。
• 数据集： 收集于第 1 次运行（2022 年 6 月–7 月）和第 2 次运行（2023 年 2 月–6 月）的数据，对应 $2.5 \times 10^{20}$ 个靶上质子（POT）。

信号定义与事例选择

• 信号拓扑： 带电流准弹性类（CCQE-like）相互作用，定义为：
  • $\nu_\mu + \text{Ar} \to \mu^- + \text{质子(s)} + \text{中子(s)}$
  • 严格约束： 末态无带电或中性π介子或其他介子。
  • 运动学截断：
    • 缪子动量 $p_\mu > 226$ MeV/c（以确保包含/径迹重建）。
    • 领头质子动量介于 400 至 1000 MeV/c 之间。
    • 对于横向运动学不平衡（TKI）变量，额外施加 $p_\mu < 0.8$ GeV/c 的截断，以确保通过连续减速近似（CSDA）进行精确的动量重建。
• 背景抑制：
  • π介子否决： 拒绝包含已识别的带电π介子或光子诱导簇射（来自$\pi^0$衰变）的事例。
  • 侧带控制： 通过反转π介子否决条件（要求标记带电π介子）来定义"π介子侧带”样本，以数据驱动的方式约束背景建模。
  • 宇宙线背景： 利用数据驱动方法和模拟（Corsika）来缓解与束流 spill 重合的宇宙线背景。

模拟与重建

• 重建： 使用 LArSoft 框架内的 Pandora 模式识别软件重建三维径迹和簇射。粒子识别（PID）依赖于 $dE/dx$ 分布和径迹拓扑。
• 模拟：
  • 通量： 使用 Geant4 和 PPFX 建模，考虑离轴效应以及在各种靶材料（C、Al、Fe）上的强子产生。
  • 相互作用： 使用 GENIE 事例生成器（版本 AR23 20i CMC）作为参考模型。它结合了用于 QE 的 Valencia 模型、用于多核子效应的 SuSAv2-MEC，以及用于 FSI 的 hA2018 模型。
  • 探测器： 对粒子传播、电离和漂移进行完整的 Geant4 模拟，并结合 Wire-Cell Toolkit 进行信号感应。

分析框架

• 提取： 使用 GUNDAM（缪子中微子全局展开与数据分析）框架。
• 拟合策略： 对信号和侧带样本进行同时拟合以约束背景。
• 变量： 测量四个运动学变量：
  1. $\cos \theta_\mu$：缪子相对于中微子方向的角度。
  2. $\cos \theta_{\mu,p}$：缪子与领头质子之间的张角。
  3. $\delta p_T$：横向动量不平衡。
  4. $\delta \alpha_T$：横向角度不平衡。
• 系统误差： 构建了一个全面的协方差矩阵，涵盖通量、探测器响应（增益、噪声、电子寿命、质子识别效率）、相互作用模型（形状因子、MEC、FSI）以及粒子传播（Geant4 重加权）。

3. 主要贡献

1. 首个 ICARUS 截面测量： 这是利用 NuMI 数据从费米实验室 ICARUS 探测器发表的首个中微子截面测量结果。
2. 离轴氩数据： 提供了一个独特的氩靶数据集，具有宽能谱（高达几 GeV），补充了低能（MicroBooNE）和更高能（ArgoNeuT）的测量。
3. 多维运动学分析： 报告了四个运动学变量的通量平均微分截面，这些变量是专门针对其对核效应（初始态动量和 FSI）的敏感性而选择的，而不仅仅是总截面。
4. 稳健的背景约束方法： 利用数据驱动的侧带方法处理π介子产生，直接从数据中约束背景不确定性，减少了对纯模拟进行背景归一化的依赖。
5. 模型比较框架： 系统地将提取的数据与广泛的理论预测进行比较，包括替代的 GENIE 配置（LQCD 形状因子、MINERvA 形状因子、hN2018 FSI）、NuWro、GiBUU 和 Neut（RDWIA 模型）。

4. 结果

• 截面提取： 论文展示了四个运动学变量的提取通量平均微分截面。信号主要由准弹性（QE）相互作用主导（约 60%），其次是多核子（MEC，约 20%）和共振（RES，约 15%）贡献。
• 模型比较：
  • 数据与参考 GENIE 模型及其他生成器进行了比较。
  • 良好一致性： 总体而言，数据与预测一致。数据与各种模型比较的 p 值，角度变量范围为 0.178 至 0.988，TKI 变量范围为 0.318 至 0.981。
  • 具体发现：
    • Neut 模型（RDWIA）在角度变量上显示出最低的一致性（$p=0.178$），但在 TKI 变量上吻合良好。
    • NuWro 和 GiBUU 在 $\delta p_T$ 分布中显示出独特的形状（费米峰和高能尾），但仍处于数据的一个标准差范围内。
    • 轴矢量形状因子的变化（LQCD 与 MINERvA 拟合）以及 FSI 模型（hA2018 与 hN2018）的变化显示出微小偏差，但并未显著改变整体的一致性。
• 不确定度预算： 总不确定度目前主要由探测器响应（量热法、质子识别效率）主导，而非通量或相互作用建模。统计不确定度也较为显著。

5. 意义

• 未来实验的验证： 这些结果为调节用于 DUNE（深地中微子实验）的中微子相互作用模型提供了必要的输入，DUNE 也将使用液态氩 TPC 和相似的能谱。
• 核物理见解： 通过测量 $\delta p_T$ 和 $\delta \alpha_T$ 等变量，该分析探测了复杂的核环境（核子关联和 FSI），而标准 QE 模型往往无法准确描述这些环境。
• 局限性与未来工作： 作者指出，当前的统计和系统不确定度还不足以明确支持某一个特定的理论模型（没有任何模型在 5% 显著性水平上被拒绝）。未来的改进需要：
  • 更多数据（可能包括反中微子束流极性）。
  • 改进的探测器模拟，以减少与质子识别和径迹重建相关的主导系统不确定度。

总之，本文建立了缪子中微子在氩上相互作用的严格基准，表明虽然当前模型与数据大体一致，但需要更高的精度才能完全解决限制下一代中微子振荡测量精度的核效应。