Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics

技术摘要：具有详细原子微观物理过程的千新星喷出物中 $\beta$ 粒子输运与热化

问题陈述
千新星（Kilonova, KN）辐射由中子星合并喷出物中 $r$-过程核的放射性衰变驱动，是理解核合成、双星演化以及致密物质状态方程的关键探测手段。精确模拟千新星光变曲线需要对衰变产物（$\beta$ 粒子、$\alpha$ 粒子、$\gamma$ 射线和裂变碎片）的热化加热率进行精确计算。虽然 $\gamma$ 射线输运通常采用全辐射输运处理，但带电粒子（$\beta$ 和 $\alpha$）的热化在历史上一直依赖于简化的描述。这些简化引入了显著的不确定性，特别是在能量沉积效率、加热的空间分布以及喷出物的电离状态方面。此外，带电粒子的输运对磁场几何结构非常敏感，而该结构目前仍缺乏约束。本文旨在解决缺乏详细原子微观物理以及需要相对论性输运框架来更好地模拟膨胀千新星喷出物中 $\beta$ 粒子热化的问题。

方法论
作者结合详细的原子物理学与相对论性蒙特卡洛输运代码，采用了一种多方面的研究方法：

1. 微观物理数据与截面：

   • 核产额与能谱： 使用 Prism 反应网络计算强、$中$、弱 $r$-过程轨迹的核产额。$\beta$ 衰变率和能谱来源于 IAEA Livechart (ENSDF)，并使用 QRPA+HF 模型计算，同时辅以 BetaShape 处理近稳定核。
   • 电子-原子相互作用： 研究利用来自 EPICS2023 数据库的评估电子数据库 (EEDL) 来确定电子碰撞电离 (EII)、电子碰撞激发 (EIE)、轫致辐射和莫特散射 (Mott scattering) 的截面。
   • 电离修正： 使用经验性的 Lotz 标度法和 Hartree-Fock 结合能对中性原子截面进行电离靶标修正。
   • 次级电子： 作者计算了“电离级联”，通过追踪初级 $\beta$ 粒子产生的次级和三级电子，以确定总电离率。
   • 其他相互作用： 对 Møller 散射（电子-电子）进行了解析计算。研究发现同步辐射和切伦科夫损失相对于碰撞损失而言可以忽略不计。

2. 全局输运模型：

   • 喷出物几何结构： 将喷出物建模为球对称、同比例膨胀的流体，具有幂律密度分布（$\rho \propto r^{-\alpha}$）。
   • 输运情景： 考虑了两种极限磁场几何结构：
     • 俘获极限 (Trapped Limit)： 对应于高度缠绕或环向磁场，粒子被限制在特定的速度坐标内。
     • 径向场极限 (Radial Field Limit)： 对应于相干径向磁场，粒子沿磁力线自由流经。作者证明，对于相干径向场，1D 输运方程在数学上等价于无磁化自由流情况，因为磁镜力与几何效应相互抵消。
   • 蒙特卡洛模拟： 一个全相对论性蒙特卡洛代码用于追踪离散的 $\beta$ 粒子包。它考虑了由于膨胀导致的绝热冷却、碰撞能量损失（阻止本领）以及散射（输运不透明度）。粒子在热化（$E < 1$ eV）、逃逸（$u > v_{ej}$）或达到最大时间（1 年）时终止。

3. 电离结构：

   • 应用简化的速率方程模型于钕 (Nd) 元素，以说明输运和电离级联对喷出物温度和电离状态的影响，同时忽略了光电离以及不同元素之间的耦合。

主要贡献与结果

• 原子微观物理： 作者提供了已评估的、按物种分类的阻止与输运截面，以及 EII 截面，并可在网上获取。研究表明，在 MeV 能级下，EII 主导能量损失；而在较低能量（热化阶段）时，与自由热电子的库仑相互作用（Møller 散射）变得占据主导地位。
• 电离级联： 研究量化了次级电子对总电离率的贡献。级联过程显著增强了电离截面，在 1 MeV 时其贡献可与初级电子相当。
• 热化效率：
  • 非局域效应： 非局域能量沉积和粒子逃逸降低了热化效率，特别是在最内层和最外层喷出物中。这使得热化效率在喷出物中分布得比局部沉积模型更均匀。
  • 输运极限： 在径向场极限下，粒子逃逸是导致低效的主要驱动因素。在俘获极限下，由于粒子在喷出物中停留时间较长，绝热膨胀损失是导致低效的原因。
  • 参数依赖性： 热化效率对喷出物质量、速度、$r$-过程丰度模式及密度剖面斜率较为敏感。它对电离状态的敏感度较低，这表明在某些语境下，无需自洽的电离计算即可对能量沉积进行建模。
• 解析描述： 作者推导出了一个空间依赖的热化效率解析公式 $f_{thm}(t) = \ln(1+\tau)/\tau$，该公式通过一个不效率时间 ($t_{ineff}$) 和一个渐近斜率 ($n$) 进行参数化。他们提供了这些参数随喷出物质量、速度和位置变化的拟合参数。蒙特卡洛结果与该描述之间的中值残差约为 $\sim 0.04$ dex（径向）和 $\sim 0.02$ dex（俘获）。
• 电离结构： 对 Nd 的演示性计算揭示了星云期存在“反转”电离结构（外层密度较低区域的电离度更高），这是由于复合效率较低所致。径向磁场会导致较低的电离态（由于粒子逃逸），而包含电离级联则会提高整个喷出物的电离态。

意义与主张
本文断言，下一代千新星模型必须纳入评估的原子数据和显式的带电粒子输运，以减少光变曲线预测的不确定性。作者主张：

1. 输运物理至关重要： 忽略非局域输运效应会导致高估内层和外层喷出物的加热量，并导致错误的电离结构。
2. 原子数据至关重要： 使用评估的原子数据 (EEDL) 并包含电离级联，对于获得准确的电离率（进而影响不透明度和冷却）是必要的。
3. 解析实用性： 提供的解析描述提供了一种计算成本极低的、空间依赖的热化效率近似方法，适用于半解析和辐射传输模型，尽管在高保真模拟中仍建议使用显式输运。
4. 局限性： 作者谦虚地指出，其结果是基于理想化的球对称和极限磁场几何结构。他们建议，现实的磁场形态、非球形喷出物以及集体等离子体效应（例如流体不稳定性）是未来研究的方向。他们还强调，其结果侧重于 $\beta$ 粒子，应与未来关于 $\alpha$ 粒子和裂变碎片的 $\beta$ 粒子研究相结合，以构建完整的热化图景。