[論文レビュー] Heavy elements and electromagnetic transients from neutron star mergers

Stephan Rosswog, Oleg Korobkin|arXiv (Cornell University)|Aug 30, 2022

Gamma-ray bursts and supernovae被引用数 4

ひとこと要約

本稿は、中性子星合体がr過程元素および電磁的過渡現象（kilonova）の源として果たす役割をレビューし、GW170817からの理論的・観測的知見を統合する。公開済みの加熱率ライブラリと、崩壊生成熱を再現するためのフィット式を提示しており、多様な噴出物質条件におけるkilonova光曲線の正確なモデリングを可能にする。

ABSTRACT

Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here we review our current understanding of the mass ejection mechanisms, the properties of the ejecta and the resulting radioactive transients. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, we also point out questions where no consensus has been reached yet, and we sketch possible directions for the future research. In an appendix, we describe a publicly available heating rate library based on the WinNet nuclear reaction network, and we provide a simple fit formula to alleviate the implementation in hydrodynamic simulations.

研究の動機と目的

• 中性子星合体からの質量噴出メカニズム、噴出物質の性質、それに続く電磁的過渡現象に関する現在の理解を統合すること。
• 特に希少で高収量な出来事の観点から、最も重いr過程元素を生成する天体的サイトが何かという未解決の問いに応えること。
• 流体力学的シミュレーションに放射性崩壊加熱を実装するための実用的ツール（フィット式と加熱率ライブラリ）を提供すること。
• r過程元素の生成パターンが電子割合 $Y_e$ にどのように依存するかを明確にすること、特に軽いr過程と重いr過程の間の遷移が $Y_e \approx 0.25$ で起こることを強調すること。
• 今後の研究における噴出物質チャネル、核合成生成物、観測的特徴に関する未解決の問題を特定すること。

提案手法

• 著者らは、WinNet核反応ネットワークを用いて、さまざまな噴出物質条件におけるr過程核の放射性崩壊からの加熱率を計算する。
• 加熱率 $\dot{\epsilon}(t)$ のパラメトリックフィット式を導出し、$Y_e$ と噴出速度 $v_{\rm ej}$ に依存する可変係数を持つ3つの指数関数の和としてモデル化する。
• 詳細な核合成シミュレーションを用いてフィットをキャリブレーションし、動的噴出、風、粘性噴出を含む複数の噴出物質設定で妥当性を検証する。
• 係数は $Y_e$ が 0.05 から 0.5 まで、$v_{\rm ej}$ が 0.05c から 0.5c までの範囲で表形式で提供され、流体力学コードへの直接実装を可能にする。
• 加熱率ライブラリは公開されており、合体シミュレーションにおけるkilonova光曲線の一貫性のあるモデリングを支援する。
• 本手法により、リアルタイムシミュレーションで完全な核反応ネットワークを解かずに、時間に依存するエネルギー供給を効率的に計算可能となる。

実験結果

リサーチクエスチョン

• RQ1 中性子星合体における主な質量噴出メカニズムは何か？ また、それらはr過程核合成にどのように影響を与えるか？
• RQ2 異なる噴出物質成分における電子割合 $Y_e$ は、軽いr過程と重いr過程元素の生成にどのように寄与するか？
• RQ3 噴出速度と $Y_e$ を変化させた場合、kilonovaの輝度の時間的変化にどのような定量的影響が生じるか？
• RQ4 流体力学的シミュレーションを効率的に行うために、放射性加熱率のシンプルで正確なパラメトリックフィット式をどのように構築できるか？
• RQ5 現在のkilonova放射のモデルにおける主な不確実性は何か？ また、合意と未解決の問いはどこに残っているか？

主な発見

• 中性子星合体は1回の出来事で太陽質量の約1%を噴出するが、異なる噴出物質成分における $Y_e$ の変動が顕著で、それにより明確に異なるr過程元素の生成パターンが生じる。
• 軽いr過程と重いr過程の間の遷移は $Y_e \approx 0.25$ で急激に起こり、$Y_e < 0.25$ の場合、$A \sim 195$（プラチナムピーク）までの元素が生成される。
• 加熱率 $\dot{\epsilon}(t)$ のパラメトリックフィット式は、多様な噴出物質条件において高い精度を達成しており、$Y_e \in [0.05, 0.5]$ および $v_{\rm ej} \in [0.05c, 0.5c]$ の範囲で係数が提供されている。
• 加熱率ライブラリにより、流体力学的シミュレーションへの実装が効率的に行えるようになり、計算コストを削減しながらもkilonova光曲線の予測精度を維持できる。
• GW170817からの観測的制約は、高 $Y_e$（軽いr過程）と低 $Y_e$（重いr過程）の両方の噴出物質成分が存在することを支持しており、マルチメッセンジャー観測データと整合的である。
• 進展は見られても、動的噴出、風、粘性噴出などの異なる噴出メカニズムの相対的寄与と、それらが核合成生成物に与える影響に関する不確実性は依然として残っている。