Population synthesis of Galactic middle-aged pulsar wind nebulae II. Observational signatures of superefficiency

技術要約：銀河系中齢パルサー風星雲の集団合成 II. 超効率性の観測的シグネチャー

問題提起
パルサー風星雲（PWN）がそのホストである超新星残骸（SNR）と相互作用する際、SNRの逆ショックによって駆動される「リバーブレーション（残響）フェーズ」に入る。このフェーズにおいて、星雲は著しい圧縮を受け、これにより磁場の増幅と粒子の断熱加熱が引き起こされる。このプロセスの特定の帰結として、「超効率性（superefficiency）」、すなわち特定の周波数帯におけるPWNの積分放射輝度が、同時期のパルサーの回転エネルギー損失率（$\dot{E}_{SD}$）を上回る状態が生じる。先行研究（Torres & Lin 2018; Torres et al. 2019など）では、この現象の物理的可能性は確立されているが、銀河系全体の集団におけるその普遍性と、詳細な観測的シグネチャーについては、依然として多くが未解明のままであった。さらに、既存のモデルの多くは、簡略化された「単一ゾーンの純粋な薄いシェル（one-zone purely thin-shell）」の処方箋に依存しており、リバーブレーションの複雑な流体力学、粒子の再処理、および磁場進化を正確に捉えきれない可能性がある。

手法
著者らは、中齢のPWN（年齢 $\le 10^5$ 年）の合成的な銀河系集団を用いて、超効率性がどのように顕現するかを調査するために、集団合成アプローチを採用している。本研究では、計算効率の高い単一ゾーンTIDEコードと、リバーブレーションフェーズ中の周囲のSNRの流体力学的構造を自己整合的にモデル化するラグランジュ・モジュールを組み合わせたハイブリッドTIDE+Lフレームワークを用いている。これにより、シェル（殻）の厚層化、遅延した圧縮、および複数のショック反射を現実的に扱うことが可能となり、これらは粒子および磁場の進化において極めて重要である。

主な手法ステップは以下の通りである：

• 集団生成： 1回の実現ごとに、合計2,400個のシミュレーションから抽出された1,600個のソースからなる合成集団を生成する。これは、銀河系のコア崩壊型超新星率、親星の質量（サルピーターIMF）、およびパルサーの誕生特性と一致している。
• 比較モデリング： 詳細なリバーブレーション物理の影響を分離するため、TIDE+Lフレームワークの結果を、標準的なTIDEモデル（ラグランジュ・リバーブレーション・モジュールを持たない純粋な薄いシェル・アプローチ）と比較する。
• 進化の追跡： 集団を3つの進化段階（自由膨張、リバーブレーションによる圧縮、およびリバーブレーションによる圧縮後）を通じて追跡する。
• スペクトル解析： $L_{band}/\dot{E}_{SD}$ を算出することで、12の周波数帯（電波からTeVまで）にわたって超効率性を評価する。解析には、個々のソースの追跡、アンサンブルSED（スペクトルエネルギー分布）、および特定の周波数帯と圧縮因子（CF）によって定義された電子スペクトルの解析が含まれる。
• 統計的堅牢性： 集団の平均期待数と確率的な分散を定量化するために、1,000回のランダムな集団実現から結果を導出する。

主要な貢献と結果

1. 超効率性の普及度：

   • TIDE+Lフレームワークを用いると、現在の銀河系における相当数の超効率的ソースが予測される。遠赤外線（FIR）帯において最も高い普及度を示し、約315個のソース（集団の19.7%）が超効率性を示す。
   • 中赤外線（MIR）、近赤外線（NIR）、可視光、UV、およびGeV帯においても、かなりの数のソースが予測される。
   • 対照的に、純粋な薄いシェルを用いたTIDEモデルは、すべての帯域において有意に少ない超効率的ソースを予測し（例：FIRにおいてTIDE+Lの315に対し約116）、X線帯域では実質的にゼロとなる。この相違は、簡略化されたモデルが磁場増幅と粒子の生存を過小評価しており、結果として超効率性を過小評価していることを浮き沢している。

2. 周波数依存のメカニズム：

   • FIRの超効率性： これは最も一般的な結果であり、長寿命の低エネルギー電子の蓄積によって駆動される。これらの粒子は、圧縮中に磁場が増幅される際に効率的に放射を行う。FIRにおける超効率性は、圧縮後フェーズまで持続することが多い。
   • X線の超効率性： これはより稀であり、かつ一時的である。これは能動的な圧縮フェーズと密接に関連している。これには、高エネルギー粒子を急速に枯渇させる致命的なシンクロトロン損失を引き起こすことなく、磁場増幅がシンクロトロン輝度を押し上げるというバランスが必要となる。
   • GeVの超効率性： 逆コンプトン（IC）散乱によって駆動されるこの現象は、蓄積された電子集団の長期的な生存を反映している。これは瞬間的な圧縮への依存度が低く、減少していく現在の $\dot{E}_{SD}$ に対する粒子注入の履歴に依存する。

3. 進化的傾向：

   • 超効率性は静的な特性ではなく、時間とともに進化する。集団がリバーブレーションフェーズに入るにつれて、超効率的なソースの割合は増加する。
   • 低周波数（FIR）の超効率性は、初期の圧縮後も長く持続するが、中間周波数（X線）の超効率性はより一時的であり、能動的な圧縮中にピークに達し、高エネルギー粒子が冷却されるにつれて減少する。
   • 本研究は、$\text{L}_0/\text{L}_{ch}$ 対 $\tau_0/t_{ch}$ 平面における2つの異なるパラメータ領域から超効率性が生じ得ることを特定した。すなわち、初期スピンダウンが低く注入期間が長いシステムと、初期スピンダウンが高かったものの大幅に減衰し、現在の $\dot{E}_{SD}$ を上回る「遺物（レリック）」の輝度を残しているシステムである。

4. 圧縮因子（CF）依存性：

   • リバーブレーションの強さ（CF）がスペクトルの特徴を決定的に制御する。緩やかな圧縮（CF < 10）は、深刻な冷却を引き起こすことなく電子の蓄積を許容するため、FIR超効率的ソースの数を最大化する。
   • 極端な圧縮（CF > 100）は、致命的なシンクロトロン冷却をもたらし、電子の貯蔵庫を枯渇させ、特にFIRにおける超効率的ソースの数を減少させる。

5. エネルギー等分（Equipartition）の安定性：

   • 粒子エネルギーと磁気エネルギーの比（$E_p/E_B$）の解析により、エネルギー等分は安定した進化の吸引子ではないことが明らかになった。システムは急速に等分状態から逸脱し、動的な再処理や逆ショックとの相互作用の間に一時的な交差が発生する。

意義
本論文は、集団合成フレームワーク内における、超効率的PWN候補を解釈するための最初の統計的ベースラインを確立するものである。超効率性は、異常な加速メカニズムを必要とする特異な現象ではなく、リバーブレーション物理学の自然かつ極端な帰結であることを示している。

本研究は、リバーブレーションフェーズの正確なモデリングが不可欠であると主張している。簡略化された薄いシェルモデルは、特にFIRおよびX線帯域における超効率性の真の普及度を捉えることができない。リバーブレーションによるエネルギー転送を無視すると、一部の高輝度なPWNが誤認されたり、その加速効率が過大評価されたりする可能性があることを示唆しており、本研究の結果は、ターゲットを絞ったマルチ波長追跡観測の指針となる。具体的には、FIRおよびGeV帯が、銀河系集団におけるこれらの現象を検出するための最も有望な領域であることを示している。