The GPSP Atlas: High-resolution Galactic Photon Survival Probability Atlas for Very-High- and Ultra-High-Energy Gamma-Ray Astronomy

技術要約：GPSP アトラス

問題提起
LHAASOやHAWCといった広視野観測装置によって推進される銀河系PeV天文学の台頭に伴い、光子–光子（$\gamma\gamma$）吸収の補正がガンマ線データ解析の極めて重要な要素となっている。ガンマ線は銀河系内物質（ISM）を伝播する際、主に銀河内放射場（ISRF）および宇宙マイクロ波背景放射（CMB）との相互作用による対生成（$\gamma\gamma \to e^+e^-$）を通じて減衰を受ける。この吸収は「銀河の霧」を生み出し、特に30–200 TeV（赤外線/可視光のISRFが支配的）およびPeVエネルギー付近（CMBが支配的）において、観測されるスペクトルを抑制する。

個別の吸収補正は特定の天体に対して適用されてきたが、急速に増加する超高エネルギー（VHE）および極高エネルギー（UHE）天体集団の解析を標準化するための、統一された高解像度のフレームワークが分野全体として欠如している。さらに、銀河系内のPeVatron（粒子をPeVエネルギーまで加速する源）を特定したり、ローレンツ不変性の破れ（LIV）のような基礎物理学を検証したりするためには、エネルギー、銀河座標、距離の全パラメータ空間にわたる生存確率 $P(E_\gamma, l, b, d)$ を正確に把握する必要がある。

手法
著者は、任意のソース構成に対して事前計算された生存確率を提供する、高解像度4Dルックアップテーブルである銀河光子生存確率（GPSP）アトラスを構築した。

• 物理フレームワーク: 減衰は、対生成のブライト・ウィーラー断面積から導出される光学的厚さ $\tau$ を用いて計算される。計算は直線伝播を仮定しており、二次的な電磁カスケード効果は無視し、一次ガンマ線成分の減衰に焦点を当てている。
• ISRFモデル: 系統的な不確実性を評価するため、アトラスは2つの独立した公開されているGALPROPベースのISRFモデルを用いて構築されている。
  • R12 (Robitaille et al. 2012): 顕著な渦状腕構造と内側のダスト「ホール」を特徴とする。
  • F98 (Freudenreich 1998): 非軸対称なバルジを持つ、より滑らかで中心に集中したダスト分布を予測する。
    両モデルは、等方的なCMBと組み合わされている。
• グリッド構成: アトラスは以下の4Dグリッドをカバーしている。
  • 銀河経度 ($l$): $0^\circ$ から $359.9^\circ$ （ステップ $0.1^\circ$）。
  • 銀河緯度 ($b$): $-5^\circ$ から $+5^\circ$ （ステップ $0.1^\circ$）。
  • 距離 ($d$): 0.1 kpc から 20.0 kpc （ステップ 0.1 kpc）。
  • エネルギー ($E_\gamma$): 1 TeV から 10 PeV （81個の対数ビン）。
    これにより、約 $5.89 \times 10^9$ 個のデータエントリーが生じる。
• LIV拡張: ローレンツ不変性の破れ（LIV）による減速的な効果を調査するために、コンパニオン・アトラスであるGPSP-LIVが構築された。このバージョンは、分散関係がLIVスケールパラメータ $\Lambda$ によって変化する、修正された対生成キネマティクスを組み込んでいる。これは、LIV強度（$\log_{10}(\Lambda/E_{Pl})$ から $-4.5$から$-2.0$）を含む5Dパラメータ空間をカバーしている。
• データ形式: アトラスはFITS形式でリリースされ、Gammapyのような標準的な解析ツールと互換性があり、32ビット浮動小数点数の生存確率を含んでいる。

主な貢献

1. 統一されたフレームワーク: 1 TeVから10 PeVまでの全パラメータ空間 $P(l, b, d, E_\gamma)$ におけるガンマ線の生存確率を記述する、初の銀河系フレームワーク。
2. 系統的不確実性の定量化: R12およびF98モデルに基づく並行アトラスを提供することで、ISRFに関連する系統的な不確実性を直接評価することを可能にした。これらは、赤外線相互作用が支配的な領域および銀河中心部において最も顕著であることが示された。
3. 斬新な可視化: 「ガンマ線生存フェーズスペース」の導入。これは、従来の1次元プロファイルを超えて、銀河の透明度と地平線を可視化する2次元のエネルギー–距離表現である。
4. LIVの統合: LIV修正された対生成断面積を取り入れた、大規模な銀河系光子生存アトラスの中で初の事例であり、銀河内伝播におけるLIV効果の系統的な研究を可能にする。

結果

• 銀河透明度マップ: マップは、銀河の不透明度が高度に構造化されていることを明らかにしている。$\sim 100–300$ TeVの領域では、高い赤外線光子密度と長い経路長により、渦状腕の接線方向（例：ノルマ腕やスカットム–クルックス腕）に沿って減衰が増強される。
• エネルギー依存性: エネルギーがPeV領域に向かって上昇するにつれ、減衰はますます等方的なCMBによって支配されるようになる。その結果、生存確率マップは、異方的で構造化されたパターンから、ほぼ等方的な分布へと遷移する。
• ガンマ線地平線: 地平線距離（$P_{surv} \approx 0.37$ と定義）は、$\sim 2–3$ PeVにおいて最小の $\sim 7$ kpcに達した後、より高いエネルギーで再び拡大する。銀河反中心方向が最も透明な方向である。
• モデル間の差異: R12とF98の間の分数偏差は、50–200 TeVの範囲で最大約9%（銀河面内では最大約13%）に達するが、CMBが支配するPeVエネルギーでは1%未満に低下する。
• LIV効果: 減速的なLIVシナリオは、対生成断面積を抑制し、運動学的閾値をシフトさせ、PeVエネルギーにおける銀河の透明度を大幅に向上させる。2–3 PeVの範囲で標準的なガンマ線地平線を超えて放出される天体は、LIVを探索するための有望なターゲットとして特定される。
• 応用:
  • 銀河中心 (GC): 銀河中心の拡散放射の解析によれば、H.E.S.S.のエネルギー域では減衰は中程度であるが、UHEエネルギーにおいては、固有のPeVatronシグネチャーを覆い隠す可能性のある有意なスペクトルの急峻化をもたらす。
  • Cygnus X-3: Cygnus X-3の固有スペクトルの再構成により、ISRFモデルの違い（R12 vs F98 vs P17）が「赤外線ショルダー」の特徴や最小生存確率のエネルギーに変動をもたらすことが示されたが、現在の測定の不確かさがこれらの系統的な差異よりも支配的である。

意義と主張
本論文は、GPSPおよびGPSP-LIVアトラスが、銀河系PeVatron、拡散ガンマ線放射、および基礎物理学の探索のための公開可能な基盤を確立すると主張している。著者は、ガンマ線天文学がPeV領域に突入するにつれ、減衰補正がUHE解析の標準的な構成要素になることを強調している。

本研究は、以下のためのツールとして提示されている：

• 増大するUHE天体集団の解析の標準化。
• 固有の加速限界と吸収によるスペクトル特徴の区別。
• 観測されたスペクトルをLIV修正された伝播モデルと比較することによる、LIVシナリオをテストするための系統的な枠組みの提供。

著者はLIV制約に関する記述については控えめであり、現在の解析は「原理の実証（proof-of-principle）」としてのデモンストレーションであると述べている。彼らは、堅牢なLIV制約には、個々の天体の検出だけでなく、よく特性化されたソース特性を持つ集団研究や、完全な尤度ベースの解析が必要になるであろうと述べている。同様に、現在のフレームワークは二次的なカスケード効果や異方的なISRFの取り扱いを省略しており、これらは将来のイテレーションで計画されていることも認めている。