[论文解读] 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的千新星AT2017gfo进行建模，揭示观测者视线与富含轻稀土元素但贫镧系元素的极向动力学喷出物成分相交，其平均电子分数为$\bar{Y}_{\rm e,dyn} = 0.22$。结果证实喷出物成分具有显著的方位角依赖性，富含镧系元素的物质被限制在赤道平面周围35°的半开角范围内，解决了与数值相对论模拟之间的先前矛盾，并支持从第三峰开始的稳定$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中喷出物特性的推断。
• 通过采用多维、依赖倾斜角的喷出物建模方法，解决先前一维千新星模型与数值相对论模拟之间的差异。
• 确定AT2017gfo中观测到的早期蓝光是否可由非球对称、成分分层的动力学喷出物解释，而无需引入额外能量源。
• 基于动力学喷出物中推断的电子分数分布，评估从第三峰开始的$r$-过程核合成的稳健性。
• 证明结合早期光变曲线与倾斜角约束，可实现对千新星中化学成分与喷出物几何结构的可靠推断。

提出的方法

• 采用基于网格的多维千新星建模框架，整合最先进的加热速率、能量沉积效率及频率依赖性不透明度。
• 通过$\bar{Y}_{\rm e,dyn}$（平均电子分数）参数化动力学喷出物，并基于数值相对论模拟假设其具有方位角依赖性：$Y_{\rm e,dyn}$从赤道处的~0.17变化至极地处的~0.41。
• 利用在合成光曲线数据上训练的前馈神经网络，实现在全参数空间内的快速贝叶斯参数估计。
• 对视线方向施加严格先验（$19^\circ < \theta_{\rm obs} < 25^\circ$），该先验来自AT2017gfo中相对论喷流的自行测量。
• 将模型拟合至AT2017gfo的多波段光度数据，以推断喷出物质量、速度及电子分数，同时对动力学喷出物与盘风喷出物成分施加约束。
• 通过与数值相对论预测对比，验证结果，并与近期基于Sr II谱线特征提出球对称喷出物的主张进行比较。

实验结果

研究问题

• RQ1 在GW170817的千新星AT2017gfo中，引入精确倾斜角约束是否显著改变了对动力学喷出物特性的推断？
• RQ2 AT2017gfo中的早期蓝光是否可由非球对称、成分分层的动力学喷出物解释，而无需引入额外能量源？
• RQ3 动力学喷出物中电子分数的方位角依赖性如何？其对$r$-过程元素生成有何影响？
• RQ4 推断的喷出物参数与数值相对论模拟预测相比如何？新分析是否解决了先前的矛盾？
• RQ5 在仅结合早期光变曲线与倾斜角约束的前提下，能多大程度上可靠推断千新星喷出物的化学成分？

主要发现

• 在采用严格倾斜角先验的最优拟合模型中，得到$\bar{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^\circ$。
• 观测者视线与贫镧系元素的极向动力学喷出物对齐，无需额外能量源即可解释AT2017gfo中的早期蓝光。
• 推断的喷出物质量与速度——$M_{\rm ej,dyn} = 1.3 \times 10^{-3}\,M_\odot$，$\bar{v}_{\rm ej,dyn} = 0.14\,c$——与数值相对论预测高度一致，解决了先前的矛盾。
• 电子分数范围可低至$Y_{\rm e,dyn} \sim 0.17$，足以实现从第三峰开始的稳定$r$-过程核素生成，与核合成模型一致。
• 结果反驳了近期关于球对称喷出物且$Y_{\rm e}$变化极小的主张，反而支持动力学喷出物中存在高度结构化、非对称的成分分布。