Long-Term Multidimensional Models of Core-Collapse Supernovae: Progress and Challenges

過去 10 年間の 3 次元シミュレーションの飛躍的進展によりニュートリノ駆動爆発メカニズムの支持や多様な予測が可能になったものの、ブラックホール形成の条件やニュートリノ振動、核状態方程式などの不確実性により、依然としてどのような星が爆発するかを決定づけるにはさらなる研究と課題の克服が必要である。

要約

星の最期の爆発：10 年間の「3 次元シミュレーション」の進歩と課題

～超新星爆発の謎を解く、壮大な宇宙のドラマ～

この論文は、天体物理学の巨人であるハンス＝トーマス・ジャンカ氏によって書かれた、「大質量星が最期に爆発する（超新星爆発）」現象についての最新の総説です。

過去 10 年間で、コンピュータの性能が飛躍的に向上し、天文学者たちは「星の爆発」を3 次元（3D）で、かつ長時間にわたってシミュレーション（計算による再現）できるようになりました。まるで、星の心臓部で起きていることを、高解像度の 3D 映画のように観られるようになったのです。

この論文の核心を、難しい数式を使わずに、日常の比喩を使って解説します。

---

1. 爆発のトリガー：「風船」と「ポンプ」の物語

星が死を迎えるとき、その中心部は重力に耐えきれず崩れ落ちます（コア・コラプス）。すると、中心には超高密度の「中性子星」という小さなボールが生まれます。

• 昔の考え方（1 次元・2 次元）：
  星の中心から「ニュートリノ」という目に見えない粒子が大量に放出されます。これを「熱風」や「ポンプ」だと思ってください。昔のシミュレーションでは、このポンプが風船（星の外層）を膨らませる様子を、単純な「円筒形」や「平面的な図」で見ていました。
• 今の発見（3 次元）：
  最新の 3D シミュレーションでは、この「ポンプ」が働く様子がもっと複雑でダイナミックであることがわかりました。
  • お風呂の泡： 中心から出る熱いニュートリノが、冷たい物質を押し上げ、大きな「泡（対流）」を作ります。
  • 揺れるシャワー： 衝撃波が止まってしまうと、物質が中心に落ち込み、それがまた跳ね返るという「揺れ（SASI：立着降着衝撃波不安定）」が起きます。
    これらの「泡」や「揺れ」が、爆発を成功させるための**「コンベアベルト」**の役割を果たしているのです。

2. なぜ「長時間」のシミュレーションが必要なのか？

爆発は、星が崩れ始めてから数秒で始まると思われがちですが、実際には**「数秒から数十秒」**にわたってエネルギーが蓄積されていきます。

• 料理の例え：
  爆発を「大きな鍋でスープを煮込む」ことに例えてみましょう。

  • 短時間シミュレーション： 火を点けてすぐの「沸騰し始めた瞬間」しか見ていません。
  • 長時間シミュレーション： 何時間も煮込んで、味がどう染み込み、具材がどう動き、最終的にどう盛り付けられるかまで追跡します。

  この論文では、爆発が完全に決まるまでの「長い煮込み時間」を追跡することで、爆発のエネルギーがどう決まるか、中性子星がどう動くか、そしてどんな元素（金や鉄など）が作られるかが詳しくわかってきました。

3. 驚きの発見：爆発しない星もあれば、ブラックホールになる星も

すべての星がきれいに爆発するわけではありません。

• 爆発の「勝敗」：
  星の重さ（質量）によって、爆発するか、あるいは中心が潰れて「ブラックホール」になるかが決まります。しかし、「どの星が爆発するか」については、使っている計算ソフト（コード）によって答えがバラバラです。

  • あるソフトでは「爆発する」と言っても、別のソフトでは「爆発しない」と言ったりします。
  • これは、星の内部の「核の硬さ（状態方程式）」や、ニュートリノの動きの計算方法が微妙に違うためです。まだ「正解」が一つに定まっていない、まさに「白熱した議論」の最中です。

• ブラックホールの誕生：
  爆発しても、その後に落ち込んだ物質が重すぎて、中性子星がさらに潰れてブラックホールになる「遅れ爆発（フォールバック）」という現象も確認されました。まるで、爆発した花火の残骸が、重力に引かれて再び地面に吸い込まれていくような現象です。

4. 中性子星の「蹴り」と「回転」

爆発が起きると、生まれたばかりの中性子星は、まるでサッカー選手に蹴られたボールのように、ものすごい速度で宇宙を飛び出します（キック）。

• なぜ飛ぶのか？
  爆発が「左右非対称」に起きるからです。一方の方向にガスを強く吹き飛ばすと、反動で中性子星は逆方向に弾き飛ばされます。
• 回転（スピン）：
  渦を巻くように物質が落ち込むことで、中性子星は高速で回転し始めます。これは、氷上でスピンするフィギュアスケート選手が腕を閉じると回るスピードが上がるのと同じ原理です。

5. 1987A 年のミステリー：「最後の 3 秒」

1987 年に観測された「超新星 1987A」は、人類が初めてニュートリノ（素粒子）として捉えた星の爆発です。

• 謎： 観測データには、爆発から**「10 秒以上経ってから」**もニュートリノが検出されたという「最後の 3 つの信号」がありました。
• 現在のシミュレーション： 多くのモデルでは、中性子星の冷却が速すぎて、10 秒以上ニュートリノを出し続けることができません。
• 解決策の候補：
  1. 中性子星の内部構造（核の状態方程式）が、今のモデルより「柔らかい」かもしれない。
  2. 爆発後に物質が再び落ち込み（フォールバック）、エネルギーを再供給したのかもしれない。
  3. 物質がクォークという別の状態に変わる「相転移」が起きて、エネルギーが放出されたのかもしれない。

この「最後の 3 秒」の謎を解くことが、今後の大きな課題です。

6. 重力波：宇宙の「震動」

爆発の瞬間、時空そのものが揺れます。これを「重力波」と呼びます。

• 比喩： 池に石を投げた時の波紋のようなものですが、宇宙規模です。
• 3D シミュレーションは、この波紋が「どんなリズムで、どのくらいの強さで」発生するかを予測しています。将来的には、この重力波を聞くだけで、「どんな星が爆発したか」「どんな元素が作られたか」がわかるようになるかもしれません。

---

まとめ：まだ見えない未来へ

この論文は、**「私たちは星の爆発を、3D 映画のように詳しく再現できるようになったが、まだ『なぜ爆発するのか』の完全な答えは出ていない」**と伝えています。

• 成功： 爆発のメカニズム（ニュートリノ加熱）は、3D 効果によって支えられていることがほぼ確実になりました。
• 課題： 計算ソフトによって結果が異なる、元素の生成量や爆発のエネルギーにばらつきがある、ブラックホールになる条件がまだ不明確など、解決すべき問題が山積みです。

しかし、スーパーコンピュータの進化と、ニュートリノや重力波の観測技術の向上により、**「次の銀河系内の超新星爆発」**が起きたとき、私たちはその瞬間を「ニュートリノ」「重力波」「光」の 3 つの目で見ることができるようになります。それが、この研究の究極のゴールであり、宇宙の秘密を解き明かすための「最終試験」になるでしょう。

一言で言えば：

「星の爆発という壮大なドラマを、3D で長時間観られるようになったが、まだ脚本（物理法則）の細かい部分に謎が多く、次回の『実演（観測）』を待ちながら、計算機という舞台で練習を続けている最中だ」ということです。

技術概要

この論文は、ハンス＝トーマス・ジャンカ（Hans-Thomas Janka）氏による、過去 10 年間の多次数元（特に 3 次元）におけるコア崩壊型超新星（CCSN）の長期シミュレーションの進展と課題に関する包括的なレビューです。以下に、問題設定、手法、主要な貢献、結果、および意義について技術的に詳細に要約します。

1. 問題設定 (Problem)

コア崩壊型超新星の爆発メカニズムの解明は、天体物理学における最も重要な未解決問題の一つです。

• ニュートリノ加熱メカニズムの限界: 球対称（1D）モデルでは爆発が起きにくいことが知られており、2 次元（2D）および 3 次元（3D）の非対称性（対流、SASI 不安定など）が爆発を助けることが示されています。しかし、異なるシミュレーションコード間や、異なる初期条件（前崩壊段階の星の構造）において、どの星が爆発し、どの星がブラックホール（BH）を形成するかという「爆発可能性（Explodability）」の予測に大きな不一致が見られます。
• 長期進化の欠如: 従来のシミュレーションは、コア崩壊から数秒程度までしか追跡できていませんでした。しかし、爆発エネルギーの飽和、核合成の完了、コンパクト天体（中性子星 NS や BH）の誕生時の特性（キック、スピン）の決定には、爆発発生後 10 秒以上、場合によっては数十秒にわたる「長期進化」の理解が不可欠です。
• 観測との整合性: SN 1987A のニュートリノ観測データ（特に最後の数イベント）や、カシオペヤ座 A などの残骸の非対称性、中性子星のキック速度などを説明するために、より精密な長期シミュレーションが必要とされています。

2. 手法 (Methodology)

本レビューは、主にドイツ・ガルヒングの Max Planck 研究所で開発されたコード群、特にPrometheus-VertexおよびPrometheus-Nemesisを用いた長期シミュレーションに焦点を当てています。

• コードと物理モデル:
  • Prometheus-Vertex: 高精度な流体力学コード（Prometheus）と、多群ニュートリノ輸送コード（Vertex）を結合したもの。ニュートリノ輸送には「Ray-by-ray-plus (RbR+)」近似を用い、一般相対論的効果、多体効果（核子相関）、ニュートリノ反応率を詳細に扱います。
  • Prometheus-Nemesis: 長期シミュレーション（10 秒以上）における計算コストの壁を突破するために開発された手法。ニュートリノ輸送を 3D で直接計算するのではなく、1D 中性子星冷却シミュレーション（MLT 対流モデルを含む）の結果を参照し、ニュートリノ加熱率を近似する「ニュートリノ外挿法」を採用します。これにより、計算資源を大幅に削減しつつ、ニュートリノ加熱による爆発の長期進化を追跡可能です。
• シミュレーション範囲:
  • 前崩壊段階（対流酸素殻燃焼の 3D 摂動を含む）から、コア崩壊、バウンス、衝撃波の再生、衝撃波の星表面からの脱出、そして残骸段階に至るまでをシームレスに追跡。
  • 様々な質量（9 M⊙〜100 M⊙）の progenitor（前駆星）モデル、異なる核物質の状態方程式（EoS: LS220, SFHo, DD2, APR など）を用いたパラメータ研究。

3. 主要な貢献と結果 (Key Contributions & Results)

A. 爆発エネルギーの飽和と時間進化

• 爆発エネルギーは爆発発生直後に決まるのではなく、ニュートリノ加熱による物質の流入・流出サイクルが数秒〜10 秒以上続くことで徐々に増加し、最終値に達します。
• 高質量鉄核 progenitor では、長時間にわたる大量の降着流（accretion downflows）が PNS（原始中性子星）近傍に存在し、ニュートリノ加熱を受けて再び噴出（outflow）します。このプロセスが爆発エネルギーを大幅に増幅させます。
• 低質量 progenitor（電子捕獲型超新星など）では、爆発直後にニュートリノ駆動風（NDW）が発生し、エネルギーが早期に飽和します。

B. 中性子星・ブラックホールのキックとスピン

• キック速度: 非対称な質量放出（流体力学的キック）と非対称なニュートリノ放射（ニュートリノキック）が、新生コンパクト天体にキックを与えます。
  • 低質量爆発ではニュートリノキック（LESA 現象など）が支配的ですが、高質量 progenitor では非対称な質量放出による流体力学的キックが支配的となり、1000 km/s 以上の速度に達する可能性があります。
  • キックは爆発発生後 10 秒以上かけて漸近値に達するため、長期シミュレーションが必須です。
• スピン: SASI の螺旋モードや非対称降着により、中性子星は爆発開始時に 10ms〜数秒の周期で回転を得ますが、その後の降着物質の角運動量により変化します。

C. 核合成（ニュートリノ加熱流出物）

• 3D 流体力学シミュレーションでは、1D 球対称モデルの「ニュートリノ駆動風（NDW）」とは異なり、降着流と流出物が衝突する高度に乱流な環境が形成されます。
• これにより、エントロピーや電子分率（Ye）の分布が広がり、二次的な衝撃波による再加熱が発生します。
• 44Ti の生成: 1D モデルでは過小評価されていた 44Ti（チタン 44）の生成が、3D 長期シミュレーションでは長時間（〜10 秒）にわたる高エントロピー流出物によって効率的に説明可能となりました。これは SN 1987A や Cas A の観測データと整合します。
• 56Ni の生成: 異なるコード間（Fornax vs Prometheus-Vertex）で 56Ni の生成量に差が見られ、これは爆発エネルギーや progenitor 構造のモデル依存性を示唆しています。

D. ニュートリノ信号と SN 1987A

• 長期シミュレーションによるニュートリノ光度曲線は、SN 1987A の観測データ（Kam-II, IMB, BUST）の全体的な特徴（総エネルギー、平均エネルギー）とよく一致します。
• 課題: 観測された最後の数イベント（Kam-II の 9.2s〜12.4s）は、現在のモデルが予測する対流による冷却時間の短さ（〜10 秒未満）では説明が困難です。
  • 解決策として、降着によるニュートリノ再放射（fallback accretion）の存在、あるいは核物質の状態方程式（EoS）の相転移（ハドロン→クォーク）による再加熱などが提案されています。

E. ブラックホール形成の多様性

• ブラックホール形成は「失敗した爆発」だけでなく、爆発を伴う場合も存在します。Burrows らは 4 つの形成チャネルを提案しました。
  1. 爆発後に降着で BH 化（強い爆発）。
  2. 弱い爆発後に降着で BH 化（fallback SN）。
  3. 脈動対不安定（PPI）を経た巨大星の BH 化（爆発なし、または弱い）。
  4. 直接 BH 化（失敗した爆発）。
• どのチャネルが実現するかは、EoS、降着率、ニュートリノ加熱効率に強く依存します。

4. 意義と今後の課題 (Significance & Future Issues)

• 理論的意義: 3D 長期シミュレーションは、ニュートリノ駆動メカニズムが観測的な超新星エネルギーや残骸の特性を説明できることを強く支持しています。また、1D モデルでは捉えきれなかった乱流混合や非対称性の重要性を定量的に示しました。
• 観測との架け橋: 重力波（GW）信号やニュートリノ信号の予測精度が向上し、将来の銀河系内超新星観測（多メッセンジャー天文学）に対する理論的基盤を提供しています。
• 残された課題:
  • コード間の不一致: 異なるコード（Fornax, Prometheus-Vertex など）間で爆発の成否や特性に差が生じており、統一的な理解のためには厳密なコード比較（Code Comparison）プロジェクトが必要です。
  • ニュートリノフレーバー変換: 超新星コア内での「高速フレーバー変換（Fast Flavor Conversion）」の影響を、大規模シミュレーションにどう組み込むかが重要な未解決問題です。
  • 回転と磁場: 回転や磁場の効果を、一貫した進化計算から得られた前駆星モデルに組み込む必要があります。
  • EoS の不確実性: 高温高密度核物質の状態方程式（特に PNS 内部）が爆発の成否やニュートリノ信号に決定的な影響を与えるため、理論・実験・観測を統合した制約が必要です。

総じて、本論文は、超新星爆発の理解が「爆発のトリガー」から「爆発後の長期進化と多メッセンジャー信号の統合的理解」へとパラダイムシフトしたことを示し、計算天体物理学の最前線における重要なマイルストーンを記述しています。