SN 1006: A Cosmic Laboratory for Investigating Shock Acceleration Physics

技術的要約：衝撃波加速物理学の宇宙的実験室としてのSN 1006

問題提起
超新星残骸（SNR）は銀河宇宙線の主要な供給源であると広く考えられているが、粒子加速、磁場増幅（MFA）、および$\gamma$線放射の起源（ハドロン成分かレプトン成分か）を支配する具体的なメカニズムについては、依然として活発な研究対象となっている。歴史的なIa型超新星残骸であるSN 1006は、その両極的な形態により、非熱的放射が周囲磁場に沿った2つの極（ポーラーキャップ）に集中していることから、ユニークな「宇宙的実験室」を提示している。この幾何学的構造により、様々な傾斜角（準平行ショック［極］対 準垂直ショック［非極領域］）における衝撃波加速の研究が可能となる。これまでの研究では、個々の象限やスペクトルエネルギー分布（SED）を個別にモデル化してきたが、全領域にわたってマルチ波長でのスペクトルおよび空間的特性を同時に再現し、かつ非線形拡散衝撃波加速（DSA）効果を考慮した自己整合的なモデルは欠けていた。

手法
著者らは、SN 1006をシミュレートするために、粒子加速の自己整合的なマルチゾーン・キネティックモデルを用いている。このフレームワークは、以下の主要な物理コンポーネントを統合している。

1. 衝撃波流体力学： モデルは、放出エネルギー $E_{SN} = 10^{51}$ erg、放出質量 $M_{ej} = 1.2 M_{\odot}$ を持つ球状のIa型SNRを仮定している。モデルは、放出物支配（ED）段階およびセドフ・テイラー（ST）段階を通じての進化を追跡する。周囲密度は、極（ポーラーキャップ）および南東（SE）領域については一様（$n_{ISM} \approx 0.02 \text{ cm}^{-3}$）として扱う一方、北西（NW）領域については、高密度の原子雲との最近の相互作用（$n_{cl} \approx 0.22 \text{ cm}^{-3}$）を考慮してステップ関数型の密度プロファイルを用いてモデル化している。
2. キネティック粒子加速： モデルの核となる部分は、非線形DSA理論を取り入れた、宇宙線（CR）輸送のための拡散・移流方程式を解くことである。ここでは、CR圧および乱流磁場によって修正された衝撃波ジャンプ条件を自己整合的に計算する。
   • 注入と効率： プロトンの注入量は衝撃波の傾斜角に基づいて調整され、準平行領域では高い加速効率（$\xi_{CR} \sim 21\%$）を、準垂直領域では低い効率（$\xi_{CR} < 1\%$）をもたらす。
   • 磁場増幅（MFA）： モデルには、ベル不安定性（非共鳴）および共鳴ストリーミング不安定性によって駆動される増幅が含まれる。また、磁気揺らぎが熱プラズマに対してドリフトすることで、圧縮比 $R > 4$ およびスペクトル指数 $q > 2$ を導く「ポストカースサー（後方追従）」効果を考慮している。
   • 拡散： X線カットオフエネルギーの方位角依存性によって制約される、ボーム拡散からより効率の低い拡散への遷移を説明するために、「ボーム因子」$\eta$ が導入されている。
3. 放射プロセス： 電子分布は、放射損失補正を用いてプロトンスペクトルから導出される。光子スペクトル（電波からマルチTeV $\gamma$ 線まで）は、シンクロトロン放射、逆コンプトン（IC）散乱、非熱的制動放射、および中性パイ中間子崩壊を含む naima パッケージを使用して生成される。
4. 空間モデリング： 非熱的粒子の不変質量殻を、断熱膨張と圧力平衡を仮定して進化させ、視線方向に沿って放射率を積分することにより、径方向プロファイルが構築される。

主な貢献と結果
本モデルは、極（ポーラーキャップ）およびNW領域の観測されたスペクトルエネルギー分布（SED）と空間プロファイルを正常に再現しており、以下の具体的な知見を得ている。

• 加速効率と傾斜角： 本研究は、準平行ショックが準垂直ショックよりも大幅に効率的にCRを加速することを裏付けている。極（ポーラーキャップ）は、$\sim 21\%$ の加速効率を示し、その結果、修正された衝撃波圧縮比 $R \approx 5.1$ と急峻なCRスペクトル $q \approx 2.2$ をもたらす。対照的に、準垂直領域は1%未満の効率と $R \approx 4.0$ を示す。
• 磁場増幅（MFA）： モデルは、極における下流磁場強度 $B_2 \approx 42 \mu\text{G}$（周囲の $3 \mu\text{G}$ に対して $\gtrsim 10$ 倍に増幅）を推定している。この単一のパラメータは、CRスペクトルの傾き、増幅された圧縮比、X線リムの厚さ、およびシンクロトロン放射の規格化という4つの独立した観測量を同時に説明できる。
• $\gamma$ 線放射の性質： 極からの$\gamma$線放射は主にレプトン的（逆コンプトン）であり、GeV帯におけるハドロン成分はわずかである。逆に、衝撃波が高密度雲と相互作用しているNW領域では、加速効率が低いにもかかわらず、ハドロン的放射（$\pi^0$ 崩壊）が支配的である。これは、「ハドロニシティ（ハドロン性）」が加速効率単独ではなく、周囲の密度によって駆動されるという仮説を支持している。
• 電波プロファイルの不一致： モデルはX線および$\gamma$線のプロファイルは良好に再現するが、観測された薄い電波リムの再現には失敗している。モデルは電波のピークを接触不連続面（CD）付近と予測するが、観測では前進衝撃波（FS）に近い位置にある。電波の幅を合わせるために磁場減衰を導入すると、SEDとの矛盾が生じる（相反する電子・陽子比を要求される）。著者らは、この不一致は、彼らの1次元流体力学近似では捉えきれない、3D流体力学効果（例：レイリー・テイラー不安定性や放出物の塊状構造）によってCDがFSに近づいていることに起因すると結論付けている。

意義
本論文は、最先端の非線形DSA理論およびキネティックシミュレーションの結果に基づいた自己整合的なキネティックモデルが、最小限の自由パラメータでSN 1006の複雑なマルチ波長形態を説明できることを示している。得られた結果は、以下の強力な観測的証拠を提供している。

1. 宇宙線加速効率の、衝撃波の傾斜角に対する連続的な依存性。
2. CRスペクトルと衝撃波圧縮を形成する上でのポストカースサー効果の決定的な役割。
3. $\gamma$ 線放射メカニズム（レプトン vs ハドロン）の環境依存性。
4. 電波形態を説明する上での1次元流体力学モデルの限界、および流体力学的不安定性がSNRシェルの構造において重要な役割を果たしていることの示唆。

本研究は、SN 1006が衝撃波加速の微視的物理学とSNR進化のマクロ物理学を理解するための主要なテストベッドであるという見解を強化し、高解像度の観測データに対して、キネティックシミュレーションから導かれた理論的予測を検証している。