[論文レビュー] Chemical Distribution of the Dynamical Ejecta in the Neutron Star Merger GW170817

Shreya Anand, P. T. H. Pang+12 more|arXiv (Cornell University)|Jul 20, 2023

Gamma-ray bursts and supernovaePhysics and Astronomy被引用数 3

ひとこと要約

本研究では、正確なハッブル宇宙望遠鏡およびVLBIの傾き制約（19–25°）を用いて、GW170817のkilonova AT2017gfoのモデルを精緻化し、観測者の視線がランタニドが乏しい極域の力学的噴出物質成分と交差していることが判明した。平均電子分率 $ar{Y}_{\rm e,dyn} = 0.22$ である。結果は、噴出物質の組成に顕著な方位依存性があることを確認し、ランタニドが豊富な物質は赤道面の周りに35°の半開口角で閉じ込められている。これは、従来の数値相対論的シミュレーションとの矛盾を解消し、第3ピークまでにまで強固な $r$-過程核合成が行われていることを支持する。

ABSTRACT

GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{ m e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{ m e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{ m e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{ m e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{ m e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{ m rd}$ peak in the tidal dynamical ejecta.

研究の動機と目的

• ハッブル宇宙望遠鏡およびVLBI観測から得られる正確な傾き角制約を組み込むことで、中性子星合体GW170817の噴出物質性質の推論を改善すること。
• 一次元kilonovaモデルと数値相対論的シミュレーションとの間にある不一致を、多次元的かつ傾き依存性のある噴出物質モデリングを用いることで解消すること。
• AT2017gfoの初期青色発光が、追加の駆動源を必要としない非球形で組成が分層した力学的噴出物質によって説明可能かどうかを検証すること。
• 力学的噴出物質における推定された電子分率分布に基づき、第3ピークまでにまで $r$-過程核合成がどれほど強固に行われているかを評価すること。
• 初期時系列光度測定と傾き制約を組み合わせることで、kilonovaの化学組成と噴出物質幾何学を信頼性高く推論できるかどうかを示すこと。

提案手法

• 最新の加熱率、熱化効率、周波数依存性を持つ透過率を組み込んだグリッドベースの多次元kilonovaモデリングフレームワークを採用した。
• 数値相対論的シミュレーションに基づき、角度依存性を仮定した $\bar{Y}_{\rm e,dyn}$（平均電子分率）を用いて力学的噴出物質をパラメータ化した。$Y_{\rm e,dyn}$ は赤道部で約0.17、極部で約0.41の範囲で変化するとされた。
• 合成光曲線でトレーニングされた順方向型ニューラルネットワークを用いて、全パラメータ空間における高速なベイズ的パラメータ推定を実現した。
• AT2017gfoにおける相対論的ジェットの固有運動測定から得られた、視線角のきつい事前分布（$19^\circ < \theta_{\rm obs} < 25^\circ$）を適用した。
• AT2017gfoの多バンド光度測定にモデルをフィットさせ、噴出物質質量、速度、電子分率（力学的およびディスク・ウィンド成分の両方）を推定した。
• 結果を数値相対論的予測と比較し、最近のSr II分光特徴に基づく球形噴出物質の主張とも照合した。

実験結果

リサーチクエスチョン

• RQ1 正確な傾き制約を組み込むことで、GW170817のkilonova AT2017gfoにおける力学的噴出物質の推定特性が顕著に変化するか？
• RQ2 AT2017gfoの初期青色発光は、追加の駆動源を必要としない非球形で組成が分層した力学的噴出物質によって説明可能か？
• RQ3 力学的噴出物質における電子分率の方位依存性は何か？ そしてそれは $r$-過程元素の生成にどのように影響するか？
• RQ4 推定された噴出物質パラメータは、数値相対論的シミュレーションの予測とどのように一致するか？ また、新しい解析は従来の矛盾を解消するか？
• RQ5 初期時系列光度測定と傾き制約を組み合わせた場合、化学組成の推定にどれほど信頼性を持って利用できるか？

主な発見

• きつい傾き事前分布を用いた最良適合モデルでは、$ar{Y}_{\rm e,dyn} = 0.22$ が得られ、傾き制約なしの従来の推定値（0.24）と比べて顕著に低減した。
• 力学的噴出物質には強い方位依存性が見られ、赤道面では $Y_{\rm e,dyn} \sim 0.17$（ランタニドが豊富、$Y_{\rm e} \lesssim 0.25$）、極部では $Y_{\rm e,dyn} \sim 0.41$（ランタニドが乏しい、$Y_{\rm e} \gtrsim 0.25$）であり、ランタニドが豊富な成分の半開口角は35°である。
• 観測者の視線はランタニドが乏しい極域の力学的噴出物質と一致しており、これによりAT2017gfoの初期青色発光が追加の駆動源なしに説明可能である。
• 推定された噴出物質質量と速度（$M_{\rm ej,dyn} = 1.3 \times 10^{-3}\,M_\odot$、$ar{v}_{\rm ej,dyn} = 0.14\,c$）は、数値相対論的予測と良好に一致しており、従来の矛盾が解消された。
• 電子分率の範囲が $Y_{\rm e,dyn} \sim 0.17$ まで下がることで、第3ピークまでにまで $r$-過程核種の強固な生成が可能となり、核合成モデルと整合的である。
• 最近の最小限の $Y_{\rm e}$ 変動を伴う球形噴出物質の主張とは矛盾し、代わりに力学的噴出物質に高密度に構造化され、非対称な組成が存在することを支持する。