Measurement of $\pi^0$ Production in $\bar{\nu}_{\mu}$ Charged-Current Interactions in the NOvA Near Detector

NOvA 近傍検出器を用いた高統計量解析により、平均エネルギー 2GeV の反ミューニュートリノ誘起チャージドカレント中性パイオン生成過程の微分断面積を測定し、GENIE モデルとの整合性を確認するとともに、他のモデルが$\Delta$(1232) 共鳴領域で断面積を過小評価している可能性を示唆した。

要約

素粒子の「お菓子作り」実験：NOvA による新しい発見

この論文は、アメリカのフェルミ国立加速器研究所（Fermilab）で行われている「NOvA（ノバ）」という実験チームによる、非常に精密な「お菓子作り」の記録です。

ここでは、専門用語を避け、誰でもイメージしやすい**「巨大なお菓子工場」と「お菓子のレシピ」**というメタファーを使って、この研究の何がすごいのかを説明します。

---

1. 実験の舞台：巨大なお菓子工場（NOvA 実験）

まず、実験が行われている場所を想像してください。
アメリカ中西部にある巨大な**「ニュートリノお菓子工場」**です。

• 原料（ニュートリノ）： 工場の奥では、120GeV という強力なエネルギーを持った「ニュートリノ」という、幽霊のように物質をすり抜ける不思議な粒子のビーム（光の束）が作られています。今回は、そのビームの向きを調整して、「反ニュートリノ（Anti-neutrino）」という、少し性格の違う原料を使っています。
• 作業台（検出器）： この原料ビームを、1 キロメートルほど離れた場所にある「近接検出器（Near Detector）」という巨大な水槽にぶつけます。この水槽は、液体が入ったプラスチックの箱が何千枚も積み重なったもので、中に入ると「光」が走る様子がカメラで捉えられます。

2. 実験の目的：「ゼロ」の粒子を捕まえる

今回の実験のゴールは、**「反ニュートリノが炭素（プラスチックの材料）にぶつかったとき、どうやって『中性パイオン（π0）』というお菓子を作るのか」**を詳しく調べることです。

• なぜ重要なのか？
  宇宙の謎を解く「ニュートリノ振動」という現象を調べるには、ニュートリノがどう反応するかを正確に知る必要があります。しかし、この「中性パイオン」は、電子ニュートリノの信号と間違えられやすい「トリックな犯人」のような存在です。
  もし、このお菓子がどうやって作られるかの「レシピ（理論）」が間違っていると、宇宙の謎を解く計算も全部間違ったものになってしまいます。

3. 実験のプロセス：お菓子の形を計る

実験チームは、反ニュートリノがぶつかった瞬間に何が起こったかを、まるで**「お菓子の断面図」を詳しく描く**ように測定しました。

• お菓子の形（運動量と角度）： できた「中性パイオン」と、一緒に飛び出した「ミューオン（もう一つのお菓子）」が、どのくらいの速さで、どの方向へ飛んでいったかを測ります。
• エネルギーの量（Q2 と W）： ぶつかり具合の強さや、できたお菓子の重さ（質量）を計算します。

これまでに、この「反ニュートリノ」を使った実験は世界で 1 回しかありませんでした。しかし、今回の NOvA 実験は、その 6 倍ものデータを集めることに成功しました。まるで、1 回しか見たことのないお菓子を、6 回分も詳しく観察して、レシピを完成させたようなものです。

4. 発見：既存のレシピは少し間違っていた！

実験結果を、世界中の科学者が使っている「お菓子作りのレシピ本（シミュレーションモデル）」と比べてみました。

• GENIE（ジェニー）というレシピ本：
  最もよく使われているレシピ本「GENIE」の予測は、全体的にはデータとよく合っていました。これは「おおむね正しいレシピ」を持っているということです。
• しかし、ある部分でズレが：
  問題が見つかったのは、「デルタ（Δ）という特別な型」を使って作るお菓子の部分です。
  実験データによると、この「デルタ型」のお菓子は、レシピ本が予想していたよりももっとたくさん作られていました。
  他の 2 つのレシピ本（NuWro や NEUT）は、この部分でデータを大きく見逃しており、「もっと少ないはずだ」と言っていたのですが、実際にはもっと多かったのです。

簡単な例え：
「このお菓子を作ると、10 個できるはずだ」というレシピがあったのに、実験では「12 個できていた！」という結果が出たのです。特に、お菓子の形が少し丸い（運動量が低い）部分で、その差がはっきりしました。

5. なぜこれがすごいのか？

この発見は、単に「お菓子の数」が違ったという話ではありません。

• 宇宙の謎を解く鍵：
  将来、ニュートリノを使って「宇宙の始まり」や「物質の成り立ち」を解明する実験では、この「中性パイオン」の反応を正確に計算する必要があります。今回の実験で「デルタ型」の作り方が少し違っていたことがわかったことで、今後の実験の計算がより正確になります。
• 史上最も精密な測定：
  これは、これまで行われた中で、最も正確で詳細な「反ニュートリノによる中性パイオン生成」の測定です。

まとめ

この論文は、**「幽霊のような粒子（ニュートリノ）を使って、お菓子（中性パイオン）を作る実験を行い、世界中のレシピ本（理論モデル）が、ある特定の作り方（デルタ共鳴）の数を少し低估していたことを、6 倍のデータ量で証明した」**という画期的な成果です。

これにより、科学者たちはより正確な「宇宙のレシピ」を完成させるための、重要なピースを手に入れたのです。

技術概要

NOvA 実験における反ミューニュートリノ誘起中性パイオンの生成測定：技術的サマリー

本論文は、フェルミ国立加速器研究所（Fermilab）の NOvA 実験のニア検出器（Near Detector）を用いて、炭化水素（CH）標的における反ミューニュートリノ（$\bar{\nu}_\mu$）の荷電流（CC）相互作用による中性パイオン（$\pi^0$）生成を高精度で測定した結果を報告するものです。平均ニュートリノエネルギー 2 GeV における微分断面積が測定され、各種ニュートリノ相互作用モデルとの比較が行われました。

以下に、問題提起、手法、主要な貢献、結果、および意義について詳細をまとめます。

1. 問題提起と背景

• ニュートリノ振動実験における課題: 将来のニュートリノ振動測定（特に電子ニュートリノの出現事象の探索）において、ニュートリノおよび反ニュートリノの原子核標的に対する断面積の精密な知識は不可欠です。断面積の不確実性、特に原子核環境に起因する効果のモデル化の誤差は、ニュートリノエネルギーの再構成精度を制限し、最終的に振動パラメータの精度を低下させる要因となります。
• 中性パイオン（$\pi^0$）の重要性: 中性カレント（NC）相互作用で生成される$\pi^0$は、$\gamma\gamma$崩壊を起こし、電子ニュートリノの信号事象と誤認される主要な背景要因となります。NC $\pi^0$生成を直接測定することは困難ですが、荷電流（CC）相互作用における$\pi^0$生成の断面積を精密に測定することで、関連する NC 過程のモデル制約を得ることができます。
• 既存データの不足: 中性子（$\nu$）誘起の$\pi^0$生成については高精度な測定例がありますが、反ニュートリノ（$\bar{\nu}$）誘起の CC $\pi^0$生成（$\bar{\nu}_\mu + CH \to \mu^+ + \pi^0 + X$）に関する測定は、平均エネルギー 3.6 GeV での 1 件のみ（MINERvA 実験）でした。本研究は、これに比べて統計量が 6 倍多く、平均エネルギーが約 2 GeV の新しい測定を提供します。

2. 手法と実験設定

• 実験装置:
  • NuMI ビーム: Fermilab のメインインジェクターから抽出された 120 GeV プロトンがグラファイト標的に衝突し、生成されたパイオン・カオンを磁気ホーンで集束・偏向させます。ホーンの極性を反転させることで、主に反ニュートリノビームを生成します。
  • ニア検出器（ND）: ビーム軸から約 14.6 mrad ずれた位置に設置された液体シンチレータ検出器です。PVC セルに鉱物油とプスコクメンを充填した構造で、全質量 290 トン（シンチレータ 130 トン）。ミューオン捕獲器（Muon Catcher）を備え、最大 2.5 GeV/c のミューオンを捕捉可能です。
• シミュレーション:
  • GENIE v3.0.6: ニュートリノ相互作用のシミュレーションに使用。G18_10j_00_000 モデル構成を採用し、Valencia 模型（QE）、Berger-Sehgal 模型（共鳴）、Bodek-Yang 模型（DIS）を組み合わせています。
  • FSI（最終状態相互作用）: 原子核内でのハドロン再相互作用を半古典的カスケード模型（hN）でシミュレート。NOvA 独自のチューニング（"NOvA Tune v2"）を適用し、パイオンの運命（吸収、電荷交換など）を調整しています。
• イベント選択と再構成:
  • 信号定義: 最終状態にミューオンと少なくとも 2 つの光子候補（$\pi^0 \to \gamma\gamma$ からの再構成）を持つ事象。
  • 選別基準: ミューオンの運動量 $0.5 \le p_\mu < 2.5$GeV/c、角度$\theta_\mu < 60^\circ$。
  • 機械学習の活用:
    • MuonID: Boosted Decision Tree (BDT) を用いたミューオン識別。
    • EMscore: 畳み込みニューラルネットワーク（CNN）を用いて、光子（EM-like）と荷電ハドロン（非 EM-like）を識別。
  • 不純物除去: 再構成された$\pi^0$の質量ウィンドウ（$90 \le m_{\gamma\gamma} < 190$MeV/$c^2$）を適用し、背景を抑制。
• 統計解析（テンプレートフィット）:
  • 信号領域（CCN$\pi^0$）と、背景支配のサイドバンド領域（CC0$\pi^0$、NC 相互作用など）を定義。
  • 4 つのカテゴリ（信号、CC0$\pi^0$、NC、その他）のテンプレート分布を用いた同時フィットにより、各運動量ビンにおける信号と背景の寄与を抽出。
• 展開（Unfolding）:
  • 検出器の分解能効果（スミアリング）を補正するため、D'Agostini 反復アルゴリズム（2 回反復）を用いて真の断面積分布を再構成。

3. 主要な貢献

• 統計量の飛躍的向上: 反ニュートリノ誘起の CC $\pi^0$生成に関するこれまでにない高統計データ（11.38 $\times 10^{20}$ POT 照射）を提供。
• 多変数微分断面積の測定: $\pi^0$の運動量・角度、ミューオンの運動量・角度、四元運動量転送 $Q^2$、ハドロン系の不変質量 $W_{EXP}$ に対する微分断面積を初めて詳細に報告。
• モデル依存性の評価: 複数のニュートリノ相互作用生成器（GENIE, NuWro, NEUT）との包括的な比較を行い、特に共鳴領域におけるモデルの精度を定量的に評価。
• FSI の影響の明確化: 原子核内での最終状態相互作用（FSI）が$\pi^0$の運動量分布に与える影響を実証的に示し、モデルへのフィッティングの重要性を強調。

4. 結果

• 断面積の形状とモデルとの比較:
  • GENIE モデル: 全体的にデータと最もよく一致するが、$\Delta(1232)$共鳴が支配的な領域（低$\pi^0$運動量、低$W_{EXP}$）では、実験データを約 20% 過小評価する傾向が見られた。
  • 他のモデル（NuWro, NEUT）: 全運動量範囲でデータを過小評価する傾向が強く、特に$\pi^0$運動量ピーク付近で 10-30% の不一致が見られた。
  • $\pi^0$運動量分布: 0.2 GeV/c 付近にピークがあり、GENIE モデル（FSI あり）は形状を良く記述するが、低運動量領域で過小評価。FSI を含まないモデルでは一致が悪化（$\chi^2/NDF$が 0.8 から 2.9 に増加）。
  • ミューオン運動量: 1 GeV/c 以上でデータがモデルより高い値を示し、より硬いミューオン分布を好む傾向が見られた。
• $\Delta(1232)$共鳴領域: $W_{EXP} < 1.4$ GeV の領域では、データが GENIE 予測を系統的に上回っており、$\Delta^0(1232) \to n + \pi^0$過程の寄与がモデルより大きいことを示唆。
• 不確実性: 統計誤差は約 5.5% だが、系統誤差（フラックス 8.9%、モデル 6.5%、検出器 7.4%）が支配的であり、総合的な不確実性は約 16.2%。

5. 意義

• ニュートリノ振動実験への寄与: 電子ニュートリノ出現事象の背景となる NC $\pi^0$生成の理解を深め、将来のニュートリノ振動実験（DUNE や Hyper-K など）における系統誤差低減に不可欠な入力データを提供する。
• 相互作用モデルの改善: 特に反ニュートリノ領域における$\Delta$共鳴生成と非共鳴過程のモデル化、および原子核内でのパイオン輸送（FSI）の記述精度向上に対する重要な制約条件となる。
• 将来の測定への基盤: 本測定で得られた共分散行列を含むデータは公開され、理論家や他の実験グループによるモデル調整の基礎資料として利用可能である。

結論として、本研究は反ニュートリノ誘起中性パイオン生成に関する最も精密な測定であり、現在のニュートリノ相互作用モデル（特に共鳴領域）が実験データを完全に記述できていないことを明らかにし、モデルの改良を促す重要な成果となりました。