Anisotropic Ejecta from Binary Neutron Star Mergers: Self-Consistent Main Thermal and Late-Time Radio Emission of NS-Powered Kilonovae

技术摘要：双中子星并合的各向异性抛射物

问题陈述
双中子星（NS）并合后遗迹的性质——具体是形成大质量中子星还是坍缩成黑洞——仍是天体物理学中的一个根本问题。虽然大质量中子星中心引擎可以通过自转减慢能量为千新星的主要热辐射和长期射电信号提供动力，但目前的射电观测未能探测到如 GW 170817 等事件所预测的晚期射电信号。先前的理论估算通常假设各向同性的抛射物分布，这可能高估了射电流量。此外，AT 2017gfo 观测到的双分量（“蓝”和“红”）热辐射表明存在复杂的各向异性抛射物几何结构。本研究旨在填补理解空白，探讨各向异性抛射物分布与中心磁星引擎相结合如何影响主要热辐射和晚期射电辐射的演化，以及该框架是否能解释 GW 170817 中射电信号的未探测现象。

方法论
作者开发了一个自洽的、角度依赖的模型，用于描述由并合后毫秒磁星驱动的千新星动力学。

• 抛射物分布： 将并合抛射物建模为双分量系统：赤道动力学抛射物（高密度、富含镧系元素、“红”分量）和极区风（低密度、无镧系元素、“蓝”分量）。质量分布使用由过渡角 $\theta_t$ 和密度比 $f_\eta$ 定义的角度依赖函数 $\eta(\theta)$ 进行参数化。
• 能量注入： 中心引擎被视为自转减慢的磁星。研究纳入了各向异性能量注入，其中自转减慢光度随纬度变化，其特征由极区与赤道能量注入之间的比率 $f_\xi$ 表征。
• 动力学演化： 作者针对每个立体角求解相对论流体力学方程，考虑了能量注入、放射性衰变、内能扩散以及与星际介质（ISM）的相互作用。考虑了两种 ISM 相互作用情景：完全扫掠（标准情形）和不完全扫掠（考虑中性氢抑制）。
• 辐射计算： 该模型通过沿观测者视线积分角度依赖的光度，并应用多普勒增亮和朗伯余弦定律，计算主要热变光曲线。晚期射电辐射源于加速抛射物与星际介质之间激波相互作用产生的同步辐射。
• 应用于 GW 170817： 该模型被应用于 AT 2017gfo 事件。作者拟合多波段热变光曲线以约束抛射物质量、速度和引擎参数，然后利用这些拟合参数预测相应的晚期射电变光曲线。

主要贡献与结果

• 双峰特征： 研究表明，各向异性抛射物分布自然地在主要热辐射和晚期射电变光曲线中产生双峰结构。这两类辐射的时间演化内在关联，因为它们都敏感地依赖于特定的抛射物几何结构。
• 热辐射与射电辐射的一致性： 在中子星驱动的情景下，主要热辐射和晚期射电辐射在变光曲线演化上表现出高度的一致性。这种相关性比放射性衰变驱动的情景更强，因为在后者中能量注入与抛射物分布耦合，而非作为外部驱动源。
• 各向异性参数的影响：
  • 赤道区域与极区之间需要显著的密度对比（$f_\eta \approx 100$），才能产生明显的双峰热轮廓。
  • 各向异性能量注入（$f_\xi > 1$）通过向极区方向注入更多能量，增强了早期蓝分量（极区风），导致蓝分量更早到达并具有更高的光度。
• GW 170817 分析：
  • 拟合 AT 2017gfo 的热数据得出抛射物质量约为 $0.0085 M_\odot$，中心中子星初始周期为 $6.6$ 毫秒，以及各向异性参数，表明存在更致密的赤道盘和更强的极区能量注入。
  • 推导出的参数意味着中心中子星的旋转能量显著减少（可能通过引力波辐射），仅剩约 $10^{51}$ 尔格用于加速抛射物。
  • 因此，GW 170817 预测的晚期射电变光曲线低于当前的观测探测极限，特别是对于星际介质密度 $n < 10^{-1} \text{ cm}^{-3}$ 的情况。这一结果与现有观测中未探测到射电信号的现象相符。

意义
本文建立了一个自洽的演化图景，通过抛射物的几何结构将千新星的主要热辐射与晚期射电辐射联系起来。其主要意义在于证明，GW 170817 晚期射电信号的未探测并不一定与磁星中心引擎模型相矛盾；相反，当考虑各向异性抛射物分布和能量损失（例如通过引力波）时，这与理论预期是一致的。

该研究提出，热辐射与射电辐射的联合分析可作为强有力的工具，用于：

1. 约束并合抛射物的几何结构。
2. 探测中心引擎的性质（例如自转周期和能量注入的各向异性）。
3. 减少多信使观测解释中的模型简并性。

最终，这项工作表明，热变光曲线与射电变光曲线之间的联系为磁星模型提供了双重一致性检验，提供了一条更好地理解并合后遗迹以及新生中子星旋转能量命运的路径。