Benchmarking neutrino-nucleus quasielastic scattering model predictions against a missing energy profile obtained using a monoenergetic neutrino beam

NEUT 事象生成器内の 3 つの原子核基底状態殻モデルを、単色ニュートリノ源を用いた JSNS$^2$実験の欠損エネルギー分布測定結果と比較検証した結果、核内カスケードと核励起チャネルを考慮した場合、相対論的平均場モデルよりもスペクトル関数モデルが基底状態および分布の尾部の記述において優れていることが示された。

要約

🎯 全体のストーリー：「幽霊の標的当てゲーム」

ニュートリノは、物質をすり抜ける「幽霊のような粒子」です。しかし、稀に原子核（ここでは炭素の原子核）とぶつかることがあります。このぶつかり方を正しく理解しないと、ニュートリノの正体（質量や振動）を調べる実験が失敗してしまいます。

これまでの実験では、ニュートリノのエネルギーがバラバラ（広帯域）だったため、原子核がどう反応したかを詳しく見るのが難しかったです。
しかし、今回の研究では、**「エネルギーが完全に一定（モノクロ）なニュートリノ」を使って、「炭素の原子核」**にぶつける実験を行いました。

これは、**「暗闇でバラバラの石を投げる」のではなく、「同じ重さの石を正確に狙って投げる」**ようなもので、原子核の内部構造を非常に詳しく観察できるチャンスでした。

---

🧩 3 つの「シミュレーションモデル」の対決

研究者たちは、コンピュータ上で原子核の動きを予測する「3 つの異なるシミュレーション（モデル）」を用意しました。これらは、原子核の中がどうなっているかという「地図」のようなものです。

1. SF モデル（スペクトル関数）: 過去の電子散乱実験のデータを元に作られた、経験則に基づいた地図。
2. SF モデル（改良版）*: 上記の地図をさらに細かく、精密に描き直したもの（特に、原子核の「基底状態」という一番安定した状態を詳しく描こうとした）。
3. ED-RMF モデル: 相対性理論を厳密に使って計算した、理論的に堅固な地図。

これら 3 つのモデルが、実際の実験データとどれくらい合致するかを競わせたのです。

---

🔍 実験の仕組み：「消えたエネルギー」を探す

ニュートリノが炭素の原子核にぶつかると、陽子や中性子が飛び出します。
ここで重要なのが**「欠損エネルギー（Missing Energy）」**という概念です。

• イメージ:
  • あなたが 100 円の硬貨（ニュートリノのエネルギー）を投げて、箱（原子核）にぶつけました。
  • 箱から飛び出してきたもの（飛び出した粒子）のエネルギーを測ると、合計 80 円しかなかったです。
  • 「残りの 20 円はどこへ行った？」
  • この「行方不明の 20 円」こそが、原子核の内部構造（中性子がどれくらい強くくっついているか）を表す鍵なのです。

JSNS2 という実験施設は、この「行方不明のエネルギー」の分布を精密に測定しました。

---

🏆 結果：どのモデルが勝った？

実験結果とシミュレーションを比較したところ、以下のようなことがわかりました。

1. 「核の脱励起（NucDeEx）」と「内部カスケード（FSI）」が重要

原子核にぶつかった後、飛び出した粒子が他の粒子とぶつかり合ったり（内部カスケード）、原子核が落ち着こうとしてエネルギーを放出したり（脱励起）するプロセスをシミュレーションに含めないと、実験データと全く合いませんでした。

• 比喩: ボウリングでピンを倒した直後、倒れたピンが他のピンに当たり、さらに転がったりする動きまで計算しないと、実際のスコアと合いません。

2. 勝者は「SF モデル」

• SF モデル: 実験データと最もよく合いました。特に、原子核の「一番安定した状態（基底状態）」の予測が正確でした。
• SF モデル & ED-RMF モデル*: これらは理論的には優れているはずでしたが、実験データとはズレがありました。
  • SF*: 精密すぎる地図だったはずですが、特定のエネルギー領域で「強すぎる」予測をしてしまいました。
  • ED-RMF: 理論的に堅いですが、原子核の「壁」の描き方が少し硬すぎて、実際の「柔らかい」動きを捉えきれていませんでした。

3. 「閾値（しきい値）」の罠

面白い発見として、**「単一の粒子を叩き出すには、ある一定のエネルギーが必要だ（しきい値がある）」というルールをシミュレーションに適用すると、「どのモデルも実験データと矛盾しなくなる」**という結果が出ました。

• 意味: 現在のモデルは、物理的にありえない（エネルギーが足りないのに粒子が出てくる）ような予測を少し含んでいました。それを修正すれば、どのモデルも「合格」ラインに入ってしまうのです。

---

💡 この研究の重要性は？

この研究は、**「ニュートリノ実験の基礎となる『原子核の地図』が、実はまだ完璧ではない」**ことを示しました。

• これまでの常識: 「広帯域のニュートリノを使っているから、細かい原子核の構造は気にしなくていい」と思われていました。
• 今回の発見: 「エネルギーを一定にして測ると、モデルの細かなズレがバレてしまう。特に、原子核がどう反応して落ち着くか（脱励起）や、粒子同士のぶつかり合い（FSI）の扱いが重要だ」ということがわかりました。

🚀 まとめ

この論文は、**「ニュートリノという幽霊を捕まえるための『網』の作り方を、新しい『モノクロの光』を使って見直した」**という話です。

• 勝者: 経験則に基づいた「SF モデル」が、今のところ最も実測に近い。
• 課題: 理論的に完璧なはずのモデルが、実際の「原子核の揺らぎ」や「粒子のぶつかり合い」をうまく表現できていない。
• 未来: この結果を元に、ニュートリノ実験（T2K や将来の DUNE など）の精度をさらに高め、宇宙の謎に迫るための「より正確な地図」を作ろうという呼びかけです。

つまり、**「原子核という箱の奥底を、もっと正確に描くための第一歩」**を踏み出した研究なのです。

技術概要

論文の技術的概要：単色中性子ビームを用いた欠損エネルギープロファイルによるニュートリノ - 原子核準弾性散乱モデルのベンチマーク

1. 背景と課題 (Problem)

ニュートリノ振動実験（T2K、NOvA、将来の DUNE や Hyper-K など）における系統誤差の主要な要因は、ニュートリノ - 原子核相互作用のモデル化、特に原子核基底状態の記述と核効果（核内効果）の理解にあります。

従来のニュートリノ事象生成器（NEUT など）は、主に包括的（inclusive）な理論モデルに依存しており、最終状態のレプトンの運動量のみを関数とする微分断面積を予測していました。しかし、近年はより排他的（exclusive）なモデル（電子散乱データに基づくスペクトル関数モデルなど）の実装が進んでいます。

主な課題:

• 加速器中性子ビームの広帯域性: 従来の加速器実験では、入射ニュートリノエネルギーの分布が広いため、欠損エネルギー（Missing Energy, $E_m$）の再構成が困難であり、原子核基底状態のモデリングと核効果（最終状態相互作用など）を分離して検証することが極めて困難でした。
• モデルの検証不足: 排他的な電子散乱データとの比較は行われていますが、ニュートリノ散乱データ、特に単一エネルギーのニュートリノを用いた精密な検証は不足していました。

2. 手法とアプローチ (Methodology)

本研究は、静止カイオン崩壊（KDAR）源から得られる単色ニュートリノ（235.5 MeV）を用いた JSNS2 実験の最新の欠損エネルギー測定結果 [13] を利用し、NEUT 事象生成器（バージョン 6.0.3）に実装された 3 つの排他的原子核基底状態モデルをベンチマークしました。

検証対象モデル:

1. Benhar Spectral Function (SF): 電子散乱データに基づく従来のスペクトル関数モデル。
2. Benhar SF (SF):** 高分解能の新しいスペクトル関数。$^{11}\text{B}^*$ の基底状態と第一・第二励起状態（2.13 MeV, 5.02 MeV）を区別し、より精密なガンマ線放出を記述可能。
3. ED-RMF (Energy-Dependent Relativistic Mean Field): 因子分解を行わない相対論的平均場モデル。ハドロンテンソルをテーブル化して実装。

シミュレーション構成:

• 相互作用: 主に荷電流準弾性散乱（CCQE: $\nu_\mu n \to \mu^- p$）をシミュレート。2p2h 過程や共鳴過程も含まれるが、寄与は小さい。
• 核内効果:
  • 核内カスケード (Intranuclear Cascade): NEUT の標準的なカスケードモデル（Bertini モデルに基づく）を用いて、非弾性の最終状態相互作用（FSI）をシミュレート。
  • 核脱励起 (NucDeEx): 専用の事象生成器 NucDeEx を用いて、核の脱励起過程（ガンマ線や軽粒子の放出）をシミュレート。
• 欠損エネルギーの定義:
  $$E_m = E_\nu - E_\mu - \sum T_p - E_\gamma$$
  ここで、$E_\nu$ は入射ニュートリノエネルギー、$E_\mu$ はミューオンエネルギー、$\sum T_p$ は陽子の運動エネルギーの和、$E_\gamma$ は核脱励起によるガンマ線エネルギー。
• 統計解析: JSNS2 の測定値と生成器の予測値を比較し、共分散行列を用いた $\chi^2$ 解析と p 値を計算。
  • 閾値処理: 単一核子放出閾値（1NKO: $^{12}\text{C}$ の中性子分離エネルギー 18.72 MeV）以下の物理的に不可能な領域を含むかどうかで 2 通りの解析を実施。

3. 主要な成果と結果 (Key Contributions & Results)

A. モデル間の比較

• FSI と脱励起なしの場合: 全てのモデル（SF, SF*, ED-RMF）が測定値と一致せず、特に基底状態のピーク（$1p_{3/2}$）とその高エネルギー側（テール）の形状が再現できていませんでした。
• FSI カスケードの導入: 非弾性 FSI を含めることで、全てのモデルにおいてピークからの強度がテール側へシフトし、$\chi^2$ 値が改善しました。
  • SF モデル: FSI 導入により、p 値が 0.78 まで向上し、測定値と統計的に有意な一致を示しました。
  • SF および ED-RMF モデル:* FSI 導入後も p 値が 0.05 未満（棄却）となりました。
    • ED-RMF: $1p_{3/2}$ ピークを過大評価し、その直後のビンで過小評価する傾向がありました。また、1NKO 閾値以下の非物理的な領域への寄与（ガウス関数の両側性による）が $\chi^2$ を増大させました。
    • SF:* 基底状態をデルタ関数で記述するため 1NKO 閾値以下の寄与はありませんが、$1p_{3/2}$ ピークの強度を過大評価し、$1p_{3/2}$ と $1s_{1/2}$ の遷移領域の強度を過小評価していました。

B. 単一核子放出閾値（1NKO）の影響

• 閾値カットなし: 1NKO 閾値以下の非物理的な領域を含めると、SF モデル以外は棄却されました。
• 閾値カットあり: 1NKO 閾値（18.72 MeV）以上のみを解析対象とした場合、FSI と脱励起を含めた全てのモデルが p 値の基準（0.05）を満たし、統計的に棄却されなくなりました。
  • ただし、$\chi^2_{N-1}$ 分析（各ビンの寄与を個別に評価）により、モデルごとの不一致の根源が異なることが明らかになりました。SF* は遷移領域、SF と ED-RMF はテール部分に不一致が残っています。

C. 追加的なパラメータ感度解析

• 核ポテンシャル: ED-RMF 以外のポテンシャル（RPWIA, EDAI, rEDAI）も検討されましたが、基本的な傾向は ED-RMF と同様でした。
• FSI 強度: NEUT カスケードの再散乱確率を±30% 変化させたところ、FSI 強度を増加させることが、特に SF モデルの適合度をさらに向上させることが示されました。

4. 結論と意義 (Significance)

1. 単色ニュートリノビームの重要性: KDAR 源を用いた単色ニュートリノ実験は、原子核のスペクトル関数（特に中性子）を直接かつ詳細に検証する強力な手段であることを実証しました。広帯域ビームでは見逃されていたモデルの微細な欠陥（閾値以下の非物理的寄与やピーク形状の不一致）を浮き彫りにしました。
2. モデルの優劣と改善点:
   • SF モデル: 現在のところ、FSI と核脱励起を適切に組み合わせた場合、JSNS2 のデータと最もよく一致するモデルです。
   • ED-RMF と SF:* 理論的に堅牢なモデルですが、現状の実装では特定のエネルギー領域（ピーク形状、遷移領域、テール）で測定値と乖離しています。特に ED-RMF における 1NKO 閾値以下の非物理的な寄与の除去や、ガウス関数による基底状態の記述方法の再検討が求められます。
3. FSI モデリングの課題: 測定データは、現在の NEUT 実装における FSI の強度をさらに増加させることを支持しています。また、異なる生成器（NuWro など）の FSI モデルとの比較や、より詳細な核効果（短距離相関 SRC など）の扱いが今後の課題です。
4. 将来のニュートリノ実験への貢献: 本研究で得られた知見は、T2K や将来の DUNE/Hyper-K などの高精度ニュートリノ振動実験において、原子核効果に起因する系統誤差を低減し、ニュートリノ質量階層や CP 位相の測定精度を向上させるために不可欠です。

総じて、本論文は単色ニュートリノビームを用いた精密測定が、ニュートリノ相互作用モデルの検証において決定的な役割を果たすことを示し、NEUT 生成器内の原子核モデルの改良に向けた具体的な指針を提供しました。