Precision three-Dimensional Atmospheric Neutrino Flux Calculation Based on Honda Flux Model

問題提起

大気ニュートリノ（$\nu_{\text{atm}}$）は、ニュートリノ振動の研究、拡散超新星ニュートリノ背景放射（DSNB）の探索、暗黒物質（ダークマター）検出、および核子崩壊探索において、極めて重要な背景事象およびプローブである。HKKMS、Bartol、FLUKAといった3次元（3D）フラックス計算フレームワークは存在するものの、現在および将来の実験の文脈において、以下の3つの具体的な制限に直面している：

1. 古い入力モデル： 従来の計算は、最新の実験測定を反映していない古い一次宇宙線モデル（例：AMS01ベース）や地磁気モデル（例：IGRF2010）に依存していることが多い。
2. 低エネルギー物理学の欠如： これまでの3Dフラックス計算では、地球内部における大気ミューオン（$\mu_{\text{atm}}$）の伝搬、その後の停止、および崩壊または原子核捕獲を考慮したものは存在しない。この欠落は、DSNBやダークマター探索に不可欠なサブ100 MeVのエネルギー領域において重要である。
3. サイト比較の不完全性： 多様な地磁気環境（南極から低緯度サイトまで）におけるフラックスの系統的な差異が、統一された最新のフレームワークを用いて包括的に比較されていない。

手法

著者らは、確立されたHKKMS15フレームワークに基づき、いくつかの主要な手法的更新を取り入れた、包括的な3D大気ニュートリノフラックス計算を提示している。

• 地球内部におけるミューオンの伝搬： 本シミュレーションでは、初めて、地球の地殻や海洋を透過する荷電二次粒子（主にミューオン）を追跡している。本フレームワークは、それらのエネルギー損失、停止、およびその後の運命をモデル化している：
  • $\mu^+$ は自由に崩壊する（$\mu^+ \to e^+ \nu_e \bar{\nu}_\mu$）。
  • $\mu^-$ は軌道内で崩壊するか、あるいは原子核捕獲を行う（$\mu^- + p \to \nu_\mu + n$）。
  • 停止媒体の特定には、CRUST1.0 グローバル地殻モデルを使用している。
  • このプロセスにより、低エネルギーニュートリノの補足的なフラックスが生成される。特に原子核捕獲は、$\sim95$ MeVまで及ぶ $\nu_\mu$ を生成する。
• 更新された入力モデル：
  • 一次宇宙線： AMS01ベースのモデルを、AMS02、BESS-polar、PAMELA、JACEE、RUNJOB、およびCREAMのデータを取り入れたAMS02ベースのモデルに置き換えた。
  • 地磁気場： サイト依存の剛性カットオフとミューオンの曲がりを正確に計算するため、IGRF2020モデルに更新した。
  • 強粒子相互作用： $\mu_{\text{atm}}$ フラックスの不一致を1–100 GeVで5%以内に抑えるよう調整された、JAMおよび修正されたDPMJET-IIIモデルの組み合わせを採用した。
  • 大気： サイトおよび時間依存の空気密度プロファイルを得るために、NRLMSISE-00 グローバル大気モデルを使用した。
• シミュレーションの範囲： 計算は、10 MeVから$10^4$ GeVのエネルギー範囲にわたって、7つの検出器サイト（JUNO、Super-Kamiokande、CJPL、KM3NeT/ORCA、IceCube、DUNE、およびTRIDENT）を対象としている。横方向の運動量と地球の曲率を保持するために、完全な球面幾何学と3D地磁気を利用している。

主な貢献

1. 地球内伝搬ミューオンの初導入： 本論文は、これまでのすべての3D計算（HKKMS11/15、Bartol、FLUKA）には存在しなかった、質的に新しいソースとして、地球内部で停止した大気ミューオンによる寄与を導入している。
2. グローバル・マルチサイト・フレームワーク： 南極（IceCube）から低緯度サイト（TRIDENT）に至るまで、地磁気効果の直接比較を可能にする、7つの主要な実験拠点に対する精密なサイト固有のフラックス予測を提供する。
3. 低エネルギーの精度： 詳細な3Dモデリングが欠けていた領域に対し、これらのサイトにおける初の精密な100 MeV以下のフラックス結果を提供する。
4. 再推定された不確かさ： ミューオン制約法（Hondaら）を用いた強粒子不確かさの推定を更新し、従来の計算と比較して系統誤差を大幅に減少させた。

結果

1. $E_\nu > 100$ MeV におけるフラックス特性：

• サイト依存性： 地磁気カットオフのため、$E_\nu < 10$ GeVにおいて顕著なサイト依存の差異が現れる。例えば、1 GeV以下では、IceCubeにおける $\nu_\mu$ フラックスはJUNOの約2倍である。10 GeV以上では、フラックスは各サイト間で収束する。
• 天頂角および方位角： フラックスは、生成幾何学と地磁気遮蔽によって駆動される強い天頂角依存性（上下および水平・垂直比）を示す。方位角の非対称性（東・西効果）はサブGeVのエネルギーで顕著であるが、高エネルギーでは消失する。
• HKKMS15との比較： 更新された一次宇宙線のスペクトル形状（40 GeV以下で高フラックス、それ以上で低フラックス）および再調整された強粒子相互作用モデルに起因して、HKKMS15と比較して2%–10%の偏差が観察される。

2. $E_\nu < 100$ MeV におけるフラックス特性（ミューオン伝搬の寄与）：

• グローバルな絶対増加量： 地球内部におけるミューオン伝搬による寄与は、グローバルな効果である。フラックスの絶対的な増加量（$\Delta \phi_\nu$）はサイトに依存せず、10–50 MeVの範囲でピークに達する。
• 相対的な増強： ベースラインのフラックスは（地磁気カットオフにより）サイトごとに異なるため、相対的な増強率は反比例する。ベースラインのフラックスが低いサイト（例：JUNO）は、高フラックスのサイト（例：IceCube）よりも大きな分率的増加を示す。
• フレーバーの階層性： 寄与は $\nu_e \approx \bar{\nu}_\mu > \nu_\mu > \bar{\nu}_e$ の階層に従う。これは、停止ミューオンの電荷非対称性（$\mu^+ > \mu^-$）および、Michelエンドポイント（$\sim52.8$ MeV）を大きく上回る$\sim95$ MeVまで及ぶ $\nu_\mu$ を生成する原子核捕獲プロセスによって駆動されている。

3. 不確かさの推定：

• ミューオン制約法を用いることで、強粒子不確かさ（$\delta_{\text{had}}$）は大幅に減少した。1–10 GeVの水平方向において、$\delta_{\text{had}}$ は $\nu_\mu/\bar{\nu}_\mu$ に対しては約3%、$\nu_e/\bar{\nu}_e$ に対しては約4%に減少した。
• 全相対不確かさ（$\delta_{\text{tot}}$）は、1–10 GeVの範囲（水平方向）で約5%–7%に減少しており、これは従来の推定と比較して著しい改善である。

意義

本論文は、これらの結果が幅広い物理プログラムに対して不可欠かつ精密な入力となることを主張している：

• ニュートリノ振動： 3Dフラックス予測と減少した不確かさは、大気ニュートリノ実験（JUNO、Super-Kamiokande、Hyper-Kamiokande、DUNE、KM3NeT/ORCA、IceCube）のための改良された入力として機能する。
• DSNB探索： 精密な低エネルギーフラックス計算（新しいミューオン伝搬成分を含む）は、10–30 MeVの窓におけるDSNB探索の支配的な不可避の背景事象を定義する。
• ダークマター探索：
  • 間接探索： 更新されたフラックスは、WIMP消滅（IceCube）および軽いダークマターの消滅（JUNO、Borexino、Super-K）の探索における背景特性の向上をもたらす。
  • 直接探索： 精密な低エネルギー予測は、直接ダークマター検出器（LZ、XENONnT、PandaX-4T）の感度を制限する「ニュートリノの霧（neutrino fog）」を制約する。
• 核子崩壊： より良く制約された大気ニュートリノ背景は、現在および次世代の検出器における核子崩壊探索の感度限界に直接影響を与える。

著者らは、現在の計算は堅牢なベースラインを提供するものであるが、局所的な山のプロファイル、時間依存の太陽変調、およびNA61/SHINEやLHCfのデータを用いた更新された強粒子モデルなどの将来的な洗練が進むことで、特に低・中間エネルギー領域における系統的な不確かさがさらに減少すると結論付けている。