Hydrodynamic simulations of the recurrent nova T Coronae Borealis: Nucleosynthesis predictions

技术摘要：T Coronae Borealis 的流体动力学模拟：核合成预测

问题与背景
复发新星（RNe）是典型的双星系统，其特征是在白矮星（WD）表面发生的热核失控（TNR），其复发周期极短，通常在 1 到 100 年之间。这种短间隔要求极高的质量吸积率、接近钱德拉塞卡极限的白矮星质量以及较高的初始白矮星光度。T Coronae Borealis（T CrB）是一个著名的银河系复发新星，曾于 1866 年和 1946 年发生爆发，预计其下一个复发周期约为 80 年。鉴于其下一次爆发迫在眉睫，迫切需要约束该系统的物理参数——特别是白矮星质量、初始光度、吸积率和金属丰度——以理解驱动此类爆发的条件。此外，复发新星是 Ia 型超新星的潜在前身，因此表征其核合成过程和质量积累对于理解其演化至关重要。

方法论
作者利用流体动力学、拉格朗日、有限差分代码 SHIVA 对 T CrB 的爆发进行了模拟。模拟采用了包含简并电子气体、多组分离子等离子体和辐射的广义状态方程，并使用了 OPAL 罗斯兰平均不透明度。核能量生成通过包含 120 种物种（从 $^1$H 到 $^{48}$Ti）并由 630 个过程链接的反应网络进行计算，该网络整合了来自 STARLIB 数据库的更新速率。

本研究提出了 11 个新的流体动力学模型，旨在通过改变不同条件来重现 80 年的复发周期（$\tau_{rec}$）：

• 白矮星质量 ($M_{WD}$): 1.30, 1.35, 和 1.38 $M_{\odot}$。
• 初始白矮星光度 ($L_{WD}$): 0.01, 0.1, 和 1.0 $L_{\odot}$。
• 吸积金属丰度 ($Z_{acc}$): 太阳组成 ($1.0 Z_{\odot}$)、次太阳组成 ($0.1 Z_{\odot}$) 以及超太阳组成 ($10 Z_{\odot}$)。
• 混合作用: 为 1.35 和 1.38 $M_{\odot}$ 的白矮星计算了另外 8 个模型，以评估吸积包层与下方白矮星层（碳-氧及氧-氖组成）之间混合对结果的影响。
• 低质量白矮星: 计算了 3 个针对 1.20 $M_{\odot}$ 白矮星的模型，以涵盖质量估计较低的复发新星情况。

对于每种配置，质量吸积率 ($\dot{M}_{acc}$) 都经过调整，以确保在经过恰好 80 年的吸积后触发爆发。

关键结果

• 吸积率与物理条件： 若要使用太阳组成物质实现 80 年的复发周期，对于质量为 $1.30 - 1.38 M_{\odot}$ 且光度为 $0.01 - 1.0 L_{\odot}$ 的白矮星，所需的质量吸积率 $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$。较低的初始白矮星光度、较低的白矮星质量或降低吸积物质的金属丰度，都需要更高的吸积率才能在 80 年的时间窗口内达到临界点火压力。
• 爆发动力学： 具有较低 $L_{WD}$ 或较低 $Z_{acc}$ 的模型会导致更高的累积包层质量和更高的点火压力。因此，这些模型会产生更剧烈的爆发，其特征是更高的峰值温度 ($T_{max}$)、更大的抛射质量 ($\Delta M_{eje}$) 和更大的动能。例如，一个 1.38 $M_{\odot}$ 的白矮星在 $L_{WD}=0.01 L_{\odot}$ 且具有太阳金属丰度时，其 $T_{max} \approx 2.92 \times 10^8$ K；而一个 1.30 $M_{\odot}$ 的白矮星在 $L_{WD}=1.0 L_{\odot}$ 时，仅达到 $2.01 \times 10^8$ K。
• 核合成与元素丰度：
  • 同位素比率： 抛射物由 CNO 循环同位素（$^{12,13}$C, $^{14,15}$N）主导。在不同模型之间观察到同位素比率（如 $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si）存在显著差异。
  • 白矮星质量依赖性： 具有相同初始光度但质量不同的白矮星模型，在 Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, 和 Sc 等元素的丰度上表现出显著差异。具体而言，1.38 $M_{\odot}$ 模型相对于太阳丰度在 S, Cl, K, Ca, 和 Sc 方面表现出显著富集，而 1.20 $M_{\odot}$ 和 1.30 $M_{\odot}$ 模型通常显示出太阳或略低于太阳的丰度（除模型 130C 外）。硅（Silicon）在 1.20 和 1.30 $M_{\odot}$ 模型中得到增强，但在 1.38 $M_{\odot}$ 模型中则呈现亏损或接近太阳水平。
  • 金属丰度效应： 在 1.38 $M_{\odot}$ 模型中，改变金属丰度会剧烈改变组成。低金属丰度 ($0.1 Z_{\odot}$) 会导致显著的 $^{40}$Ca 生成，而高金属丰度 ($10 Z_{\odot}$) 则使主导成分转向 N, C, Ne, 和 Si。
• 混合情景： 与非混合模型相比，纳入了包层与白矮星基底（CO 或 ONe）之间混合的模型缩短了吸积阶段并减少了总抛射质量。至关重要的是，混合显著增强了伽马射线发射体（$^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al）的产生。虽然非混合模型预测的伽马射线发射微乎其微，但混合模型可能会产生可探测的特征，从而可能区分 CO 型和 ONe 型白矮星。
• 伽马射线发射： 对于标准的非混合模型，由于放射性同位素产量较低，作者预测伽马射线谱线（511 keV 和 1275 keV）及其相关的连续谱将比当前仪器的探测极限低几个数量级。

意义与主张
本文断言，抛射物的成分差异，特别是诸如 Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, 和 Sc 等物种的元素丰度，可作为一种潜在的诊断工具。这些变化可以约束参数空间，并区分本研究提出的各种 T CrB 模型。具体而言，作者声称，对即将到来的 T CrB 爆发进行的早期光谱测量可以确定白矮星质量和混合程度，从而阐明控制该系统爆发的物理条件。研究结论指出，尽管在没有显著混合的情况下检测到伽马射线是不太可能的，但元素丰度模式为识别底层白矮星性质及所发生的混合过程提供了一种稳健的方法。