A Multiwavelength Interpretation of HESS J1857+026 Emission Using the Fermi-LAT, VERITAS, and HAWC Observatories

技術要約：HESS J1857+026放射のマルチ波長解釈

問題提起
HESS J1857+026は、高エネルギーなパルサーPSR J1856+0245と一致する、広がった未同定の超高エネルギー（VHE）$\gamma$線源である。このソースはパルサー風星雲（PWN）の候補であるが、そのマルチ波長的な性質は複雑なままである。先行研究では、GeV放射とTeV放射の広がりにおける空間的な不一致や、確認されたX線対応天体の欠如が明らかになっており、支配的な放射メカニズムに疑問を投げかけている。具体的には、その放射が純粋にレプトン的（PWNからの逆コンプトン散乱）なのか、純粋にハドロン的（分子ガスと相互作用する超新星残骸からのパイ中間子崩壊）なのか、あるいはその両方の組み合わせなのかが不明である。さらに、進化したPWNの直近環境における粒子拡散特性は、平均的な星間物質（ISM）と比較して、依然として十分に制約されていない。

手法
本研究は、15年間のFermi-LATデータ（300 MeV–2 TeV）、VERITASデータ（0.3–10 TeV）、および新しいHAWCデータ（1–316 TeV）を用いた、HESS J1857+026の包括的なマルチインストルメント解析を提示する。解析は以下の段階を経て進行する：

1. データ解析:

   • Fermi-LAT: Pass 8 SOURCEクラスのデータに対して、ビン化尤度解析を行った。関心領域（ROI）は、4FGL-DR3カタログとディフューズ・テンプレートを用いてモデリングされた。空間的な広がりのテストでは、点源、放射状ディスク、および放射状ガウス型テンプレートを比較した。
   • VERITAS: Gammapyフレームワークを用い、30時間の露光データに対して3Dビン化最大尤度解析を実施した。ソースは、対称なガウス形状を持つ冪関数（パワーロー）としてモデル化した。
   • HAWC: Pass 5データセット（2860日）に対して、フォワードフォールディング最大尤度解析を適用した。複数の空間テンプレート（ガウス型 vs 拡散ベース）およびスペクトル形状（冪関数、ログパラボラ、指数関数的カットオフを持つ冪関数）をテストし、最良の適合モデルを選択するために、ベイズ情報基準（BIC）および赤池情報量基準（AIC）を用いた。

2. 放射モデリング:

   • 定常的モデル: 純粋にハドロン的な（SNR）シナリオおよびレプト・ハドロン的な（PWN/SNR）シナリオを検証するため、NAIMAパッケージを用いて広帯域スペクトルエネルギー分布（SED）のフィッティングを行った。
   • 非定常的モデル: 進化したPWNの膨張と冷却をシミュレートするために、進化モデル（Gelfand et al. 2009と同様）を採用した。マルコフ連鎖モンテカルロ（MCMC）法を用いて、系の年齢、磁場強度、および粒子注入スペクトルのような物理パラメータを制約した。

3. 拡散解析:

   • 全てのインストルメントから放射輝度プロファイル（表面輝度プロファイル）を抽出した。これらのプロファイルは、粒子拡散係数 $D(E)$ とISMに対する拡散の抑制率を制約するために、解析的な拡散ベースのモデル（式16）を用いてフィッティングされた。

主要な結果

• スペクトルおよび空間的特性:

  • Fermi-LATデータは、半径 $0.38^\circ \pm 0.04^\circ$ の放射状ガウス型で最も良くフィットする広がったソースを示しており、スペクトル指数は $\Gamma \approx 2.07$ である。スペクトル破断が観測されており、5 GeV以下では軟らかい指数（$\Gamma \approx 2.63$）、5 GeV以上では硬い指数（$\Gamma \approx 1.99$）となっている。
  • VERITASは有意度 $15.2\sigma$ でこのソースを検出し、スペクトル指数 $\Gamma = 2.48 \pm 0.06$ およびガウス幅 $0.236^\circ \pm 0.017^\circ$ を測定した。
  • HAWCデータは、ソースが $\sim 37$ TeVまで広がっていることを示している。最良適合の空間モデルは、磁場 $1 \mu$Gにおいて、単純なガウス型よりも拡散テンプレートを支持しており、カットオフエネルギーは $E_c \approx 14$ TeV である。
  • ソースのサイズはエネルギーの増加とともに減少しており（GeVエネルギーでの $0.38^\circ$ からHAWCエネルギーでの $0.26^\circ$ まで）、進化したPWNと一致するエネルギー依存の形態を示している。

• 放射の起源:

  • レプト・ハドロン vs ハドロン: 低エネルギー（$<10$ GeV）成分については純粋なハドロンモデルを完全に排除することはできないが、レプト・ハドロンモデルの方が広帯域データに対して優れたフィットを示す。検出されたSNRシェルが存在しないこと、および進化モデルから導出された周囲のガス密度が低いこと（$n_0 \lesssim 0.13$ cm$^{-3}$）は、TeV放射における支配的なハドロン起源を否定するものである。
  • PWNの進化: 進化モデルは、この系が進化したPWNであり、年齢 $\tau = [16, 21]$ kyr、低磁場強度 $B = [0.4, 1.6] \mu$Gであることを示唆している。これらのパラメータは、フェーディングするX線放射（シンクロトロン冷却による）とエネルギー依存の形態を説明する。この低い磁場も、TeV放射における局所的な光子場への逆コンプトン散乱の優位性を支持している。

• 粒子拡散:

  • HESS J1857+026周辺の拡散係数は、平均的な銀河系ISMと比較して抑制されていることが分かった。50 TeVにおいて、拡散係数は $D_{PWN} \sim 2.14 \times 10^{28}$ cm$^2$ s$^{-1}$ （$B=1 \mu$Gの場合）であり、これはISMの値（$D_{ISM} \sim 1.2 \times 10^{30}$ cm$^2$ s$^{-1}$）よりも約2桁低い。
  • この抑制は、剛性依存性（スペクトル指数 $\delta$）ではなく、主に正規化（$D_0$）において生じており、$\delta$ はコルモゴロフ（$\delta=1/3$）またはクライコネン（$\delta=1/2$）の乱流モデルと一致している。これは、ソース付近における局所的な磁気乱流の強化というシナリオを支持している。

意義および主張
本論文は、結合されたマルチ波長データが、進化したPWNの近傍における粒子拡散特性を制約するための強力なツールとなることを主張している。主要な結論は、HESS J1857+026はPSR J1856+0245によって駆動されるPWNである可能性が極めて高く、観測された$\gamma$線放射は相対論的電子による逆コンプトン散釈によって支配されているということである。

本研究は、低エネルギーのスペクトル成分（$E < 10$ GeV）には、依然としてホストまたは近傍のSNRからのハドロン的寄与、あるいは、同一のPWN内における異なる電子集団が含まれている可能性があることを強調している。しかし、高エネルギー放射は明確にレプトン的である。導出された拡散抑制因子は、他のTeV PWN（例：Geminga, Monogem）と一致しており、これらのソースの近くでは、一般的なISMと比較して粒子の輸送が著しく妨げられているという仮説を補強している。著者らは、低有意度でニュートリノの過剰が報告されているものの、$\gamma$線領域と比較してニュートリノ放射領域の大きな上限値は、そのような潜在的なニュートリノ信号に対して複雑な、あるいは無関係な起源を示唆していると述べている。