Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

技術要約：超巨大ブラックホールコロナにおける衝突無衝突衝撃波による陽子加速

問題の所在
IceCubeニュートリノ観測所による最近の観測では、近傍の電波静穏な活動銀河核（AGN）、特にセイファート銀河NGC 1068から、高エネルギーニュートリノの顕著な過剰が検出されている。これらの観測結果は、一つの不一致を提示している。すなわち、高エネルギーニュートリノは検出されている一方で、対応する高エネルギーガンマ線のフラックスは、標準的なハドロン相互作用モデルによる予測よりも著しく低い。これは、ニュートリノがガンマ線に対して不透明な領域、おそらく中心の超巨大ブラックホール（SMBH）を取り囲む高温のコロナ内で生成されており、そこではガンマ線がX線光子との$\gamma\gamma$対生成によって減衰していることを示唆している。

観測されたニュートリノスペクトルを説明するためには、親陽子がコロナ内で$\sim 100$ TeVのエネルギーまで加速される必要がある。しかし、SMBHコロナにおける主要な加速メカニズムは依然として未確立である。拡散衝撃波加速（DSA）は有力な候補であるが、その生存可能性は不透明である。具体的には、コロナの衝撃波は低い音速マッハ数（$M_s \sim 2-10$）および準相対論的な速度を持つことが予想されており、これらは伝統的に粒子加速には非効率とされる領域である。さらに、二温度プラズマ（$T_i \gg T_e$）の挙動や、このような環境における陽子対電子の加速効率の具体的な制約についても、第一原理に基づく制約が欠如している。

手法
これらの不確実性に対処するため、著者らはSMILEIコードを用いた、第一原理に基づく一次元粒子インセル（PIC）シミュレーションを用いている。本研究は、SMBHコロナの物理的パラメータに適合させた準平行衝突無衝突衝撃波（$\theta = 30^\circ$）に焦点を当てている。

シミュレーション・サーベイは、セイファート銀河の観測的制約に基づいた広範なパラメータ空間を系統的に探索している：

• 音速マッハ数 ($M_s$): 遷速（trans-sonic, $M_s \approx 2$）からより高い値（$M_s \approx 30$）まで。
• アルヴェンマッハ数 ($M_A$): $4$から$46$ まで。
• 衝撃波速度 ($v_{sh}$): 準相対論的速度（$v_{sh}/c \approx 0.17 - 0.66$）。
• イオン対電子温度比 ($T_i/T_e$): 衝突無衝突コロナで予想されるデカップリングした熱状態を反映し、$0.01$から$50$ まで。

シミュレーションでは、計算の実現可能性を確保しつつ、主要な無次元パラメータを保持するために、低減されたイオン対電子質量比（$m_i/m_e = 100$）を使用している。実行は、電磁不安定性が飽和し、非熱的な冪乗則の裾（power-law tails）が発達するのに十分な時間（$t \gtrsim 400 \Omega_{ci}^{-1}$）にわたって行われた。

主要な結果

1. 堅牢な陽子加速: 本研究は、DSAがSMBHコロナにおける陽子加速の堅牢かつ効率的なメカニズムであることを示している。低い音速マッハ数にもかかわらず、衝撃波は一貫して衝撃波の運動エネルギーの約10%を非熱的イオンへとチャネル化している。この効率は、全パラメータ・サーベイを通じて安定しており、$M_s \approx 2$ のケースを含めても安定していることは、低$M_s$の衝撃波が非効率な加速器であるという概念に異を唱えるものである。
2. 加速ダイナミクス: 加速される陽子の最大エネルギーは、衝撃波構造が準定常状態に達した後、時間に対して線形に増加する。加速率は、DSA理論が予測するスケーリング、$t_{acc} \propto \gamma \beta_{sh}^{-2} \Omega_{ci}^{-1}$ に従う。著者らは、正規化された加速係数 $C^{-1} \equiv d(\gamma \beta_{sh}^{-2})/d(t\Omega_{ci})$ を導出した。これは普遍的な値 $C^{-1} \approx 0.1$ を示し、テストされたすべてのパラメータにおいて2倍未満の変動しか示さない。
3. 電子加速の非効率性: 陽子とは対照的に、電子加速は非常に非効率であり、非熱的電子エネルギー分率は1%未満（典型的なコロナ的 $T_i/T_e$ 比では $\sim 10^{-3}$）にとどまる。電子加速の効率は、$T_i/T_e$ 比が増加するにつれて低下する。特定の強磁場増幅を伴う準相対論的なケースでは、電子の冪乗則の裾が高エネルギーへと伸び続けることができず、増幅された磁場が衝撃波を電子に対して事実上超光速（superluminal）にする状況が断続的に発生することで、持続的なDSAが阻害されることが示唆されている。
4. 不安定性メカニズム: 散乱に必要な上流の磁気乱流は、ほとんどのパラメータ空間において、主に非共鳴不安定性（NRI、またはBellの不安定性）によって駆動されている。共鳴不安定性（RI）は、非常に低いアルヴェンマッハ数（$M_A \lesssim 8$）においてのみ支配的となる。

NGC 1068への示唆とニュートリノ生成
これらの結果をNGC 1068に適用することで、著者らは以下を実証している：

• タイムスケール: 陽子は、関連する冷却（ベーテ・ハイラー、フォトメソン、$pp$）および動力学的（自由落下）タイムスケールよりも短いタイムスケールで、$\sim 100$ TeVまで加速可能である。
• エネルギー予算: 観測されたニュートリノ光度を説明するために必要な非熱的陽子パワーは、全衝撃波パワーの $\sim 3\%$ と推定される。これは、シミュレーションで測定された $\sim 10\%$ の加速効率の範囲内に十分に収まっている。
• ガンマ線欠損: 電子加速の低い効率と、ガンマ線に対するコロナの高い不透明度は、観測された高エネルギーガンマ線の欠損を自然に説明し、同時にニュートリノが脱出することを可能にしている。

意義と主張
本論文は、SMBHコロナという特定の文脈におけるDSAに対する、初の第一原理に基づいた自己整合的な制約を提供したと主張している。主な意義は、低音速マッハ数や高イオン・電子温度比といった困難な条件下においても、衝突無衝突衝撃波が高エネルギー陽子を生成するための実行可能かつ効率的なメカニズムであることを確立した点にある。

著者らは、本研究の結果がニュートリノ生成のハドロンモデルを強く支持するものである一方で、局所的なキネティック・シミュレーションの結果をAGNコロナのような複雑な天体に適用するには、慎重な検討が必要であることを明示している。エネルギー予算の分析は、本論文で確立された微視的物理学に基づくオーダー・オブ・マグニチュード（桁違い）の評価として提示されており、微視的なプラズマ物理学と巨視的な高エネルギー現象の間の架け橋となっている。本研究は、最大陽子エネルギーに関するすべての不確実性や、3Dにおける不安定性の詳細な遷移を解決すると主張しているわけではなく、それらにはさらなる調査が必要であると述べている。