[论文解读] Crustal heating in accreting neutron stars from the nuclear energy-density functional theory. I. Proton shell effects and neutron-matter constraint

A. F. Fantina, J. L. Zdunik+4 more|Lirias (KU Leuven)|Jun 11, 2018

Pulsars and Gravitational Waves Research参考文献 85被引用 3

一句话总结

本研究采用核能量密度泛函（EDF）理论模拟吸积中子星的壳层加热，明确引入质子壳效应与核物质状态方程（EoS）约束。结果表明，壳层闭合——尤其是Z=14处——主导了内壳层的热量沉积，导致总壳层加热达1.5–1.7 MeV，显著高于以往模型，归因于更真实的核结构效应。

ABSTRACT

Observations of soft X-ray transients in quiescence suggest the existence of heat sources in the crust of accreting neutron stars. Heat is thought to be released by electroweak and nuclear processes triggered by the burying of ashes of X-ray bursts. The heating is studied using a fully quantum approach taking consistently into account nuclear shell effects. We have followed the evolution of ashes made of $^{56}$Fe employing the nuclear energy-density functional theory. Both the outer and inner crusts are described using the same functional, thus ensuring a unified and thermodynamically consistent treatment. To assess the role of the neutron-matter constraint, we have employed the set of accurately calibrated Brussels-Montreal functionals BSk19, BSk20, and BSk21 and for comparison the SLy4 functional. Due to nuclear shell effects, the fully accreted crust is found to be much less stratified than in previous studies. In particular, large regions of the inner crust contain clusters with the magic number $Z=14$. The heat deposited in the outer crust is tightly constrained by experimental atomic mass data. The shallow heating we obtain does not exceed $0.2$~MeV and is therefore not enough to explain the cooling of some soft X-ray transients. The total heat released in the crust is very sensitive to details of the nuclear structure and is predicted to lie in the range from $1.5$~MeV to $1.7$~MeV. The evolution of an accreted matter element and therefore the location of heat sources are governed to a large extent by the existence of nuclear shell closures. Ignoring these effects in the inner crust, the total heat falls to $\sim 0.6$~MeV. The neutron-matter constraint is also found to play a key role. The large amount of heat obtained by Steiner et al. (2012) could thus be traced back to unrealistic neutron-matter equations of state.

研究动机与目标

• 理解吸积中子星在宁静期壳层加热的起源与分布。
• 评估核壳效应——尤其是质子壳层闭合——对内壳层热量沉积的影响。
• 评估中子物质状态方程（EoS）约束对总壳层加热的影响。
• 调和早期模型（如HZ、Steiner 2012）与软X射线暂现源观测约束之间的差异。
• 使用同一EDF提供外层与内层壳层的热力学一致、统一的描述。

提出的方法

• 采用核能量密度泛函（EDF）理论，使用布鲁塞尔-蒙特利尔泛函（BSk19、BSk20、BSk21）并以SLy4作为对比。
• 通过自洽EDF计算追踪56Fe残余物在壳层各层的演化，包含配对与介质效应。
• 沿吸积路径计算电子捕获、中子发射/吸收及聚变核反应的热量释放。
• 利用实验原子质量数据约束外层壳层中的反应，确保与观测核质量一致。
• 应用Chamel & Fantina（2016a）的解析公式计算中子滴落密度与压强，以定义壳层边界。
• 与Mackie & Baym（1977）的可压缩液滴模型（CLDM）比较，以分离壳效应的影响。

实验结果

研究问题

• RQ1 质子壳效应，特别是Z=14处，如何影响内壳层中壳层加热的分布与大小？
• RQ2 中子物质状态方程在多大程度上影响壳层中总热量的沉积？
• RQ3 为何早期模型（如HZ、Steiner 2012）高估了壳层加热，而与观测约束不符？
• RQ4 在EDF理论中引入核壳效应后，与冷催化壳层相比，吸积壳层的分层结构与组成有何变化？
• RQ5 实验原子质量数据在多大程度上约束了外层壳层中的热量沉积？

主要发现

• 核壳效应，尤其是幻数Z=14，导致内壳层中形成大范围、扩展的簇集区域，相比以往模型降低了分层程度。
• EDF计算将总壳层加热严格约束在1.5–1.7 MeV之间，BSk19–BSk21与SLy4结果一致。
• 外层壳层加热受限于≤0.2 MeV，不足以解释某些软X射线暂现源中的浅层加热现象。
• CLDM（MB）模型中缺乏壳效应，使总加热降低至约0.6 MeV，凸显了壳层闭合的关键作用。
• Steiner等人（2012）报告的高加热值源于不现实的中子物质EoS与经验壳参数化，而非物理现实。
• EDF方法得到的壳层加热在密度约1.6×10¹² g cm⁻³时与Lau等人（2018）结果良好一致，但在更高密度下因一致处理介质修正壳效应而出现偏离。