An 18-year Fermi-LAT stacking limit on GeV $γ$-ray emission from particle-accelerating colliding-wind binaries

技術要約：粒子加速を行う衝突風連星系からのGeV $\gamma$線放射に関する、Fermi-LATによる18年間のスタッキング限界

問題提起
大質量衝突風連星系（CWB）は、超音速の星風が衝突する風衝突領域（WCR）を有しており、そこでは強力な衝撃波が発生し、粒子を相対論的なエネルギーへと加速させることが可能である。非熱的電波シンクロトロン放射はこれらの系における粒子加速を裏付けており、$\eta$ CarもGeV $\gamma$線源として明確に検出されているが、CWB集団レベルでのGeV放射については未解明のままである。限られた露光量による過去の専用サーベイでは、個々の系からの有意な放射は見出されていない（$\gamma^2$ VelorumはVela領域の混雑により$\eta$ Carほど確実ではないが、例外として挙げられる）。本研究が取り組む中心的な問いは、粒子加速を行うCWB（PACWB）の集団全体が、個別の閾値下にあるGeV信号を集合的に生成しているのかどうかである。

手法
著者らは、1 GeV以上のFermi-LATデータ（2008年8月から2026年5月まで）の17.8年分を解析した。エネルギー閾値は、1 GeV以下での広範な応答とは異なり、点広がり関数（PSF）が著しく絞り込まれる（68%包含角で$\lesssim 0.8^\circ$）ことから選択された。これにより、連星を銀河面近傍の他の天体から分離することが可能となる。

1. サンプル選択: 本研究では、De Becker & Raucqのカタログに含まれる61個の確定した（「リストA」）PACWBを利用した。

   • 除外: 個別に検出されている2つの系（$\eta$ Carおよび$\gamma^2$ Vel）は、未検出の集団による集合的な信号を探索するために、対象から除外された。
   • クリーニング: 明るい天体（FL16YのSignif$_{Avg} \ge 5$）の周囲に厳格な1.0$^\circ$のマスクを適用した。これらの系は、PSFが他の明るい天体（星形成領域や銀河中心の源など）と混ざり合うため、個別の源としてのフラックス測定が信頼できないことから、マスク内で除外された。
   • 結合: 相関のあるスタッキングを避けるため、1.5$^\circ$未満の隔たりがある系は、最も高いテスト統計量（TS）を持つものを選んで結合した。
   • 最終サンプル: このプロセスにより、互いに孤立した高緯度領域の「クリーン」なサンプル $N=6$ が得られた。

2. スタッキング解析:

   • 光子指数（$\Gamma \in [1.0, 3.5]$）および積分フラックス（$\log_{10} F \in [-13, -6]$）に関する2次元尤度スキャンを実行した。
   • 結合テスト統計量（$TS_{stack}$）は、$N$ 個のソースの個別の尤度比の総和として計算された。
   • コントロールフィールド: 帰無分布を確立するために、200個のコントロールフィールドを生成した。これらは、銀河緯度（$|b|$）、局所的なソース密度（2$^\circ$以内のFL16YのSignif$_{Avg}$の総和）、および銀河拡散強度（$> 1$ GeV）においてターゲットサンプルと一致させた。極めて重要な点として、偶然の明るい天体との一致を帰無分布に反映させるため、コントロールフィールドに対しても、ターゲットと同じ1.0$^\circ$のカタログマスクを適用した。

主要な結果

1. 非検出: 6つの系のクリーンなスタックは、累積TS 3.91を示した。これは、環境が一致したコントロールフィールドの帰無分布（中央値 TS = 10.0）と一致しており、p値は0.83（$-0.9\sigma$）であった。PACWB集団からの集合的なGeV放射の証拠は見出されなかった。
2. 上限値: 光子指数 $\Gamma = 2$ と仮定した場合、平均的な1ソースあたりのフラックスの95%上限値は以下の通りである：
   $$F(> 1 \text{ GeV}) \lesssim 1.1 \times 10^{-11} \text{ ph cm}^{-2} \text{ s}^{-1}$$
   これは、光度限界 $L_\gamma \lesssim 1.3 \times 10^{31} (d/\text{kpc})^2 \text{ erg s}^{-1}$ に相当する。
3. 効率の制約: 代表的な風の運動エネルギー（$L_{wind} \sim 3.3 \times 10^{36} \text{ erg s}^{-1}$）に対して、$\gamma$線生成効率は以下のように制約される：
   $$\eta = L_\gamma / L_{wind} \lesssim 4 \times 10^{-6} (d/\text{kpc})^2$$
   サンプルの平均距離（$\sim 1.7$ kpc）において、これは$\eta$ Carの効率（$\eta \sim 10^{-3}$）よりも2桁ほど低い。
4. アーティファクトの特定: カタログマスクを除去して全集団（59系）をスタックしたところ、$\gtrsim 3.5\sigma$ の偽の過剰（excess）が現れた。著者らは、これが銀河面におけるPACWBと明るいカタログ天体との間の偶然の一致によるアーティファクトであることを示しており、混雑した領域におけるスタッキングにおける重要な方法論的注意点を浮き彫りにしている。

粒子加速への制約
単一ゾーンモデルを用いて、著者らはこのフラックス限界を粒子加速の制約へと変換した：

• 電子: この限界は、電波データから推測される磁場（数mGから$\sim 1$ G）において、電子加速効率（$\eta_e$）を$\sim 10^{-3}$に制限する。標準的な効率（$\eta_e \sim 0.01–0.1$）を達成するには、磁場が磁気・光子エネルギー等分配値（$B \gtrsim 3–7$ G）と同等かそれ以上である必要があるが、これはCWBで通常推測される値よりも高い。
• 陽子: 陽子は $pp \to \pi^0$ 崩壊による放射を行う前にWCRから脱出する可能性が高いため、この限界は陽子の加速をより弱く制約する（$\eta_p \lesssim 0.11$）。

意義と主張
本論文は、$\eta$ Carが、GeV $\gamma$線を放出する共通の集団の中で最も明るいメンバーであるというよりは、極端な風の出力と離心率を持つ特異な天体であると結論付けている。この非検出は、大多数のPACWBが、風のエネルギーをGeV $\gamma$線へと変換する効率が著しく低いか、あるいはLATのバンドにおける放射を抑制するスペクトル特性（例：カットオフ）を持っていることを示唆している。

本研究は、銀河面でのスタッキングにおいては、緯度だけでなく、局所的なソース密度や拡散強度についても一致させたコントロールフィールドを使用する必要性を強調している。ターゲットとコントロールに対して同一のマスキングを行わない限り、スタッキング解析は、カタログの偶然の一致によって偽陽性を生み出す可能性があることを確立した。得られた制約は、高緯度の、汚染のないサブセットに対しては強固なものであるが、最も有望な候補（密な星形成領域にあるもの）は、必要なマスキングのために依然として未調査のままである。