Radiation-hydrodynamics of star-disc collisions: From system parameters to outflows and lightcurves

技术摘要：恒星-吸积盘碰撞的辐射流体动力学：从系统参数到流出物与光变曲线

问题与背景
准周期爆发（QPEs）是与超大质量黑洞（SMBHs）相关的核瞬变现象，其特征是重复出现的软X射线耀发。虽然存在多种模型，但“恒星-吸积盘碰撞”情景——即一颗处于紧凑、偏心且倾斜轨道的恒星反复穿过吸积盘——是一个领先的候选模型。该模型自然地解释了如 GSN 069 等源中所观察到的复现时间、交替的“长-短”间隔以及“强-弱”耀发模式。然而，要将这些碰撞的局部物理参数与可观测的耀发特性（亮度、持续时间、不对称性）联系起来，需要理解在膨胀且光学厚的介质中，激波、流出物几何形状与辐射传输之间的耦合演化。以往的研究利用 1D、2D 或非辐射方法探索了该问题的部分内容，但一直缺乏一项专门的、系统的、将系统参数与光变曲线联系起来的 3D 辐射流体动力学研究。

方法论
作者使用平滑粒子流体动力学（SPH）代码 Phantom 进行了一系列三维局部辐射流体动力学模拟。模拟是在以吸积盘中平面为中心的局部笛卡尔坐标系中进行的，忽略了全局轨道动力学、广义相对论以及气体的自引力，因为碰撞时间尺度远短于轨道时间尺度。

研究采用了一个基准模型，该模型描述一颗主序星（$R_\star \approx R_\odot$）与围绕 $10^6 M_\odot$ 超大质量黑洞的、光学厚且由辐射压主导的 $\alpha$-吸积盘发生碰撞。恒星被建模为一个以 $v_\star = 0.1c$ 运动的刚性球体。作者系统地改变了相对于基准模型的五个关键参数，每次仅改变一个变量：

1. 恒星速度 ($v_\star$)： 在基准值的 $0.5$到$1.5$ 倍之间变化。
2. 吸积盘表面密度 ($\Sigma_d$)： 在基准值的 $1$到$100$ 倍之间变化，涵盖了从标准薄盘到致密的潮汐瓦解事件（TDE）后吸积盘。
3. 恒星半径 ($R_\star$)： 在基准值的 $0.5$到$4$倍之间变化，有效地改变了盘标高与恒星半径的比值 ($H/R_\star$)。
4. 垂直吸积盘密度分布： 建模为高斯分布 $\rho(z) \propto e^{-z^2/2\sigma^2}$，改变向中平面集中的质量浓度 ($\sigma/R_\odot$)。
5. 局部碰撞角 ($i$)： 相对速度矢量与吸积盘中平面之间的夹角，范围从 $90^\circ$（垂直）到 $30^\circ$（斜向）。

辐射传输采用基于电子散射不透明度的通量限制扩散近似处理。通过动态确定的光球层积分扩散通量和平流通量，计算得出涌现的玻色洛度（bolometric luminosity）。

主要贡献与结果
本研究提供了峰值亮度 ($L_{\text{peak}}$) 和耀发持续时间 ($\Delta t$) 作为系统参数函数的经验定标关系。主要发现如下：

• 碰撞动力学： 恒星驱动一个弓形激波，在吸积盘两侧产生两个流出物（前向和后向）。前向流出物通常质量更大且能量更高。尽管在 $v_\star$ 和 $\Sigma_d$ 的变化下，弓形激波和流出物的定性形态基本保持不变，但热力学归一化发生了变化。
• 恒星速度 ($v_\star$)： 较快的恒星产生更亮的耀发（$L_{\text{peak}} \propto v_\star^2$），这是由于更高的激波加热率。它们也产生更长的耀发（$\Delta t \propto v_\star^{1.4-2.2}$），因为注入了更多能量，尽管由于更强的绝热损失，其亮度衰减更快。
• 吸积盘表面密度 ($\Sigma_d$)： 更致密的盘增加了流出物的质量和动量，但对峰值亮度的影响较弱（前向流出物 $L_{\text{peak}} \propto \Sigma_d^{-0.1}$，后向流出物 $L_{\text{peak}} \propto \Sigma_d^{0.2}$）。这是由于增加的激波加热与光学厚度更高的喷出物中更强的绝热损失之间的竞争。然而，致密的盘显著延长了耀发持续时间（$\Delta t \propto \Sigma_d^{0.5-0.6}$），这是由于更长的光子扩散时间。
• 恒星半径 ($R_\star$)： 更大的恒星与更多气体发生作用，产生更亮（前向 $L_{\text{peak}} \propto R_\star^{1.5}$，后向 $\propto R_\star^{0.7}$）且更长（$\Delta t \propto R_\star^{1.5}$）的耀发。增加 $R_\star$ 会增强前向与后向流出物之间的动量不对称性。
• 垂直密度分布 ($\sigma$)： 更中心集中的盘（较小的 $\sigma$）产生更亮且更短的耀发。致密的中心层吸收了更多激波能量，而低密度的外层则允许辐射以更低的绝热损失更有效地逃逸。
• 碰撞角 ($i$)： 更斜的穿过（较小的 $i$）会降低两个流出物之间的动量和亮度不对称性。后向流出物变得显著更亮（$L_{\text{peak}} \propto (\sin i)^{-0.8}$），且耀发变得更长（$\Delta t \propto (\sin i)^{-0.9 \text{ 至 } -1}$），这是由于更长的相互作用时间和更对称的受激气体重新分布。

对 GSN 069 的应用
作者将这些经验定标应用于特征明确的源 GSN 069。通过采样参数空间以匹配观测到的峰值亮度（$\sim 3 \times 10^{42}$ erg s$^{-1}$）、持续时间（$\sim 1$ 小时）以及强-弱振幅比，他们确定了两个可行的区域：

1. 首选解： 一个逆行轨道，碰撞于一个致密且垂直集中的 TDE 后吸积盘（$\Sigma_d \sim 45 \times$ 基准值，$\sigma/R_\odot \sim 0.3$）。该模型很好地重现了观测到的耀发特性。
2. 替代解： 一个具有标准薄盘密度但具有极度掠过式碰撞几何形状（$i \sim 6^\circ$）的逆行轨道。该解准确性较低，且需要极端参数。

作者指出，对致密且垂直集中的盘的偏好，与 GSN 069 的 TDE 类活动观测证据是一致的。

意义与局限性
本文声称提供了第一个将局部碰撞参数与 QPE 光变曲线联系起来的系统性 3D 辐射流体动力学研究。所推导的经验定标关系作为局部模拟与全局 QPE 模型之间的桥梁，允许根据系统参数预测耀发特性。作者强调，其发现表明特定的参数组合可以重现 QPE 特有的振幅、持续时间和占空比。

研究也适度地承认了局限性：模型假设参数定标是可分离的，忽略了盘的旋转与剪切，将恒星视为刚性球体，并且未考虑观测视角效应或光谱演化（例如康普顿化）。因此，结果是以物理动机驱动的参数空间区域呈现，而非唯一的、决定性的最佳拟合模型。作者建议，未来的工作应纳入真实的恒星结构、吸积盘动力学以及依赖于观测视角的亮度，以全面量化该模型重现多样化观测 QBE 源的能力。