Measurements of transverse kinematic imbalance in muon neutrino charged-current interactions on argon with a single proton and no pions in the final state with SBND

技術要約：SBNDによるアルゴン上でのミューオン・ニュートリノ荷電流相互作用における横方向運動学的不均衡の測定

問題と動機
加速器ニュートリノ実験は、検出されたニュートリノ相互作用の頻度とエネルギー分布から振動パラメータを抽出する。数GeV領域におけるこれらのスペクトルの解釈は、重い標的であるアルゴンに固有の競合する相互作用チャネルや核効果によって複雑化している。これらの核効果（核基底状態および最終状態相互作用（FSI）を含む）は、最終状態粒子の運動学を変化させ、長基線プログラム（T2K、NOvAなど）に重大な系統誤差をもたらし、DUNEやHyper-Kamiokandeのような次世代実験への課題を提起している。

$\nu_\mu$ CCQE-like（準弾性様）最終状態（単一のミューオンと単一の陽子からなり、パイ中間子が介在しないものとして定義される $\nu_\mu$CC1p0$\pi$）は、明快な実験的シグネチャを提供するが、核モデルを検証するためにはその運動学の精密な測定が必要である。MicroBooNEやICARUSによるアルゴン上での、またMiniBooNE、T2K、MINERv4による他の標的を用いた先行測定はベンチマークを提供してきたが、相互作用モデルをさらに制約するためには、アルゴンにおける高統計かつ多変数微分的な測定が必要である。

手法
本研究は、フェルミ研究所のShort-Baseline Near Detector (SBND) による、初のニュートリノ・アルゴン断面積測定を提示する。解析には、Booster Neutrino Beam (BNB) からの $8.8 \times 10^{19}$ プロトン・オン・ターゲット (POT) に相当するデータセットを利用している。平均ニュートリノエネルギーは $\langle E_\nu \rangle = 0.75$ GeV である。

• 検出器とシミュレーション： SBND検出器は112トンの液体アルゴン時比例室（LArTPC）である。シミュレーションおよび再構成はLArSoftフレームワークを用いて行われた。ニュートリノ相互作用はGENIE v3.4.0 (AR23 20i 00 000) を用いて、宇宙線はCORSIKAを用い、粒子の伝搬にはGeant4を用いてシミュレートされた。
• イベント選択： シグナルは、正確に1つのミューオン ($0.22 < p_\mu < 1.0$ GeV/c) と正確に1つの陽子 ($0.3 < p_p < 1.0$ GeV/c) を生成し、0.3 GeV/cを超える追加の陽子、70 MeV/cを超える荷電パイ中間子、および中性パイ中間子やより重いメソンが存在しない $\nu_\mu$ 荷電流相互作用として定義される。宇宙線除去、トポロジー選択（2トラック）、およびカロリメトリと運動学に基づく粒子識別（PID）を含む逐次的な選択手順により、純度91%で12,804個のシグナル候補が得られた。
• 運動学変数： 解析は、ニュートリノビームの横方向に垂直な平面における運動量の不均衡を定量化する横方向運動学的不均衡（TKI）変数に焦点を当てている：
  • $\delta p_T$: 横方向の欠損運動量ベクトル ($\vec{p}_\mu^T + \vec{p}_p^T$) の大きさ。
  • $\delta \alpha_T$: ミューオンの横方向運動量と横方向の欠損運動量との間の角度。
  • $\delta \phi_T$: ミューオンと陽子の横方向運動量間の方位角。
• アンフォールディングと系統誤差： 微分断面積は、二階微分正則化行列を用いたWiener-SVDアンフォールディングを用いて抽出された。系統誤差は、フラックス、相互作用モデル（GENIE）、二次相互作用（Geant4Reweight）、検出器応答、および正規化因子に関するシミュレーションの再重み付けによって評価された。モデルによる陽子運動学の予測を検証するために、データ駆動型の条件付き制約が適用された。

主要な貢献と結果
本論文は、アルゴン上での $\nu_\mu$CC1p0$\pi$ 相互作用のフラックス平均単一微分断面積を、$\delta p_T$、$\delta \alpha_T$、および $\delta \phi_T$ の関数として報告している。

• Data-MC比較： 選択されたデータは、評価された不確かさの範囲内で、正規化および形状の両方において、公称GENIE予測と一致している。フラックス平均断面積の全系統誤差は8.7%であり、ニュートリノフラックス（6.5%）によって支配されている。
• モデル識別： 高統計な測定は、イベントジェネレータの予測を強く識別する。比較されたモデル（GENIE AR23、GiBUU 2025、NEUT v6.1.4）のいずれも、3つのTKI変数を同時にうまく記述することには成功していない。
  • $\delta p_T$: すべてのモデルが満足な記述を提供できていない。GENIE AR23は、ピークおよび中間領域でデータを過小評価し、高エネルギー側では過大評価している。GENIE内のハートリー・フォック（HF）に基づくモデルは、この変数において改善した一致を示すが、他の変数の記述を悪化させる。
  • $\delta \alpha_T$: ほとんどのジェネレータは測定値を合理的に記述しているが、GiBUUは観測よりも急峻な勾配を予測し、HFモデルは勾配が緩やかすぎると予測している。
  • $\delta \phi_T$: ほとんどのジェネレータはこの変数を良好に記述しているが、GENIE AR23の軸性形式因子（axial form-factor）のバリアントは顕著な不一致を示している。
• 核物理学的知見： 不一致、特に $\delta p_T$ は、現在のジェネレータが最終状態の運動学全般の不一致というよりも、基礎となる核物理学（例：核子相関、FSI）のモデリングに欠陥があることを示唆している。結果は、異なるTKI変数が、現在のシミュレーションではまだ十分に捉えられていない相互作用モデルの補完的な側面を調査していることを示している。

意義
本研究は、ニュートリノ・アルゴン相互作用モデリングのための新しいベンチマークを確立するものである。現在のイベントジェネレータがTKI変数の全セットを同時に記述できないことを示すことで、本論文は核物理学のモデリングにおいて改善が必要な特定の領域を浮き彫りにしている。これらの結果は、振動測定における系統誤差を減少させ、将来の高精度実験の発見能力を高めるために極めて重要である。論文では、予想されるデータセットサイズ $10 \times 10^{20}$ POT を伴い、SBNDはニュートリノ・原子核相互作用モデリングの詳細を研ぎ澄ますために必要な、高統計かつ多変数微分的な測定を提供し続けると述べている。