On the impact of the carbon fusion rate over the properties of superbursts -- Numerical simulations of superbursts with MESA

技术摘要：关于碳聚变速率对超爆（superbursts）性质影响的研究

问题陈述
超爆是发生在低质量X射线双星（LMXBs）中中子星壳层内的剧烈热核爆炸，由温度 $T \le 10^9$ K 下的不稳定碳燃烧触发。虽然之前的多区模拟（例如 Keek & Heger 2011, 2012）已成功模拟了超爆的复发时间和光变曲线，但这些研究依赖于经典的碳聚变速率（Caughlan & Faber 1988, CF88）。近期的核物理研究对该速率提出了挑战，表明在相关温度（$T \le 10^9$ K）下，实际速率可能显著不同。具体而言，来自 LUNA 实验的数据（Tumino et al. 2018）表明，速率可能会提高约 $\sim 10^3$ 倍；而阻碍模型（Monpribat et al. 2022）则表明速率可能会降低约 $\sim 10^3$ 倍。这些巨大的变化对超爆可观测物理量（如复发时间、衰减时间、点火深度和核合成产物丰度）的影响，在现代恒星演化代码的背景下仍有待探索。

研究方法
作者利用公开的恒星演化代码 MESA v.24.08.1 模拟了经历质量吸积过程的中子星包层演化。模拟使用了缩减的核反应网络 Approx_A56（覆盖核素至 $A=56$），并使用 Approx_A64 进行部分运行以验证截断效应。

关键参数与配置包括：

• 吸积率： 测试了两种典型速率：一种是“中等”速率（$\dot{M}_{moderate} \approx 5.88 \times 10^{-9} M_\odot \text{yr}^{-1}$）和一种是“高”速率（$\dot{M}_{high} \approx 2 \times 10^{-8} M_\odot \text{yr}^{-1}$）。
• 组分： 吸积物质被建模为 20% 的 $^{12}\text{C}$ 和 80% 的 $^{56}\text{Fe}$，这符合触发超爆的文献要求。初始包层为纯 $^{56}\text{Fe}$。
• 底部加热： 在包层底部施加了一个固定的底部光度（$L_b$）以确保不稳定的碳点火，该值经过校准以匹配之前的 KEPLER 模拟。
• 碳聚变速率： 对四种不同的反应速率进行了比较：
  1. CF88： 经典速率（基准模型）。
  2. LUNA： 基于 LUNA 测量，预测在 $T \sim 2\text{--}3 \times 10^8$ K 时，其速率比 CF88 高约 $\sim 10^3$ 倍。
  3. HIN： 基于阻碍模型，预测其速率比 CF88 低约 $\sim 10^3$ 倍。
  4. HIN-RES： 包含低能级共振态的 HIN 模型，其结果介于 HIN 和 CF88 之间。
• 物理处理： 模拟包含了辐射不透明度（Wolf 方案）、对流的混合长度理论以及固定的表面压力边界条件。省略了广义相对论修正。

主要结果
本研究量化了碳聚变速率的变化如何改变超爆的性质：

1. 时间尺度与点火深度：

   • 增强速率（LUNA）： 导致更短的复发时间和衰减时间，并在更浅的柱密度（$y_{ignition}$）处触发点火。对于 $\dot{M}_{moderate}$，复发时间从 $\sim 241$ 天（CF88）降至 $\sim 100$ 天。
   • 降低速率（HIN）： 导致更长的复发和衰减时间，且点火发生在更深的柱密度处。对于 $\dot{M}_{moderate}$，复发时间增加至 $\sim 364$ 天。
   • 中间速率（HIN-RES）： 产生的结果与 CF88 基准模型几乎无法区分。

2. 温度与光变曲线：

   • 爆炸期间达到的最高温度对速率具有敏感性。HIN 速率产生最高的峰值温度（$\sim 5 \times 10^9$ K），而 LUNA 产生的峰值温度较低（$\sim 3 \times 10^9$ K）。
   • 光度衰减剖面随之移动：LUNA 模拟显示快速衰减，而 HIN 模拟显示缓慢衰减。
   • 在 $\dot{M}_{high}$ 时，峰值光度对速率的选择基本不敏感，但衰减时间和复发间隔仍有显著差异。

3. 核合成与灰烬组分：

   • 所有模型的灰烬分布均表现出特征性的双峰结构（分别围绕 $^{28}\text{Si}$ 和 $^{56}\text{Fe}$）。
   • 速率敏感性： 对于 $A \le 44$，$\alpha$-核素的丰度对速率高度敏感。LUNA 速率（更高温/更高效率）产生的 $\alpha$-核素显著多于 CF88/HIN-RES，而 HIN 速率产生的较少。
   • 高吸积率效应： 在 $\dot{M}_{high}$ 时，在 $28 \le A \le 50$ 的核素合成方面有显著增强，并且出现了 $\sim 0.1\%$ $^{16}\text{O}$ 的特定积累，这一特征在中小吸积率下并不显著。
   • 光解作用： 通过将惰性 $^{56}\text{Fe}$ 替换为 $^{72}\text{Ge}$ 的模拟证实，铁峰灰烬的分布主要由 $^{56}\text{Fe}$ 的光解作用驱动，而非仅仅由碳聚变驱动。移除光解通道会阻止超过 $^{40}\text{Ca}$ 的核合成。

4. 与底部加热的简并性：

   • 研究表明，改变碳聚变速率对超爆时间尺度和点火深度的影响，与改变包层的底部光度（$L_b$）相当。例如，降低速率（HIN）并增加 $L_b$ 会产生相似的复发时间。这表明存在一种简并性，即仅通过拟合光变曲线无法唯一确定融合速率和底部加热。

意义与结论
本文得出结论，在 $T \le 10^9$ K 下的碳聚变速率是决定超爆物理性质的关键参数。作者认为：

• 碳聚变速率的不确定性可以解释观测到的复发和衰减时间的差异，而不必必然引入质量吸积率或底部加热的变化。
• 基于之前使用 CF88 的 KEPLER 模拟所得出的“典型”灰烬组分，在双峰结构方面是稳健的，但在 $\alpha$-核素丰度方面对特定的反应速率具有敏感性。
• 碳聚变速率与底部加热之间的简并性意味着，仅通过拟合单个超爆光变曲线，无法明确判定底层的核物理过程或中子星壳层的热状态。
• 本研究验证了使用 MESA 进行超爆建模的有效性，显示其与之前的 KEPLER 结果保持一致，同时提供了一个用于探索核物理不确定性的公开代码替代方案。

作者谦虚地指出，虽然他们的发现表明，如果 LUNA 速率是正确的，那么碳聚变速率的变化可以减轻模型对高底部加热的需求，但由于特定源（如 4U 1636-53）质量吸积率的不确定性，目前尚无法得出关于哪种速率在物理上实现的定论。此外，他们明确指出，由于密度和温度处于热核反应范畴，将其处理为纯热核过程是合理的，这与压致核反应（pycnonuclear）机制不同。