Analysis and implications of the spatio-spectral morphology of the Fermi Bubbles

技術要約：フェルミ・バブルルの空間・スペクトル形態の解析と示唆

問題提起
フェルミ・バブルル（FB）（銀河面の上方および下方へ約50°にわたって広がる巨大な$\gamma$線構造）の起源と動力源は依然として不明である。その形態は銀河中心（GC）を起源としていることを示唆しているが、$\gamma$線放射を駆動する微視的なメカニズムについては、ハドロン・チャネル（周囲のガスとの陽子の非弾性衝突）とレプトン・チャネル（宇宙線電子による星間放射場光子の逆コンプトン散乱）の二つの主要な経路の間で議論されている。従来の解析は、粗い空間表現（例：緯度平均のスライス）に依存しており、位置依存的なスペクトルの変化を見落としている可能性がある。さらに、粒子輸送と冷却に要するタイムスケールは、バブルルが活動銀河核（AGN）活動によって駆動されているのか、あるいは星形成バーストによって駆動されているのかについて、重要な制約を与える。本論文は、観測された$\gamma$線スペクトルを再現するために必要な親宇宙線（CR）集団を制約するため、高解像度のピクセル単位のスペクトル解析を行うことで、これらの不確実性を解決することを目的とする。

手法
著者らは、Platzら（2023）によって提供された10年間の全天Fermi/LATデータのテンプレートフリーな再構成を利用している。具体的には、「非ダスト」（$I_{nd}$）成分であるM2モデルを用いる。これにより、熱的ダストマップと相関のない放射を分離し、銀河前景に対してFBを明確に抽出する。空はHEALPixスキーム（$N_{side}=128$）を用いてピクセル化されており、FBをカバーする約20,000個のピクセルが得られる。

放射をモデリングするために、著者らは以下の手順を踏む：

1. 3D幾何学の構築： 目視でトレースされた境界を銀河の$z$軸の周りに回転させることで、FBを回転体として近似する。これにより、視線方向の深さと、各光線に沿った物理量の積分が可能となる。
2. 星間放射場（ISRF）のモデリング： Popescuら（2017）による軸対称放射伝達モデルを用い、各視線方向における位置依存のISRF（恒星、ダスト、およびCMB成分）を決定する。モデルのグリッド外（$z > 10$ kpc）の領域については、べき乗則または$z^{-1}$のスケーリングを用いて外挿する。
3. フォワード・モデリング： Naima Pythonパッケージを使用し、各ピクセルの観測された$\gamma$線スペクトルを、3つの親集団モデル（べき乗則（PL）、指数関数的カットオフを持つべき乗則（EPL）、および折れ曲がりを持つべき乗則（BPL））に対してフィッティングする。
4. ハドロン対レプトン： ハドロンモデルでは、基準となる一様な水素数密度 $n_H = 10^{-3} \text{ cm}^{-3}$ を仮定し、結果をそれに応じて正規化する。レプトンモデルでは、各ピクセルについて導出された平均ISRFに対する逆コンプトン（IC）散乱を計算する。
5. 適合度： 簡約$\chi^2$を用いてモデルを評価する。データ分解における相関のある不確実性を考慮し、この指標は絶対的な統計的有意性ではなく、相対的なモデル性能を比較するためのものである。

主な結果

• スペクトル形状： 陽子および電子の両方に対する単純なべき乗則（PL）モデルは、データへの適合が悪く、棄却される。EPLおよびBPLモデルは、ハドロンおよびレプトンの両方のシナリオにおいて、等しく良好な適合を示す。
• スペクトル指数と硬化： 低エネルギーにおける最良適合スペクトル指数は、両方の宇宙線陽子（CRp）および宇宙線電子（CRe）において、バブルルの南端に向かって顕著な硬化を示す。この傾向は先行研究の結果と一致しているが、以前に解析されたよりも高い緯度まで及んでいる。
• 折れ曲がり／カットオフエネルギー：
  • EPLモデルでは、カットオフエネルギーは数TeVのオーダーであり、緯度による変化はほとんど見られない。
  • BPLモデルでは、折れ曲がりエネルギーは約100 GeVであり、緯度による変化は最小限である。
  • BPLモデルにおける折れ曲がり前後のスペクトル指数の変化は $\lesssim 0.5$ である。これは、シンクロトロン放射やIC損失による冷却ブレイクで期待される $\Delta \alpha \approx 1$ よりも大幅に小さく、スペクトルの折れ曲がりが冷却による特徴ではないことを示唆している。
• 放射場の寄与： レプトンモデルでは、全Fermi/LATエネルギー帯域において、ISRFのダスト成分が$\gamma$線生成を支配しており、次いで恒星成分（高エネルギーではクライン＝ニシナン効果により抑制される）、そしてCMBが続く。したがって、IC損失はバブルル全体を通してトムソン領域で発生していると近似できる。
• エネルギー密度と圧力：
  • ハドロン： 総宇宙線エネルギーは $\sim 1.7 \times 10^{57} / n_{H,-3}$ erg と推定される。宇宙線圧力は熱ガス圧よりも支配的である可能性が高いが、仮定された標的ガス密度に強く依存しており、数倍程度の差にとどまる。
  • レプトン： CReエネルギー密度は銀河面からの距離とともに増加し、バブルルの縁付近でピークに達する。これは、リム・ブライトニング（縁輝光）の現象と定性的に一致している。
• 冷却タイムスケール：
  • ハドロン： 冷却タイムスケールは極めて長く（$n_H = 10^{-3} \text{ cm}^{-3}$ の場合 $\sim 100$ Gyr）、これはバブルルが非常に古いか、定常状態から遠いことを意味する。
  • レプトン： 高エネルギー電子の冷却時間は極めて短く（バブルルのキャップ部で $\lesssim 1$ Myr、銀河中心付近では $\sim 0.1$ Myrまで低下する）。

意義と示唆
本論文は、スペクトル形状のみに基づけば、ハドロンモデルとレプトンモデルのどちらも決定的に排除することはできないと結論付けている。なぜなら、EPLおよびBPLモデルはどちらの粒子種に対しても等しくデータをよく説明するからである。しかし、その結果は物理的シナリオに対して厳格な制約を課している：

1. レプトンの輸送 vs. その場加速： 短い冷却時間と、バブルルの縁に向かってCReエネルギー密度が増加するという観察結果の組み合わせは、電子が銀河中心で注入され、外側へと輸送されるというモデルに対して深刻な課題を突きつけている。電子を、その冷却時間である $\sim 1$ Myr以内にバブルルのキャップ（$\sim 8$ kpc）まで輸送するには、X線現象論と矛盾する膨張速度（$\sim 8000$ km/s）が必要となる。著者らは、**その場での再加速（in-situ re-acceleration）**が、観察された形態およびエネルギー密度分布に対するより自然な説明であると示唆している。
2. ハドロンの年齢制約： 長いハドロン冷却タイムスケールは、もしバブルルが若い（例：$< 100$ Myr）場合、要求されるエネルギー予算は星形成よりもAGN起源を必要とすることを意味する。逆に、もしバブルルが古い（数Gyr）のであれば、宇宙線集団が蓄積するのに十分な時間があったという条件下で、星形成起源の可能性が残る。
3. スペクトルの硬化： 南端に向かうスペクトルの硬化と、冷却ブレイク（$\Delta \alpha < 0.5$）の欠如は、スペクトルの特徴が放射冷却損失ではなく、注入または加速メカニズムによって駆動されていることを示唆している。

要約すると、本研究はフェルミ・バブルルに対する初のピクセル単位のスペクトル制約を提供しており、単純なべき乗則は不十分であるが、指数関数的カットオフまたは折れ曲がりを持つモデルがデータをよく説明することを実証している。得られた知見は、核からの相対論的輸送、あるいは特にレプトン・チャネルにおいては、バブルルの縁における分散的なその場加速のシナリオを強く支持している。