Chandra X-ray Observations of the Pulsar Wind Nebula within CTA 1

CTA 1 内のパルサー風星雲に対するチャンドラ X 線観測の技術的要約

問題と動機
パルサー風星雲（PWN）は、磁化された極相対論的風における粒子加速を理解するための重要な実験場を提供する。チャンドラ X 線天文台（CXO）による高解像度 X 線イメージングは、ジェットやトーラスといった形態的特徴を明らかにし、パルサーの自転軸と磁気軸の向きを制限してきたが、磁気圏放射理論におけるモデル依存性のため、観測角度（$\zeta$）と磁気傾斜角（$\alpha$）に対する独立した制限を得ることは依然として困難である。さらに、電波静穏なパルサー PSR J0007+7303 を含む複合超新星残骸（SNR）CTA 1 は、パルサー風、逆衝撃波、および周囲の媒質との相互作用が完全に特徴付けられていない複雑な系である。具体的には、パルサーの横方向速度、自転軸の 3 次元方位、および電波から PeV $\gamma$ 線までの広帯域スペクトルエネルギー分布（SED）は、深層 X 線データと現代の $\gamma$ 線解析を用いた更新された制限を必要とする。

手法
著者らは、2003 年のアーカイブ観測（エポック 1）1 回と 2023–2024 年の新規観測 6 回を合わせた、総露出時間約 194 ks の PSR J0007+7303 に対する深層チャンドラ観測を提示する。データは CIAO v4.15 を用いて再処理された。

1. 天体測位と固有運動：20 年間のベースラインにわたるパルサーの固有運動を測定するため、エポック間の X 線源を照合して相対天体測位を行った。
2. イメージングと形態：統合画像を解析し、構造成分（ジェット、反ジェット、トーラス）を同定した。放射強度の半径方向プロファイル分析により、観測された放射を MARX でシミュレートした点像関数（PSF）と比較し、トーラスの広がりを定量化した。トーラスの頂部への楕円フィットを用いて観測角度 $\zeta$ を推定した。
3. スペクトル解析：パルサー、トーラス、ジェット、反ジェット、および拡大星雲のスペクトルを、星間吸収（$tbabs$）で修正されたべき乗則モデルを用いて抽出し、Sherpa を用いてスペクトルフィットを行った。
4. $\gamma$ 線解析：Fermi-LAT データ（16 年分）を 50 GeV–1 TeV 範囲で解析し、空間モデル（点源対拡大円盤）をテストし、高エネルギースペクトルを導出した。
5. 広帯域モデリング：電波から PeV までのマルチウェーブ長 SED を、2 つのアプローチでモデル化した：
   • 指数関数的カットオフを伴う破れたべき乗則（ECBPL）電子分布を用いた、1 領域レプトンモデル（naima 使用）。
   • 角サイズ、電波の上限、X 線/$\gamma$ 線フラックスによって制限された、SNR と PWN の時間進化を追跡する 1 領域動力学進化モデル（Gelfand et al. 2009）。

主要な結果

• 形態：統合チャンドラ画像は、南へ伸びて南西へ曲がる約 20 秒角のジェット、北側の faint な反ジェット、およびジェット軸とほぼ垂直に配向したコンパクトなトーラスを明らかにした。ジェットは、ラム圧縮と一致する曲がり部分で明るくなるなど、顕著な表面輝度の変動を示す。
• 固有運動：距離 1.4 kpc において、パルサーの横方向速度は 95% 信頼区間で $\lesssim$200 km s$^{-1}$ と制限された。この上限値は、SNR 中心からのパルサーのオフセットから導かれた以前の推定値よりも著しく低く、SNR がより古い、あるいは爆発中心が以前想定されていたよりも現在のパルサー位置に近い可能性を示唆している。
• 幾何学：トーラスの観測角度は $\zeta \approx 50^\circ$ と決定された。パルスプロファイルモデリング（外部ギャップモデルと 2 極カウスティックモデル）と組み合わせると、磁気傾斜角 $\alpha \sim 20^\circ$–$70^\circ$ が示唆される。トーラスは南東側で明るく、これは接近側のドップラー増幅と一致する。
• スペクトル特性：
  • コンパクト領域（トーラス、ジェット）は、光子指数 $\Gamma \approx 1.2$–$1.4$ の硬いスペクトルを示し、放射冷却が最小限であり、加速スペクトルを直接観測していることを示している。
  • 拡大星雲は、シンクロトロン冷却と一致する軟らかいスペクトル（$\Gamma = 1.85 \pm 0.11$）を持つ。
  • コンパクト PWN の X 線放射効率は低い（$\sim 3.5 \times 10^{-5}$）。
• 広帯域 SED と動力学：
  • 1 領域モデリングは、低い磁場（$B \approx 1.4$–$3.2$ $\mu$G）と高い電子カットオフエネルギー（$E_{cut} \sim 0.2$–$0.3$ PeV）をもたらす。
  • 動力学進化モデル（モデル A：圧縮後；モデル B：圧縮前）の両方が、硬い低エネルギー指数（$p_1 \approx 0.7$）と急峻な高エネルギー指数（$p_2 \approx 2.7$–$2.9$）を持つ破れたべき乗則注入スペクトルを支持する。
  • モデル A（圧縮後）は、年齢 15.5 kyr および $B \approx 1.4$ $\mu$G で良好なフィット（$\chi^2_{red} = 1.65$）を提供し、PWN が約 2.7 kyr 前に逆衝撃波に遭遇したことを示唆する。
  • 低い磁場は、SNR 内でのパルサー風の急速な断熱膨張に起因すると考えられる。

意義と主張
著者らは、CTA 1 を PeV エネルギーまで粒子を加速可能な、若く X 線効率が低い PWN として確立した。硬い注入スペクトル（$\Gamma \approx 1.2$–$1.4$）と LHAASO による超高エネルギー（UHE）放射の検出は、縞状風領域における磁気リコネクションが支配的な粒子加速メカニズムであることを示唆する。低い放射効率と硬いスペクトルは、低い対多重性を示す可能性があり、これはパルサーの磁場（$1.1 \times 10^{13}$ G）と表面温度が最大対生成の範囲からわずかに外れていることと一致する。本研究は、特定の磁気圏放射モデルを支持する独立した幾何学的制限（$\zeta, \alpha$）を提供し、PWN の動力学状態（おそらく圧縮後）が導き出された磁場と粒子分布に大きく影響することを示している。