A Three-Dimensional Exploration of Magnetic Fields, Rotation, and Shock Revival in a $39 M_\odot$ Core-Collapse Supernova Progenitor

技術的概要：39 M⊙ 核崩壊型超新星の progenitor における磁場、回転、および衝撃波再生の 3 次元探求

問題と動機
核崩壊型超新星（CCSNe）における停滞した衝撃波の再生は、ニュートリノ加熱、回転、磁場、および流体力学の相互作用を含む複雑な問題である。ニュートリノ機構は多くの progenitor、特に高コンパクトネスを持つものにおいて堅牢であるが、磁気回転（MR）機構は、衝撃波再生と高エネルギー遷移現象（例えば、ハイパーノヴァ、Type Ic-BL 超新星）の駆動のための代替経路を提供する。この機構は、回転自由エネルギーをねじれや不安定性を通じて磁気エネルギーに変換し、双極性アウトフローを駆動するために、急速な回転と強力な磁場を必要とする。しかし、過去の研究はしばしば軸対称（2D）に依存しており、それはジェットを人工的に安定化させ、キント不安定性のような非軸対称不安定性を抑制する可能性がある。さらに、高コンパクトネスかつ急速に回転する 3 次元環境において、磁場が衝撃波再生の時間スケール、アウトフローの形態、およびブラックホール形成に対する原始中性子星（PNS）の安定性に与える具体的な影響については、さらなる調査が必要である。

手法
著者らは、高コンパクトネス（$\xi_{2.5} = 0.75$）および急速な回転（$\Omega_c \approx 0.58$ rad s$^{-1}$）を持つ 39 M⊙ 零歳主系列（ZAMS）progenitor の 3 次元流体力学（HD）および磁気流体力学（MHD）シミュレーションを実施した。Progenitor は、準化学的に均一な進化を用いて MESA で進化させた。3 つの異なるモデルを、適応的メッシュ細分化（AMR）および M1 閉鎖を用いた多次元・エネルギー依存型 2 モーメントニュートリノ輸送スキームを採用した Spark ソルバーを備えた FLASH コードでシミュレートした：

1. B1R1： 回転し、磁化されたモデル。回転軸と整列した極性セeding磁場（$B_0 \approx 5 \times 10^{11}$ G）を初期化。
2. B0R1： 回転するが非磁化の対照モデル（HD のみ）。
3. B0R0： 非回転かつ非磁化の基準モデル（HD のみ）。

シミュレーションは、崩壊前からコアバウンスを経て、初期のバウンス後段階（約 1 秒まで）へと進化し、コア付近で最高約 0.61 km の分解能を持つ直交格子を利用した。一般相対論的効果は、有効ポテンシャル（Marek et al. 2006 の Case A）を用いて近似された。

主要な結果
3 つのモデルすべてがシミュレーション期間内で暴走する衝撃波膨張を達成したが、形態、時間スケール、および残骸の運命において著しく異なっていた：

• 衝撃波再生と形態：

  • B1R1（磁化/回転）： 最も早い衝撃波再生（バウンス後 $t_{rev} \approx 150$ ms）を示し、最も明確な双極性アウトフローを発達させた。アウトフローは最初はジェット状であったが、非軸対称不安定性によって広がり、歪んだ。
  • B0R0（非回転）： 強力なニュートリノ駆動乱流を介して遅く（$t_{rev} \approx 250$ ms）再生した。衝撃波は、好ましい極方向を持たず、大規模で無秩序なプラムとして膨張した。重要なのは、このモデルが回転支持の欠如により、バウンス後約 1 秒後にブラックホールへ崩壊した点である。
  • B0R1（回転/非磁化）： 最も遅い衝撃波再生（$t_{rev} \approx 370$ ms）を経験した。回転は対流を抑制しニュートリノ光度を低下させ、B0R0 に比べてよりコンパクトな進化をもたらした。明確な双極性アウトフローは発達しなかった。

• 角運動量とエネルギー輸送：

  • B1R1 において、マクスウェル応力が効率的に角運動量を内部 PNS から外側へ再分配し、内部コアの自転を減速させ、回転自由エネルギーを抽出した。この過程は、極方向のアウトフローを駆動するために微分回転の貯蔵庫を利用した。
  • 角運動量の予算は、移流が内部コアへ角運動量を供給する一方で、磁気トルクが主要な吸収源となり、スピンの増加を相殺したことを示した。ニュートリノ放射は二次的かつ小さな吸収源として作用した。
  • B1R1 の診断エネルギー（$E_{diag}$）は、回転エネルギーの抽出と一致して、最も早くかつ最も急速に増加した。磁気エネルギー予算は、抽出されたエネルギーの相当部分が、解像された磁気エネルギーとして保持されるのではなく、おそらく未解決の再結合およびサブグリッド移動に関連する実効的な散逸項（$\dot{Q}_{diss}$）を通過したことを示した。

• 不安定性と 3 次元力学：

  • B1R1 の磁化されたアウトフローは、非軸対称キント不安定性（$m=1$ モード）の影響を受けた。診断は、磁気エネルギー重心の変位の早期の指数関数的成長（e 倍時間 $\sim$1.5 ms）を示し、理論的なキント時間スケールと一致した。
  • 安定した狭いジェットを予測する可能性がある軸対称シミュレーションとは異なり、3 次元シミュレーションでは、アウトフローの漏斗が時間とともに数度傾き、広がることが示された。アウトフローは生存し、大半径まで伝播したが、その狭いコリメーションを失い、ジェット安定性を正確にモデル化するには 3 次元効果が決定的であることを示唆した。

• 核合成：

  • 放出物のポストプロセッシングは、明確な組成の違いを示した。B0R0 はニュートリノ駆動爆発に典型的な陽子過剰な放出物（$Y_e > 0.5$）を生成し、鉄族元素（$A \sim 56$）まで合成した。
  • B1R1 は、急速なアウトフローの発達により、$Y_e$ が $\sim 0.33$ まで低下する顕著な中性子過剰成分を生成し、$A \sim 90$ までの弱い r 過程を可能にした。しかし、初期のセeding磁場とキント誘起による広がりは、より極端な磁場を持つモデルで見られる極めて中性子過剰な放出物（$Y_e < 0.32$）の生成を防いだ。

意義と主張
本論文は、極端な高コンパクトネス progenitor において、急速な回転と強力な磁場の組み合わせが、本質的に 3 次元的な早期の磁場支援衝撃波再生経路を開くことを主張する。主要な結論は以下の通りである：

1. 駆動源としての磁場： 磁気応力は単にアウトフローを形成するだけでなく、回転自由エネルギーを積極的に抽出し、内部コアの自転を減速させ、非磁化回転の場合よりも早く衝撃波再生を促進する。
2. 3 次元の必要性： 磁場駆動アウトフローの進化は本質的に非軸対称である。キント不安定性はアウトフローを拡大・傾斜させ、持続的な狭いジェットの形成を防ぐ。軸対称シミュレーションは、そのようなアウトフローの安定性とコリメーションを過大評価する可能性がある。
3. 残骸の安定性： 回転は、高コンパクトネス progenitor における即座のブラックホール形成に対して重要な支持を提供するが、残骸の長期的な安定性はシミュレーション期間を超えて依然として不確かである。
4. エネルギー経路： 微分回転からアウトフローへのエネルギー転移には、重要な未解決の散逸が関与しており、MHD CCSN シミュレーションにおけるサブグリッド物理学の重要性を浮き彫りにしている。

著者らは、未解決の小規模ダイナモ増幅を補うための強力なセeding磁場の使用、散逸診断における分解能の感度、および回転モデルの究極の運命を決定するためのより長いシミュレーションの必要性を含む限界を指摘している。彼らは、シミュレーションが磁気回転衝撃波再生の図を支持する一方で、具体的な結果はアウトフローの完全な 3 次元力学に大きく依存すると強調している。