Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

技術要約：ミラージュ・ソースと大規模なTeVハロー・パルサー・オフセット

問題提起
LHAASO（Large High Altitude Air Shower Observatory）およびHAWCによる最近の観測により、中年期パルサーの周囲に広がったTeVガンマ線放射（パルサー・ハロー）が特定されている。これらの天体を解釈する上での大きな課題は、パルサーの位置とTeV放射の重心との間に、頻繁に大きな空間的オフセット（ずれ）が観測されることである。電子・陽電子（$e^\pm$）が星間物質（ISM）中で等方的に拡散するという従来のモデルでは、このオフセットを自然に説明することができない。なぜなら、等方拡散は物理的スケールが複数の磁気コヒーレンス長（磁場の相関長）を超える場合にのみ有効な近似となるからである。しかし、超高エネルギー電子（$>100$ TeV）は冷却速度が極めて速いため、その伝播スケールは乱流的なISM磁場のコヒーレンス長よりも小さくなる可能性がある。その結果、粒子輸送は一様に拡散するのではなく、局所的な磁力線に沿って高度に非対称になる可能性がある。本論文では、このような非対称な伝播が、加速器から離れた位置に現れる「ミラージュ・ソース（蜃気楼的な光源）」や、顕著な空間的オフセットを生じさせ得るかどうかを調査する。

手法
著者らは、GPU加速された第一原理に基づいたテスト粒子シミュレーションを用い、100 TeVの電子が乱流磁場中に注入された際の軌跡を追跡する。

• 磁場生成： 乱流磁場は、フーリエ空間法を用いて生成された、4つの独立したコルモゴロフ・スペクトル（$n_{in} = -11/3$）を持つ磁場の重ね合わせとしてモデル化されている。このマルチスケール・アプローチ（グリッド間隔 $4 \times 10^{-9}$ から $5 \times 10^{-4}$ kpc）により、現実的な慣性領域と発散のない場（$\nabla \cdot B = 0$）を確保している。全磁場は、規則的な成分（$B_r$）と乱流成分（$B_t$）の両方を含む。
• 粒子伝播： 電子は、ローレンツ力の下での相対論的な運動方程式を解くためにボリス・プッシャー（Boris pusher）アルゴリズムを用いて伝播される。シミュレーションでは、シンクロトロン放射による放射エネルギー損失、および宇宙マイクロ波背景放射（CMB）と星間放射場（ISRF）に対する逆コンプトン（IC）散乱を考慮している。
• 統計的サンプリング： 初期エネルギー・スペクトル（$dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$）を効率的にサンプリングするため、重み付きマクロ粒子法を使用している。8,192個のマクロ粒子を注入し、各粒子は物理的な電子数に比例した重みを持つ対数エネルギービンを表す。
• 観測合成： シミュレートされた電子分布は、IC散乱行列を介して合成ガンマ線マップ（25–40 TeV帯）に変換される。これらのマップは、ガウス型点広がり関数（PSF）$\sigma_{PSF} = 0.43^\circ$（LHAASO-KM2Aを模したもの）で畳み込みが行われ、ポアソン統計的な背景ノイズモデルが加えられる。
• 光源特定： マルコフ連鎖モンテカルロ（MCMC）を用いたベイズモデル比較フレームワークが適用される。ガウス型ソースの数を増やしながらモデルをフィッティングする。コンポーネントを追加した際にテスト統計量（TS）が25以上増加した場合（$\Delta \log L > 12.5$）、そのソースは「特定された」とみなされる。

主要な貢献および結果
本研究は、パルサーの距離（$d$）、電子注入スペクトル指数（$\alpha$）、全磁場強度（$B$）、乱流コヒーレンス長（$L_c$）、および規則的磁場と乱流磁場の強度比（$B_r/B_t$）を変化させながら、パラメータ空間を系統的に探索している。

1. ミラージュ・ソースの形成： シミュレーションにより、非対称な伝播が単一の加速器から複数の明確なガンマ線源（「ミラージュ・ソース」）を作り出し得ることが確認された。これは、電子が磁力線に沿って伝播する際、投影効果によって、注入点から離れた場所に起源があるように見える場合に発生する。
2. ソースの多重性（$N$）：
   • 特定されたソースの数は、磁場強度（$B$）に最も強く影響される。高い $B$ は冷却を速め、電子をより小さなスケールに閉じ込めるため、複数の明確なソースを形成する可能性を低下させる（低い $N$）。
   • ソース数は、距離（$d$）およびコヒーレンス長（$L_c$）と負の相関があり、$B^2$ と正の相関がある。
3. ソース間の分離とフラックス比：
   • 分離： 特定されたソース間の最大分離距離は、$L_c$ と正の相関があり、$B^2$ と負の相関がある。
   • フラックス比（$\chi$）： メイン・ソースのフラックスとミラージュ・ソースのフラックスの比は、$B^2$ の増加とともに減少する。強い磁場においては、ミラージュ・ソースがメイン・ソースよりも明るくなることもある。
   • 規則的磁場 vs 乱流磁場： 分離距離は $B_r/B_t$ に対して非単調な依存性を示す。最大分離は $B_r/B_t \approx 1$ のときに発生する。このとき、規則的な磁場が方向性を与えつつ、乱流がオフセットを生むための十分な散乱を許容する。純粋な乱流、あるいは支配的な規則的磁場の場合、分離は小さくなる。
4. 空間的オフセット：
   • パルサーと最も近い特定されたソースとのオフセットは、$B_r/B_t \approx 1$ かつ $L_c$ が大きいときに最大となる。
   • 大きなコヒーレンス長を持つ領域（例：腕間領域や銀河ハロー）では、オフセットは最大50 pcに達する場合がある。

意義と主張
本論文は、TeVハローが関連するパルサーと空間的に切り離されているように見えるLHAASOの天体（1LHAASO J0216+4232u, J0212+4254u, J0206+4302uなど）に対し、定量的な説明を提供すると主張している。具体的には、これらの天体が、個別の加速器ではなく、単一のミリ秒パルサー（PSR J0218+4232）によって駆動されるミラージュ・ハローである可能性を示唆している。

本研究は、高エネルギーのレプトンにとって冷却時間が磁気コヒーレンス時間よりも短い場合、ISMにおける加速器付近の粒子輸送をモデル化する上で、等方拡散近似が不十分であることを強調している。投影効果が、乱流磁場においてオフセットや複数のソース形態を自然に生み出すことを示すことで、本論文は、観測されたTeVの空を標準的なパルサー風星雲モデルと整合させるメカニズムを提示している。著者らは、想定している合成乱流が（MHDにおける異方的なアルヴェン型乱流とは異なり）等方的であるものの、非対称な放射を引き起こすISM乱流の大規模なフィラメント構造を定性的に捉えていると述べている。