Bayesian Inference of Stellar r-Process Abundances

技术摘要：恒星 r-过程丰度的贝叶斯推断

问题陈述
快速中子捕获过程（r-过程）是宇宙中重元素合成的主要机制。虽然在贫金属、r-过程增强星以及太阳 r-过程残余物中，重元素（A > 130）的丰度模式表现出显著的稳健性，但在较轻重元素（Sr, Y, Zr，直至第二峰）方面存在显著的星际差异。以往解释这些模式的研究提出了多种机制，包括原生过程、裂变、弱 r-过程或多个天体物理场所的贡献。然而，确定能够重现观测到的恒星丰度的具体天体物理条件（电子分数、熵、膨胀时间尺度）仍然具有挑战性。现有的研究通常依赖于拟合固定数量的组分或特定的轨迹，而缺乏一个严谨的统计框架来量化不确定性或确定所需贡献条件的必要数量。

方法论
作者提出了一个利用马尔可夫链蒙特卡洛（MCMC）采样进行贝叶斯推断的框架，用以推断 r-过程的天体物理条件。该方法集成了以下组成部分：

• 核合成网格： 分析依赖于一个大规模、与场所无关的调查（KAR25），其中包含 120,000 个参数化核合成计算。这些计算使用 WinNet 反应网络和 FRDM2012 质量模型，并由初始电子分数 ($Y_e$)、初始熵 ($s$) 和膨胀时间尺度 ($\tau$) 进行参数化。
• 统计模型： 观测到的恒星丰度模式被建模为 $N_{comp}$ 个核合成组分的线性叠加。每个组分由四个参数定义：$Y_{e,i}$、$s_i$、$\tau_i$ 以及权重 $w_i$。总丰度由下式给出：$Y_{sample}(Z) = \sum w_i \cdot Y_{model,i}$。
• 推断设置： 作者使用 emcee Python 包对参数的后验概率分布进行采样。似然函数通过比较模型与观测丰度（归一化至锶）的简化卡方度量来定义。
• 先验分布： 研究测试了两种不同的先验：
  1. 不可知论均匀先验（Agnostic Uniform Prior）： 在物理边界内赋予所有参数相等的概率（$Y_e \in [0.005, 0.5]$，特定的 $s$ 和 $\tau$ 范围）。
  2. 基于模拟的先验（Simulation-Informed Prior）： 利用基于 22 种不同模拟（涵盖中子星合并、吸积盘喷发和磁旋转超新星）的 29,000 多条流体动力学轨迹构建的高斯核密度估计（KDE）生成。该先验根据当前模拟中条件的出现频率对条件进行加权。
• 可识别性： 为防止重复计算相同的组分，MCMC 采样器实施了一个排序约束，即组分的电子分数必须递增（$Y_{e,1} \le \dots \le Y_{e,N_{comp}}$）。

关键结果
该框架被应用于近乎完整的 r-过程模板星 HD222925、有限-r 星（HD128279, HD122563）以及太阳 r-过程残余物。

1. 两组分主导地位： 对于 HD222925，仅需两个组分即可高保真地重现观测到的丰度模式。

   • H-组分（高 $Y_e$ 组分）： 较低的 $Y_e$ 条件（$Y_e \approx 0.16-0.18$），产生从第二峰到第三峰的元素（重 r-过程）。
   • L-组分（低 $Y_e$ 组分）： 较高的 $Y_e$ 条件（$Y_e \approx 0.32-0.33$），产生从第一峰到第二峰的元素（轻 r-过程）。
   • 增加组分数量（$N_{comp} = 3, 4, 5$）并不能显著改善与观测值的平均一致性，因为额外的组分主要是在与两组分解相同的参数空间区域进行采样。

2. 条件的普遍性： 相同的两组天体物理条件（H 和 L 组分）主导了 HD222925、HD128279、HD122563 和太阳 r-过程残余物的 r-过程模式。这些恒星之间的差异主要归因于这两个组分的相对权重（$w_2/w_1$），而非存在根本不同的核合成场所。

   • HD222925: $w_2/w_1 \approx 1.7$
   • HD128279: $w_2/w_1 \approx 2.7$
   • HD122563: $w_2/w_1 \approx 7.4$
   • 太阳 r: $w_2/w_1 \approx 4.9$

3. 太阳 r-过程差异： 虽然两组分模型很好地拟合了 HD222925，但它无法重现太阳 r-过程残余物中较轻的中子捕获元素（$Z \approx 30-35$），这些元素的丰度在太阳中显著高于 HD222925。若要拟合太阳第一峰，需要第三个具有更高 $Y_e$（$\approx 0.42$）和更长膨胀时间尺度的组分。

4. 剩余差异： 模型与观测之间的剩余不匹配集中在特定元素上：

   • 观测系统误差： 如 Ag 和 Cd ($Z=47, 48$) 的差异可能源于 3D 效应和非局部热力学平衡（NLTE）效应。
   • 核物理： 第三峰（Re, Os, Ir）附近以及 Th 的差异表明其对核质量模型和裂变势垒具有敏感性。
   • Ga, Ge, As： 这些元素的太阳丰度高于 HD222925，表明在太阳中存在额外的产生通道，或者在提取太阳 r-过程成分时存在问题。

5. 先验比较： 使用不可知论均匀先验和基于模拟的先验所得到的结果几乎相同。这表明重现 HD222925 所需的条件在当前的流体动力学模拟中已是常见情况，并不需要超出模拟参数空间的奇异条件。

意义与主张
本文声称展示了通过贝叶斯推断将恒星丰度观测与大规模核合成计算相结合，以约束 r-过程条件的效力。其主要意义在于：

• 证明了仅通过两个截然不同的天体物理条件（一个轻组分和一个重组分）的叠加，即可解释在贫金属星和太阳系中观察到的稳健 r-过程模式，从而支持了 r-过程的“普遍性”。
• 提供了一个统计框架，能够在不进行任意假设的情况下，量化不确定性并确定必要的组分数量。
• 强调了剩余的差异很可能源于观测系统误差（NLTE、3D 效应）和核物理输入（质量、裂变），而非对基本天体物理场所理解的缺乏。
• 指出了对更大规模、同质分析的 r-过程增强星样本，以及改进的核物理/原子数据进行进一步约束的需求。

作者并未声称已经确定了特定的天体物理场所（例如中子星合并与超新星），而是确定了所需的热力学条件，并指出这些条件存在于各种当前的模拟情景中。