Internal Heating in Magnetars: Role of Electron Captures

技术摘要：磁星内部加热：电子俘获的作用

问题陈述
磁星以极强的磁场（$B \gtrsim 10^{14}$ G）为特征，表现出持续的 X 射线光度（$\sim 10^{33} - 10^{35}$ erg/s）以及偶尔超过其自转能量输出的巨型耀发。虽然由磁应力引起的壳层形变是公认的加热机制，但它们在内壳层（此处熔化温度较高）最为有效，然而在如此深层的区域并不预期存在热源。另一种机制假定，通过磁场衰减驱动的物质压缩，磁能通过外壳中原子核的电子俘获转化为热能。此前的加热估算依赖于 HFB-24 原子质量模型；本研究通过使用更新的核数据和 HFB-27 原子质量表来研究电子俘获的作用，旨在精确预测热释放量及其来源位置。

方法论
作者将磁星的外壳建模为由原子核 $(A, Z)$ 和自由电子组成的电中性等离子体。研究假设壳层最初处于强磁场（$B_\star \equiv B/B_{rel} = 2000$，对应于 $B \approx 8.8 \times 10^{16}$ G）下的完全热力学平衡状态。关键方法步骤包括：

• 量子化效应： 对电子应用朗道-拉比（Landau-Rabi）量子化，假设由于高磁场和相对较低的温度（$T \sim 10^8 - 10^9$ K），电子占据最低朗道能级。
• 平衡组成： 通过迭代方法，利用最小化每核子的吉布斯自由能来确定初始组成。
• 反应阈值： 通过解析法计算电子俘获的起始点。考虑了磁场强度和原子核的空间排列（维格纳-塞茨近似），推导出了原子核 $(A, Z)$ 俘获电子的阈值压力 $P_\beta$。
• 热量计算： 通过考虑亲核与子核之间的能量差（包括子核态的激发能）来估算每个原子核释放的热量。研究计算了基态到基态以及基态到激发态跃迁的热释放。
• 核模型： 计算利用了 2016 年原子质量评估中的实验核质量以及来自 HFB-27 模型（基于自洽变形哈特里-福克-博戈留波夫计算）的理论预测。这些结果与之前使用 HFB-24 模型获得的结果进行了对比。
• пи克诺核聚变（Pycnonuclear Fusion）： 考虑了轻元素聚变（如碳、氧）作为最稠密区域外壳热释放的上限。

关键结果

• 主导加热反应： 研究确定了释放最多热量的特定电子俘获反应。基态跃迁包括 $^{12}\text{C} \to {}^{12}\text{Be}$ ($\sim 0.14$ MeV/nucleon)、$^{16}\text{O} \to {}^{16}\text{C}$ ($\sim 0.15$ MeV/nucleon) 和 $^{62}\text{Cr} \to {}^{62}\text{Ti}$ ($\sim 0.09$ MeV/nucleon)。基态到激发态跃迁，例如 $^{82}\text{Ge} \to {}^{82}\text{Zn}$ 以及各种锶/氪跃迁，贡献了额外的热量 ($\sim 0.07 - 0.1$ MeV/nucleon)。
• 模型对比： 虽然 HFB-27 和 HFB-24 模型预测的主要热源相似，但在次要反应方面存在显著差异。具体而言，由于计算出的 $Q_{EC}$ 值存在差异，HFB-27 模型预测的 $^{56}\text{Ti} \to {}^{56}\text{Ca}$ 和 $^{80}\text{Zn} \to {}^{80}\text{Ni}$ 跃迁的热量分别高出约 2 倍和 3 倍。
• 热量量级与位置： 每个核子的最大热释放量估计为：电子俘获约为 $\sim 0.1$ MeV，而皮克诺核聚变约为 $\sim 1 - 2$ MeV。这些反应发生在密度 $\rho_\beta \sim 10^{10} - 10^{11}$ g cm$^{-3}$ 和压力 $P_\beta \sim 10^{29} - 10^{30}$ dyn cm$^{-2}$ 的区域。
• 时间尺度与功率： 不同壳层中电子俘获的时间尺度（$\tau$）与磁星的运动学年龄（几千年）处于同一数量级。由此产生的总加热功率估计为 $W_\infty \sim 10^{35} - 10^{36}$ erg/s，与观测到的磁星光度和之前的冷却模拟一致。

意义与主张
论文得出结论，由磁场衰减驱动的外壳电子俘获和皮克诺核聚变构成了磁星内部加热的稳健来源。一个核心发现是，最大热释放量本质上是由核数据（质量和激发能）决定的，而非取决于特定的磁场强度或电子-离子等离子体的相态（固态或液态）。因此，即使对于磁场较低的中子星，这些加热估算仍然是有效的近似，尽管热源的物理位置会发生移动。作者指出，即使在壳层部分熔化的情形下，该机制仍然有效，这为依赖于壳层破裂弹性能量耗散的模型提供了明显的优势。结果与中子星冷却数据以及经验确定的热沉积密度范围保持一致。