Pulsar based modeling of point spread function of Fermi Large Area Telescope

技術要約：パルサーに基づくFermi Large Area Telescopeの点広がり関数（PSF）のモデリング

問題提起
銀河外天体周囲の広がったガンマ線放射の探索における感度は、本質的に望遠鏡の点広がり関数（PSF）の精度によって制限される。Fermi Large Area Telescope（LAT）のインストゥルメント・レスポンス関数（IRF）で提供されているPSFモデルの不正確さは、光源の広がりに関する偽陽性の検出を招く可能性がある。過去の研究では、モンテカルロ・シミュレーションに基づく初期のPSFモデリングが、観測されたイベント分布を再現できなかったことが示されており（例：Ando & Kusenko 2010）、これがブラザー（blazar）周囲の広がり放射の誤った主張につながった。Fermi/LATコラボレーションは、改良されたモンテカルロ・シミュレーションを用いてPSFモデルを更新しているが、検証研究は現在、主に68%包含半径に依存している。著者らは、広がった放射の探索という文脈において、この指標は不十分であると主張している。なぜなら、異なるPSFプロファイル（例：ガウス型と、べき乗則の裾を持つ指数関数型）は、同じ包含半径を共有しながらも、その形状は大きく異なる可能性があるからである。

手法
著者らは、解析的なモデルをパルサー（Vela、Geminga、およびCrab）からの信号と比較することにより、Fermi/LATのPSFモデルを評価し改善するためのデータ駆動型アプローチを提案している。

• データ選択: 残留電荷宇宙線の背景放射が低いP8R3_SOURCEVETO_V3イベントクラスを使用する。データは0.1–300 GeVのエネルギー範囲における17年間（2008年3月～2026年1月）にわたる。
• 信号の分離: 広がったパルサー風星雲や拡散銀河放射から点光源信号を分離するために、著者らはパルス位相解析を用いている。パルサーのエフェメリス（暦）を使用して、「オンパルス」フラックスから「オフパルス」フラックスを減算する。これにより、非周期的な背景成分を除去しつつ、点光源の形態を保持する。
• モデル比較: 著者らは、解析的なPSFモデル（「コア」成分と「テール」成分を記述する2つのキング関数の和）を、様々なエネルギービン、オフアクシス角、およびイベントタイプ（PSF0–PSF3、FRONT/BACK）にわたるスタックされたパルサーデータに適合させる。
• パラメータ最適化: 標準的なIRFパラメータと、PSFパラメータ（$f_c, \sigma_c, \gamma_c, \sigma_t, \gamma_t$）を自由変数とした再適合セットとの適合度（$\chi^2$）を比較する。

主要な貢献と結果

1. PSFの不一致の特定: 解析の結果、Fermi/LATのIRFファイルで提供されているPSらパラメータが、観測されたパルサーデータと一致しないことが明らかになった。エネルギーおよびオフアクシス角の半数以上のビンにおいて、IRFパラメータはパルサーデータから導出された90%信頼区間の外側に位置している。これらの不一致は、PSFのコアとPSFのテールの両方に現れ、残差は広い角度範囲にわたって10%以上に達する。
2. 改良されたPSFモデルの開発: キング関数のパラメータをパルサーデータに直接適合させることで、著者らは観測と一致する改良されたPSFパラメータのセットを導出した。元のモデルと改訂モデルの68%包含半径は数パーセント以内で一致しているが、形状については、特にテールの部分で改訂モデルは大きく異なっている。
3. Mrk 501の再解析: 著者らは、改良されたPSFを、以前に（Webar et al. 2025）強さ $B \sim 10^{-15}$ G の固有間磁場（IGMF）を示す広がった放射があるとして報告されたブラザー、Mrk 501の再解析に適用した。
   • 標準的なIRF PSFを用いた場合、解析は $\sim 4.57\sigma$ ($TS = 20.92$) の有意度で広がった放射の検出を再現する。
   • 改訂された、パルサーに基づくPSFを用いた場合、広がった放射の証拠は消失する。テスト統計量は $TS = 4.08$($\epsilon = 0.40 \pm 0.12$) まで低下し、統計的に有意ではない。
   • 改訂されたPSFは、点光源信号自体の適合度も向上させた（点光源の $TS$ を49,375から49,579へ増加させた）。
   • その結果、Webar et al. (2025) が主張したIGMFパラメータに対応する広がった放射モデルは、改訂されたPSFを用いると95%信頼レベルで否定される。

意義と限界
本論文は、真の広がった放射と、PSFモデルの不一致によるアーティファクトを区別するために、精密なPSFモデリングがいかに重要であるかを実証している。著者らは、以前に報告されたMrk 501周囲の広がった放射の検出は、実際のデータを正確に記述していないPSFパラメータを使用したことによる偽陽性である可能性が高いと結論付けている。

著者らは、自身のデータ駆動型アプローチにおける特定の限界についても述べている。それは、高統計のパルサーのパルス信号に依存している点であり、これらは $\sim 10$ GeV以下でのみ利用可能である。このエネルギーを超えると、パルサーのフラックスは抑制されるため、高エネルギー探索においてはモンテカルロに基づくPSF決定に戻る必要がある。しかし、0.1–10 GeVの範囲については、本研究は現在のIRFパラメータ化が、広がった放射の探索における系統誤差を避けるために改訂を要することを確立している。