Probing Lorentz Invariance Violation in Cosmogenic Neutrino Propagation with KM3-230213A

技術要約：KM3-230213Aを用いた宇宙線起源ニュートリノ伝播におけるローレンツ不変性破れの探索

問題提起
KM3NeTコラボレーションによって検出された超高エネルギー（UHE）ニュートリノイベントKM3-230213A（推定エネルギー $220^{+570}_{-110}$ PeV）は、観測されたスペクトルと組成を記述することに成功しているUHE宇宙線（UHECR）モデルから導かれる標準的な宇宙線起源ニュートリノフラックスの予測に対し、テンション（乖離）を提示している。KM3-230213Aの起源については議論が分かれているものの、このイベントは非標準的なニュートリノ物理学のユニークなプローブとして機能する。具体的には、本研究では、ニュートリノセクターにおける超光速的なローレンツ不変性破れ（LIV）を調査している。超光速シナリオ（$v > c$）では、ニュートリノは運動学的に不安定となり、標準模型では禁止されている真空ペア放出やニュートリノ分裂（$\nu \to 3\nu$）といった崩壊チャネルが開かれる。これらの崩壊は、地球における観測されるニュートリノフラックスを大幅に変化させる可能性があり、これによりKM3-230213Aのイベントと、Pierre Auger ObservatoryおよびIceCubeによって確立された宇宙線起源ニュートリノのフラックス上限値との間のテンションを解消できる可能性がある。

手法
本研究は、UHECRの伝播とそれに続く宇宙線起源ニュートリノの生成をモデル化するために、CRPropa 3.2ツールキットに基づいたシミュレーションフレームワークを採用している。

• ソースモデル: 5つの代表的な一次核（H, He, N, Si, Fe）を、Pierre Auger Observatoryの最良適合組成パラメータに一致するように重み付けしてシミュレートしている。2つの主要なシナリオがテストされている：「ベストフィット」組成（中間質量核が支配的）と、「プロトン強化」シナリオ（より軟らかいスペクトルと高い剛性カットオフを持つ10%の純粋な陽子成分を加えたもの）。ソースの進化は $\Psi(z) \propto (1+z)^m$ とパラメータ化され、$m$ は 0.0 から 2.0 の範囲をとる。
• LIVの実装: 超光速ニュートリノ（$\delta_{\nu, n} > 0$）に対する修正分散関係（MDR）を実装している。光子や宇宙線に必要な複雑なカスケード計算とは異なり、ニュートリノの伝播は解析的に扱われる。生存確率 $P_{surv}$ は、支配的なニュートリノ分裂チャネル（$\nu \to 3\nu$）のエネルギー依存崩壊率を宇宙論的な距離にわたって積分することによって計算される。
• フラックス再構成: 標準的なフラックス（$\Phi_{std}$）に対して、LIVを考慮するためのリウェイト・スキームを適用している。修正されたフラックス（$\Phi_{LIV}$）は、生存した一次フラックスと、崩壊による二次フラックスの総和として構築される。ここで、娘ニュートリノは親のエネルギーの約 $1/3$ を保持する。この手法は、高エネルギーにおけるスペクトルの抑制と、低エネルギーにおける「蓄積（pile-up）」による増強をモデル化している。
• ベンチマーキング: KM3-230213Aイベントは、感度を評価するためのベンチマークとして使用される。解析では、LIVによって修正された予測フラックスを、Pierre Auger Observatory、IceCube（EHE限界）、およびグラショフ共鳴イベントからの観測制約と比較する。

主な結果

• スペクトル修正: 超光速LIVは特徴的なスペクトル歪みを誘起する。すなわち、最高エネルギー領域（ニュートリノが不安定になる領域）でのフラックスの抑制と、それに伴う中間エネルギー（PeVスケール）での増強（蓄積）である。
• パラメータ感度: 本研究は、線形LIV係数 $\delta_{\nu, 1}$ に対する特定の観測感度領域を特定している。係数が $10^{-24} \text{ eV}^{-1} \lesssim \delta_{\nu, 1} \lesssim 10^{-22} \text{ eV}^{-1}$ の範囲にある場合、KM3-230213Aの再構成エネルギー範囲（72 PeV – 1 EeV）内での期待イベント率が最大化され、かつ高エネルギー側でのIceCubeやPierre Augerからの上限値を破ることなく成立する。
• 組成依存性: 結果は、仮定されたUHECRソースの組成に強く依存する。最も保守的な「プロトン強化」シナリオは、より高いベースライン・フラックスをもたらす。保守的な「ベストフィット」（Augerのみ）シナリオでは、6 PeVにおける標準フラックスは無視できるほど小さく、KM3-230213Aのレベルまでフラックスを増強するのに十分な大きさのLIV係数は、おそらく他の制約と矛盾することになる。
• 検出確率: プロトン強化シナリオでは、超光速LIVの導入により、KM3NeT/ARCAにおける少なくとも1つの宇宙線起源ニュートリノイベントの検出確率が系統的に増加する。最適な検出確率は $\delta_{\nu, 1} \sim 10^{-24} - 10^{-22} \text{ eV}^{-1}$ のバンドで発生する。
• 統計的有意性: 著者らは、様々なLIV仮定の下でのKM3-230213Aイベントの統計的有意性（$\sigma$）を算出している。LIVは有意性を高めることができるが、利用可能な統計量（単一のイベント）は、形式的な制約を与えるには不十分である。有意性は、特定の係数値で低下した後、プラトーに達するという非単調な挙動を示す。

意義と主張
本論文は、ニュートリノセクターにおける超光速LIVに対する「観測感度領域」を特定したと主張しており、将来のニュートリノ望遠鏡へのテスト可能な予測を提供している。

• クリーンなプローブ: 著者らは、宇宙線起源ニュートリノはさらなる相互作用なしに伝播するため、複雑なカスケードシミュレーションを必要とせず、分散関係の直接的な解析的修正によってLIV効果を組み込むことができるため、「クリーン」な基礎物理学のプローブであると強調している。
• エネルギー・スケール: 制約されたLIVエネルギー・スケール（$E_{LIV} \sim 10^{22} - 10^{24}$ eV）は、プランクスケールよりも約5〜6桁低い。著者らは、これはLIV効果が追加の抑制因子を伴って現れる可能性のある、特定の量子重力理論（例：ループ量子重力、弦理論のコンパクト化）と一致していると述べている。
• 控えめな結論: 著者らは、KM3-230213Aイベントはベンチマークとして用いられているものの、LIVパラメータ空間を形式的に制約するための統計的有意性が欠如していること（変動はわずか $\sim 0.5\sigma$ である）を明示的に述べている。主な貢献は、その手法の感度を実証し、このアプローチの可能性を強調することにある。結果は、将来の施設（例：IceCube-Gen2, GRAND, POEMMA）とともに統計量が増加するにつれて、宇宙線起源ニュートリノがローレンツ不変性のより厳格なテストを提供するようになることを示唆している。
• 他セクターとの矛盾なし: 本研究は、これらの制約がニュートリノセクターに特有のものであり、UHECRスペクトルやガンマ線透過率から導かれるハドロンおよび光子セクターにおけるより厳しいLIVの境界値と矛盾しないことを指摘している。