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

本論文は、オフアックスNuMIビームデータを用いたICARUS検出器による最初のニュートリノ断面積測定を提示し、複雑な原子核効果に対するニュートリノ事象発生器を検証するために、荷電流準弾性類似事象のフラックス平均微分断面積を様々な運動量変数について報告する。

要約

イリノイ州の地下深くに設置された巨大なハイテク水中カメラを想像してください。このカメラは「ICARUS」と呼ばれ、760 トンの液体アルゴン（凍結したネオンガス）で満たされています。その仕事は、宇宙から、そして近くの粒子加速器から絶えず地球に降り注ぐ「ゴーストのような粒子」、すなわち「ニュートリノ」の「写真」を撮影することです。

この論文は、この特定のカメラが初めてニュートリノがアルゴンとどのように相互作用するかを詳細に測定した際の成績表です。以下に、彼らが何を行い、何を発見したかを、簡単なアナロジーを用いて解説します。

設定：ゴーストとのビリヤードゲーム

ニュートリノは目に見えないゴーストのようなものです。彼らはめったに何とも衝突しません。もし衝突したとしても、それはゴーストのビリヤード玉が本物の玉に当たるようなものです。

• 源: 科学者たちは、フェルミ国立加速器研究所（巨大な粒子加速器）から発射されたニュートリノビームを使用しました。カメラはビームの中心からわずかにずれた位置にあるため、カメラに衝突するニュートリノは、特定の低いエネルギーの「速度」を持っています。
• 標的: 標的はカメラ内部の液体アルゴンです。
• 目的: 彼らは「準弾性散乱」と呼ばれる特定の衝突タイプを研究したかったのです。ニュートリノが原子の構成要素である陽子に衝突し、それを弾き飛ばす一方で、ニュートリノ自体が電子の重い従兄弟であるミューオンに変わる様子を想像してください。ここでの重要なルールは、パイオン（π中間子）は禁止ということです。もし衝突によってパイオン（別の種類の粒子）が生成された場合、それは別のゲームとなります。彼らが望んだのは、クリーンな「ノックアウト」ヒットだけでした。

課題：「核の霧」

この論文は、アルゴンの原子核が単一の陽子ではなく、陽子と中性子が混雑した部屋であることを示しているため、これらの衝突を研究することが難しいと説明しています。

• アナロジー: 暗く混雑した部屋で、ビリヤード玉が別の玉に衝突する様子を目撃しようとするのを想像してください。部屋にある他の玉が、移動している玉にぶつかり、その方向を変えたり、部屋から出る前に吸収したりする可能性があります。
• 問題点: 科学者たちは、この混雑した部屋の挙動を予測するための異なる「ルールブック」（コンピュータモデル）を持っています。あるモデルは玉同士が頻繁に跳ね返ると言い、他のモデルはくっつくと言っています。この不確実性が、宇宙の秘密を測定しようとする将来の実験にとって最大の頭痛の種となっています。

彼らが行ったこと：「写真アルバム」

研究者たちは、標的に衝突した 2.5 × 10²⁰個の陽子からのデータ（膨大な量のデータ）を収集しました。その後、コンピュータプログラムを使用して数百万の事象を分類し、以下の条件を満たす特定の「クリーン」な衝突を見つけました。

1. ミューオンが放出された。
2. 陽子が放出された。
3. その他何も（パイオンも、余分な破片も）放出されなかった。

彼らは、衝突後のビリヤード玉の測定を行うように、これらの衝突について 4 つの特定の事項を測定しました。

1. ミューオンの角度: ミューオンはどの方向に飛んだか？
2. ミューオンと陽子の間の角度: 彼らは互いからどのくらい離れて飛んだか？
3. 2 つの「不均衡」測定: 運動量は完全にバランスしていたか、それとも混雑した部屋（原子核）からの「蹴り」によって何かが乱されたか？

結果：ルールブックは一致するか？

測定値を取得した後、彼らはそれらを様々なコンピュータモデル（「ルールブック」）の予測と比較しました。

• 判決: 彼らが収集したデータは予測と一致しました。モデルが間違っているわけではありません。単に、どれが現実の最良の説明であるかを特定するには、まだ精度が十分ではないのです。
• 限界: この論文は、彼らの「不確実性予算」（測定値の誤差範囲）が現在広すぎることを述べています。それは、ぼやけたカメラで非常に似た 2 種類の青の微妙な違いを区別しようとするようなものです。彼らは青を見ることができますが、それが具体的にどの青であるかを確定的に言うことはまだできません。
• 主な犯人: 誤差の最大の源はニュートリノそのものではなく、検出器でした。カメラの感度と、粒子の「写真」を記録する方法が、最も大きな不確実性を引き起こしました。

結論

この論文は、この特定のカメラ（ICARUS）がアルゴン上の特定のニュートリノ相互作用を測定した初めての事例であるため、一つのマイルストーンです。

• なぜ重要か: 将来の実験（DUNE など）は、同様の検出器と標的を使用します。宇宙を理解するためには、ニュートリノがアルゴンに衝突したときにどのように振る舞うかを正確に知る必要があります。
• 教訓: 科学者たちは新しい「真の基準（グラウンド・トゥルース）」データセットを提供しました。現在のモデルはテストを合格していますが、データは異なる理論の中から勝者を選ぶにはまだ十分に精密ではありません。それを行うためには、より多くのデータと、カメラの仕組みについてのより鋭い理解が必要です。

要約すると：彼らはハイテクカメラを構築し、ニュートリノの衝突の 100 万枚の写真を撮影し、これらの粒子の振る舞いに関する現在の地図が概ね正しいことを確認しました。しかし、未来を航海するためには、より良い地図が必要です。

技術概要

論文「ICARUS におけるオフアックス NuMI ビームを用いたミューオン（反）ニュートリノの荷電流準弾性類似断面積の測定」の詳細な技術的要約を以下に示す。

1. 問題提起

現在および将来の長基線ニュートリノ振動実験（例：DUNE）の主な目標は、質量順序や CP 対称性の破れを含む振動パラメータを精密に測定することである。しかし、これらの測定は現在、ニュートリノ - 原子核相互作用のモデル化に支配される系統誤差によって制限されている。

• 課題: 原子核（この文脈では特にアルゴン）とのニュートリノ相互作用は複雑である。核効果（核子の初期状態の運動量分布や、粒子が原子核内で再散乱または吸収される最終状態相互作用（FSI）など）は、最終状態の運動学と粒子構成を歪ませる。
• ギャップ: 多くの事象生成器は、頂点相互作用とハドロン伝播を別々に扱う簡略化された因数分解スキームに依存している。理論モデルを検証し調整するために、特にこれらの核効果に敏感な運動量変数における、アルゴン標的に対する高精度な断面積測定が緊急に必要とされている。

2. 手法

実験設定

• 検出器: フェルミ国立加速器研究所のビーム発生源から 600 メートルの地点に設置された、760 トンの液体アルゴン時間投影室（LArTPC）であるICARUS T600検出器。
• ビーム: ニュートリノモードで運用されているNuMI（Main Injector におけるニュートリノ）ビーム。検出器はオフアックス（100 ミリラジアン）に位置しており、オンアックスビームと比較してニュートリノエネルギー分布がより低いエネルギー側へシフトし（ピークは 1–2 GeV 付近）、低エネルギー領域に集中している。
• データセット: ラン 1（2022 年 6 月–7 月）およびラン 2（2023 年 2 月–6 月）中に収集されたデータに相当する標的プロトン数（POT）$2.5 \times 10^{20}$。

シグナル定義と事象選択

• シグナルトポロジー: 以下のように定義される荷電流準弾性類似（CCQE-like）相互作用:
  • $\nu_\mu + \text{Ar} \to \mu^- + \text{陽子 (複数)} + \text{中性子 (複数)}$
  • 厳密な制約: 最終状態に荷電または中性のピオン、その他のメソンが存在しないこと。
  • 運動量カット:
    • ミューオン運動量 $p_\mu > 226$ MeV/c（閉じ込め/追跡を確保するため）。
    • 主要な陽子の運動量を 400 から 1000 MeV/c の間とする。
    • 横方向運動量不均衡（TKI）変数については、連続減速近似（CSDA）による精密な運動量再構成を確保するため、追加のカット $p_\mu < 0.8$ GeV/c を適用する。
• 背景抑制:
  • ピオン veto: 特定された荷電ピオンまたは光子誘起シャワー（$\pi^0$崩壊由来）を含む事象を棄却する。
  • サイドバンド制御: ピオン veto 条件を反転させた（タグ付き荷電ピオンを要求する）「ピオンサイドバンド」サンプルを定義し、背景モデル化をデータ駆動型で制約する。
  • 宇宙線背景: ビームスプリームと一致する宇宙線に対して、データ駆動型手法とシミュレーション（Corsika）を用いて軽減する。

シミュレーションと再構成

• 再構成: LArSoft フレームワーク内のPandoraパターン認識ソフトウェアを使用して、3D 軌跡とシャワーを再構成する。粒子識別（PID）は $dE/dx$ プロファイルと軌跡トポロジーに依存する。
• シミュレーション:
  • フラックス: 各種標的材料（C, Al, Fe）におけるハドロン生成およびオフアックス効果を考慮し、Geant4 および PPFX を用いてモデル化。
  • 相互作用: 参照モデルとしてGENIE事象生成器（バージョン AR23 20i CMC）を使用。QE にはバレンシアモデル、多核子効果を SuSAv2-MEC、FSI を hA2018 モデルとして組み込んでいる。
  • 検出器: 粒子の伝播、電離、ドリフトの完全な Geant4 シミュレーションを、信号誘起のための Wire-Cell Toolkit と結合。

解析フレームワーク

• 抽出: GUNDAM（Muon-neutrino 用のグローバル展開およびデータ解析）フレームワークを使用。
• フィッティング戦略: 背景を制約するために、シグナルおよびサイドバンドサンプルに対して同時フィッティングを行う。
• 変数: 4 つの運動量変数を測定する:
  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. オフアックスアルゴンデータ: 低エネルギー（MicroBooNE）および高エネルギー（ArgoNeuT）の測定を補完する、広範なエネルギー分布（数 GeV まで）を持つアルゴン標的におけるユニークなデータセットを提供する。
3. 多次元運動量解析: 単なる総断面積ではなく、核効果（初期状態の運動量および FSI）に対して特に敏感に設計された 4 つの変数におけるフラックス平均微分断面積を報告する。
4. 堅牢な背景制約手法: ピオン生成に対してデータ駆動型のサイドバンド手法を用いて、背景の系統誤差をデータから直接制約し、背景正規化における純粋なシミュレーションへの依存度を低減する。
5. モデル比較フレームワーク: 抽出されたデータを、代替 GENIE 設定（LQCD 形状因子、MINERvA 形状因子、hN2018 FSI）、NuWro、GiBUU、およびNeut（RDWIA モデル）を含む広範な理論予測と比較する体系的なアプローチを確立する。

4. 結果

• 断面積抽出: 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$ 分布において独特の形状（フェルミピークおよび高エネルギーテール）を示すが、データの 1 標準偏差の範囲内に留まっている。
    • 軸性形状因子の変動（LQCD 対 MINERvA フィット）および FSI モデルの変動（hA2018 対 hN2018）はわずかな逸脱を示すが、全体的な一致を著しく変えるものではない。
• 不確かさ予算: 現在の総不確かさは、フラックスや相互作用モデル化ではなく、検出器応答（カロリメトリー、陽子識別効率）によって支配されている。統計的不確かさも有意である。

5. 意義

• 将来実験への検証: この結果は、液体アルゴン TPC と同様のエネルギー分布を使用するDUNE（Deep Underground Neutrino Experiment）で用いられるニュートリノ相互作用モデルの調整に不可欠な入力情報を提供する。
• 核物理への洞察: $\delta p_T$ や $\delta \alpha_T$ などの変数を測定することで、標準的な QE モデルが正確に記述できない複雑な核環境（核子相関および FSI）を探求する。
• 限界と今後の課題: 著者は、現在の統計的および系統的不確かさは、特定の理論モデルを他よりも決定的に支持するには不十分であり（5% の有意水準でどのモデルも棄却されていない）、今後の改善には以下が必要であると指摘している:
  • より多くのデータ（反ニュートリノビーム極性を含む可能性）。
  • 陽子識別および軌跡再構成に関連する支配的な系統誤差を低減するための、洗練された検出器シミュレーション。

要約すると、この論文はアルゴンに対するミューオンニュートリノ相互作用のための厳密な基準を確立し、現在のモデルがデータと広範に一致している一方で、次世代ニュートリノ振動測定の精度を制限する核効果を完全に解明するには、より高い精度が必要であることを示している。