An Accretion Flare Interpretation for the Ultra-High-Energy Neutrino Event KM3-230213A

問題提起
KM3NeTによって検出された、エネルギー約220 PeVの超高エネルギー（UHE）ニュートリノイベント「KM3-230213A」は、その天体物理学的起源を特定する上で大きな課題を提示している。99%の包含領域内に17個のブレーザーが特定されているものの、より大きな有効面積と長い運用履歴を持つIceCubeによる同様のイベントの非検出との間の緊張関係が、その解釈を複雑にしている。具体的には、定常的な拡散宇宙線起源は高い信頼度（$\sim > 3\sigma$）で否定されているが、過渡的な点源シナリオは、この緊張を緩和できる可能性がある。主要な候補であるMRC 0614-083は、ニュートリノ検出と一致するX線フレアや、変動する可視光/赤外線ライトカーブを含む多波長放射を示す。しかし、この天体の赤方偏移は未知であり、降着活動、X線フレア、およびUHEニュートリノ生成を結びつける物理的メカニズムはまだ確立されていない。

手法
著者らは、MRC 0614-083に対して、可視光ライトカーブを超巨大ブラックホール（SMBH）の降着率の時間進化のプロキシとして扱う、自己整合的な「降着フレア」モデルを提案している。その手法は以下の通りである：

1. 多波長データ解析： アーカイブおよび新規の観測（電波、可視光、赤外線、X線、$\gamma$線）を組み合わせ、スペクトルエネルギー分布（SED）およびライトカーブを構築する。可視光フレアは、増強された降着活動として解釈される。
2. ダスト・エコー・モデリング： 遅延して発生する赤外線（IR）放射を「ダスト・エコー」としてモデル化する。すなわち、降着フレアからの可視光放射が周囲のダストトーラス（$R_{IR} \sim 10^{17}–10^{18}$ cm）を加熱し、それが赤外線として再放射されるというものである。時間広がり関数を用いて、視線方向の初期ダスト成分と球状トーラス成分の両方を考慮しながら、赤外線ライトカーブをフィッティングする。
3. 粒子加速と相互作用： プロトンを等方的な放射領域（ダストトーラス）に、降着パワーに対する効率 $\epsilon_p = 0.2$ で注入する。モデルは超エディントン降着フェーズを仮定している。これらのプロトンは、光子メソン（$p\gamma$）過程を通じて熱的な可視光および赤外線光子と相互作用し、パイオンを生成し、それがニュートリノおよび二次的な電磁（EM）粒子へと崩壊する。
4. カスケード・シミュレーション： 著者らは、プロトン、電子/陽電子、光子、およびニュートリノの結合輸送方程式を解くために、オープンソースソフトウェアAM3を利用している。この完全な時間依存モデルは、$p\gamma$ 相互作用、 $pp$相互作用、$\gamma\gamma$ 消滅、ベテ・ハイティラー対生成、およびそれに続くシンクロトロン/逆コンプトン・カスケードを考慮に入れている。
5. 赤方偏移への感度： 赤方偏移（$z$）が未知であるため、モデルは $z = 0.5$ から $2.0$の範囲でテストされており、詳細な結果は基準となる$z=0.5$ のケースについて示されている。

主な貢献

• 統一的な解釈： 本論文は、単一の超エディントン降着フレアが、可視光の変動、遅延した赤外線ダスト・エコー、およびX線フレア、そしてUHEニュートリノイベントを同時に説明するという、統一的な枠組みを提供している。
• UHEニュートリノのメカニズム： ダストトーラスを $p\gamma$ 相互作用のターゲットとして特定し、降着フレア中に加速されたプロトンがUHE（$\sim 10^9–10^{10}$ GeV）に達し、観測されたエネルギーレベルのニュートリノを生成できることを実証した。
• 緊張の解消： 降着フレアの過渡的な性質を持つモデルは、ニュートリノフラックスが定常的な拡散背景放射ではなく特定の時間窓に集中しているため、IceCubeの非検出限界との間の緊張を緩和するのに役立つ。
• X線の起源： モデルは、観測されたX線フレアを $p\gamma$ 過程によって開始された電磁カスケードに帰属させており、$p\gamma$ 相互作用のタイムスケールと赤外線エコーの遅延の関数として、可視光のピークとX線ピークの間の時間差（$\sim 700–1000$ 日）を自然に説明している。

結果

• フレア特性： 基準となる赤方偏移 $z=0.5$ において、モデルはピーク降着率 $\dot{m}_{pk} \approx 1.3 M_9^{-1} \eta_{rad, -1}$ （超エディントン）を必要とし、静穏期には $\sim 0.74$ まで減少する。
• ニュートリノフラックス： 予測されるニュートリノスペクトルは $\sim 100$ PeV でピークを持ち、フラックスは $4 \times 10^{-12}$ erg s$^{-1}$ cm$^{-2}$ である。これは、低エネルギーにおけるIceCubeの測定値の結合 $E^{-2}$ フィットと一致しており、KM3-230213Aの $E^{-2}$ フィットの $3\sigma$ 不確定範囲内に収まっている。
• 電磁カスケード： モデルは観測されたX線ライトカーブとスペクトル（0.2–2.3 keV）を再現することに成功し、GeV領域におけるFermi-LATの上限値を満たしている。カスケードは $\gamma\gamma$ 消滅による電子・陽電子対によって支配されており、これが10 GeV以上の高エネルギー$\gamma$線を抑制しており、これは非検出の結果と一致している。
• プロトン加速： 解析により、磁場が $\sim 0.1$ Gであれば、ウィンド・ショック加速または再加速によって、プロトンが放射領域内でUHEまで加速可能であることが確認された。

意義
本論文は、MRC 0614-083が、特定の降着駆動型の過渡的事象をUHEニュートリノ検出に結びつけることで、超高エネルギー宇宙線の起源に関する最初の直接的な証拠となる有望な候補であることを主張している。自己整合的な記述は、ダストトーラスを粒子加速と相互作用の重要なサイトとして導入することで、標準的なブレーザー・ジェットのシナリオを補完するものである。著者らは、未知の赤方偏移が依然として主要な不確実性であることを強調しており、したがって、精密な光度スケーリングとモデルの妥当性を確認するために不可欠な赤方偏移を決定するための、可視光/紫外線分光測定を強く推奨している。本研究は、超エディントン降着フレアが、ブレーザーにおける、重要ではあるが過渡的なUHEニュートリノ集団の源となり得ることを示唆している。