Multi-Dimensional MHD simulations of young Core-Collapse Supernova Remnants

技術的サマリー：若年コア崩壊型超新星残骸の多次元MHDシミュレーション

問題と動機
若年コア崩壊型超新星残骸（SNR）は、特にペタ電子ボルト（PeV）エネルギーにおける宇宙線（CR）スペクトルの「膝」構造に関して、銀河宇宙線の加速の主要な候補である。これらのシステムのモデリングにおける重要な不確実性は、爆発前の進化星による質量放出によって形成される星周物質（CSM）にある。過去の研究では、定常状態の静的な風と平均質量放出率（$\dot{M}$）および終端速度（$v_\infty$）を仮定した、一次元（1D）または簡略化された多次元モデルに依存することが多かった。これらの簡略化は、赤色超巨星（RSG）対ウォルフ・ライエ（WR）相など、異なる進化段階間の相互作用によって生じる高密度殻、希薄なバブル、多重衝撃波など、時間依存型の恒星進化に起因する複雑なCSM構造を捉え損なう。さらに、粒子加速に不可欠な磁場配置は、恒星の自転と風力学に依存するが、既存のモデルではこれらが流体力学的進化から切り離されることが多い。

手法
著者らは、MHDコードpionを用いて、詳細な恒星進化軌跡から自己整合的に生成されたCSM中へ膨張する若年SNRの2次元円筒座標および3次元直交座標シミュレーションを行った。

• 恒星進化: CSMは、広域連星系における$31.6 M_\odot$の主星（単星として進化）についてMESAで計算された1次元恒星進化軌跡から導出された。軌跡は主系列（MS）、RSG、およびWR相をカバーする。質量放出則は相によって異なる（MSにはVinkら、RSGにはNieuwenhuijzen & de Jager、WRにはNugis & LamersおよびYoon）。
• MHD設定: シミュレーションには、理想MHD、静的メッシュ細分化、および放射加熱・冷却が含まれる。重要な仮定として、風と星間物質（ISM）の完全な光電離（$\sim 8300$ Kへの加熱）があり、これにより衝撃波のマッハ数が低下し、殻の形成が変化する。恒星磁場は分裂モノポールとして注入され、自転によって螺旋状に巻き上げられる。著者らは、銀河系O型星の観測と整合する、控えめな表面磁場（$B_\star = 10$ G）と低い初期自転率（臨界値の20%）を仮定した。
• シナリオ: 2つの標準的な進化星シナリオがモデル化された。
  1. RSG: ISM中へ膨張する遅く高密度の風。
  2. WR: 先行するRSG物質を掃き集め、複雑で構造化されたCSMを形成する、高速かつ高質量放出の風。
• 超新星の実装: 人工的な対称性を破るため、ガウス分布の摂動（$\pm 40\%$）を付与した球対称爆発（$5 M_\odot$、$10^{51}$ erg）を、超新星爆発前のCSM中に挿入した。放出物の磁場は、周囲のCSMのプラズマ$\beta$に一致するように設定された。

主要な結果

1. CSM構造と光電離: 光電離の仮定は、以前の非電離モデルと比較してCSM構造を著しく変える。
   • RSGの場合: 風バブルは、風終端衝撃波（WTS）において薄く高密度の殻を形成しない。代わりに、光電離された風とISMにより、$\sim 3$ pcにわたって密度が増加する広域領域が生じ、ISM中への亜音速膨張となる。
   • WRの場合: 高速風がRSG物質を掃き集め、高密度で放射性、かつレイリー・テラー不安定な殻を形成する。WR風は、SNRがこの殻と相互作用する際に、一貫した高速反射衝撃波群を生成する。
2. 衝撃波力学:
   • 前方衝撃波: 両シナリオにおいて、前方衝撃波の速度は、静的CSMを仮定する解析理論（例：Truelove & McKee 1999）の予測よりも速い。RSGの場合、これは光電離による熱圧力からの追加的な風加速に起因する。WRの場合、これはCSM自体の膨張速度に起因する。
   • 反射衝撃波: WRシミュレーションでは、掃き集められたRSG殻との相互作用により、反射衝撃波が発生し、内側へ伝播して逆衝撃波と相互作用し、 significantなエネルギー（$\sim 3 \times 10^{50}$ erg）を処理する。
3. 磁場:
   • 磁場強度は、圧縮・冷却領域（殻）において数十$\mu$Gのピーク値に達する。
   • Axford-Cranfill効果（方位角磁場の蓄積$B_\phi \propto r$）は顕著に再現されない。恒星自転率が低い（$v_{rot}/v_{crit} \sim 1\%$）ため、磁場はWTSの下流で依然として放射方向に支配的である。
   • 磁場構造は高度に非対称であり、特に固有運動効果によりRSGの場合に顕著である。
4. 粒子加速の可能性:
   • 前方衝撃波: 爆発後最初の10年間の最大粒子エネルギーは、RSGで$\sim 100$ TeV、WRで$\sim 30$ TeVと推定される。これらの値は、宇宙線の膝を説明するために必要なPeVエネルギーよりも著しく低い。
   • 反射衝撃波: WRの場合の反射衝撃波は大量のエネルギーを処理するが、内部磁場は弱い（$B \ll \mu$G）。その結果、最大電子エネルギーは数十TeVに制限される。ただし、著者らは、これらの電子がSNR内に閉じ込められている場合、逆コンプトン散乱による検出可能なガンマ線放射を予測している。

意義と結論
本論文は、詳細な多次元MHD処理と自己整合的な恒星進化を結合することが、SNRの初期進化を正確にモデル化するために不可欠であることを示している。

• 解析近似を超えて: 定常状態の静的CSMを仮定する解析モデルは、前方衝撃波の速度を過小評価し、時間依存型の質量放出と光電離に起因する複雑な密度構造（例、薄殻対拡張密度勾配）を捉え損なう。
• 恒星自転と磁気: 進化星モデルにおける急速自転の欠如は、WTS下流で顕著に蓄積しない放射方向に支配的な磁場をもたらし、孤立したSNRにおけるPeV粒子加速の可能性を制限する。
• PeVatron: 結果は、標準的なRSGまたはWR進化星に由来する孤立した若年SNRが、陽子をPeVエネルギーまで加速する可能性は低いことを示唆している。そのようなエネルギーを達成するには、高密度の星団内での超新星、スーパーバブル内での発生、または爆発直前に放出された高密度物質との相互作用など、極端な条件が必要である。
• 観測的シグネチャ: 本研究は、WR進化星における反射衝撃波がPeV加速器ではないものの、逆コンプトン散乱を介して検出可能なガンマ線放射を生成する可能性を強調しており、これらの特定の進化チャネルに対する潜在的な観測的シグネチャを提供する。

著者らは、若年SNRの衝撃波特性と非熱放射を正しく予測するためには、将来のモデルがCSMの動的な不安定性と時間依存性を考慮しなければならないと結論づけている。