[論文レビュー] SPH Methods in the Modelling of Compact Objects

Stephan Rosswog|arXiv (Cornell University)|Jun 17, 2014

Gamma-ray bursts and supernovae参考文献 197被引用数 4

ひとこと要約

本稿は、中性子星、白色矮星、ブラックホールの合体をシミュレートするための滑らか化粒子流体力学（SPH）手法についてレビューする。質量・運動量・エネルギー・角運動量の保存性、柔軟性、複雑なミクロ物理学の実装を強調する。SPHは、短時間ガンマ線バースト、キロノーバ、重元素の核合成の理解に不可欠な、中性子星、白色矮星、ブラックホールの合体の高精度な3次元モデリングを可能にする。最近のシミュレーションでは、ニュートン近似、後ニュートン近似、共形平坦性近似の枠組みにおいて、高精度な結果が得られている。

ABSTRACT

We review the current status of compact object simulations that are based on the Smooth Particle Hydrodynamics (SPH) method. The first main part of this review is dedicated to SPH as a numerical method. We begin by discussing relevant kernel approximation techniques and discuss the performance of different kernel functions. Subsequently, we review a number of different SPH formulations of Newtonian, special- and general relativistic ideal fluid dynamics. We particularly point out recent developments that increase the accuracy of SPH with respect to commonly used techniques. The second main part of the review is dedicated to the application of SPH in compact object simulations. We discuss encounters between two white dwarfs, between two neutron stars and between a neutron star and a stellar-mass black hole. For each type of system, the main focus is on the more common, gravitational wave-driven binary mergers, but we also discuss dynamical collisions as they occur in dense stellar systems such as cores of globular clusters.

研究の動機と目的

• SPHがコンパクト天体の合体をシミュレートする数値的手法として現在どのように位置づけられているかを評価すること。
• ニュートン的、特殊相対論的、一般相対論的流体力学をモデル化する際の、さまざまなSPH定式化の精度と頑健性を評価すること。
• 白色矮星、中性子星、ブラックホールの二重星合体および動的衝突をシミュレートするSPHの応用を要約すること。
• SPHが、降着円盤内のニュートリノ冷却、核燃焼、粘性といった複雑なミクロ物理学をどのようにモデル化できるかを強調すること。
• 2014年時点で未解決の課題、特に完全に自己自己一貫的な動的一般相対論的SPHシミュレーションの欠如を特定すること。

提案手法

• 質量・運動量・エネルギー・角運動量の正確な保存を保証するため、非対称カーネル勾配を用いたSPHを採用する。
• SPH離散化における精度を向上させるために、さまざまなカーネル関数と近似手法を用いる。
• SPHシミュレーションにおいて、一般相対論のニュートン近似、後ニュートン近似、共形平坦性近似を適用する。
• 多圧縮的および現実的な状態方程式、ニュートリノ放出、降着円盤内のα粘性を含むミクロ物理学モデルを実装する。
• シャクーラ＝スニャエフ粘性によって駆動される、合体後の長時間スケールの円盤進化を2次元軸対称SPHシミュレーションでモデル化する。
• ブラックホールを、イベントホライズン境界条件を満たすニュートン的点質量として扱う。

実験結果

リサーチクエスチョン

• RQ1 異なるSPH定式化およびカーネル関数は、コンパクト天体シミュレーションの精度にどのように影響するか？
• RQ2 SPHは、二重系における角運動量輸送および質量移動をどの程度正確にモデル化できるか？
• RQ3 ニュートリノ冷却やレプトン数勾配といった、合体後の降着円盤を形作る主要な物理過程は何か？
• RQ4 中性子星–ブラックホールおよび中性子星–中性子星合体のSPHシミュレーションは、短時間ガンマ線バーストやキロノーバの観測的特徴とどの程度一致するか？
• RQ5 現在のSPH手法には、コンパクト天体合体における動的一般相対論的記述にどのような限界があるか？

主な発見

• SPHシミュレーションは、中性子星–ブラックホール合体の円盤で、約10^53 erg/sの高ニュートリノ放射度と、約0.4秒間にわたって合計10^52 ergのニュートリノエネルギー放出を正確に再現した。
• 円盤内にレプトン数勾配の逆転が発生し、遷移半径（約10^7 cm）近辺でY_e ≈ 0.05、ブラックホール近辺でY_e ≈ 0.10となる。この勾配が強力な対流運動を駆動する。
• 大きな半径部では陽子が豊富な状態（Y_e > 0.5）が見られ、r過程核合成の可能性がある場所であることが示唆された。
• SPHは、幾何学的柔軟性と保存性の両方を備えているため、コンパクト天体合体のような複雑で非対称的かつ高動的な系をモデル化するのに適していることが確認された。
• 進展は見られたが、2014年6月現在、完全に自己自己一貫的な動的一般相対論的SPHによるコンパクト天体合体シミュレーションは未解決の課題のままである。
• SPHは、短時間ガンマ線バースト、キロノーバ、重元素の合成を結びつける中心的役割を果たしており、特にGRB 130603BのnIRトランジット観測との整合性からその重要性が強調されている。