[论文解读] AlterBBN v2: A public code for calculating Big-Bang nucleosynthesis constraints in alternative cosmologies

B. Bilki, Jérémy Auffinger+2 more|arXiv (Cornell University)|Jun 28, 2018

Particle physics theoretical and experimental studies参考文献 19被引用 15

一句话总结

AlterBBN v2 是一款公开的开源 C 语言代码，用于在标准与替代宇宙学模型下精确计算原初核合成的同位素丰度。该代码支持快速、精确且用户友好的轻元素产额（D、He-4、Li-7）计算，内置误差与相关性分析功能，可处理多种宇宙学情景，如修改的哈勃膨胀率、额外中微子、衰变标量场及 WIMPs 等，所有计算通过 OpenMP 实现并行化，并统一采用能量单位为 GeV。

ABSTRACT

We present the version 2 of AlterBBN, an open public code for the calculation of the abundance of the elements from Big-Bang nucleosynthesis. It does not rely on any closed external library or program, aims at being user-friendly and allowing easy modifications, and provides a fast and reliable calculation of the Big-Bang nucleosynthesis constraints in the standard and alternative cosmologies.

研究动机与目标

• 提供一款快速、可靠且用户友好的开源工具，用于在标准与替代宇宙学模型下计算原初轻元素丰度。
• 通过在修改的哈勃膨胀率、额外中微子种类及熵注入等条件下计算核合成产额，实现对宇宙学参数的精确约束。
• 采用先进的数值方法求解刚性微分方程，并整合误差与相关性分析，以支持对观测数据的约束评估。
• 支持新型宇宙学情景，如有效再加热、衰变标量场及 WIMP 相互作用，扩展标准 BBN 的适用范围。
• 统一使用 GeV 为单位，并通过 OpenMP 实现并行计算，提升多参数扫描中的计算性能。

提出的方法

• 该代码通过一组耦合的常微分方程模拟大爆炸核合成过程中核素丰度的演化，采用自适应的 Runge-Kutta 方法（阶数 2–5）并配合步长控制求解。
• 采用线性化与隐式积分技术实现刚性方程求解器，以应对 BBN 中反应时间尺度跨度极大的问题。
• 核反应速率数据来源于 REACLIB 数据库，反应网络完整扩展并经过验证，包含 22 种核素与 102 个反应。
• 误差与相关性矩阵的计算方法与文献 [1998PhRvD..58f3506F, 2016JHEP…11..097A] 类似，支持基于 χ² 的观测数据对比。
• 代码支持多种宇宙学拓展情景：修改的哈勃参数（EoS）、额外中微子种类、熵注入、衰变标量场及 WIMP 相互作用。
• 单位统一为 GeV 以与 SuperIso Relic 软件包保持一致，同时支持 OpenMP 并行化，提升多参数研究中的性能扩展能力。

实验结果

研究问题

• RQ1 在重子-光子比、中子寿命及中微子简并度发生偏离时，对预测的原初轻元素丰度有何影响？
• RQ2 与标准 BBN 相比，替代宇宙学模型（如修改的哈勃膨胀率或有效再加热）在多大程度上改变了 D、He-4 与 Li-7 的预测产额？
• RQ3 AlterBBN v2 在基于 χ² 统计检验的情况下，对 4He、D 与 7Li 丰度的观测约束再现精度如何？
• RQ4 引入误差与相关性矩阵对 BBN 参数空间中约束评估的可靠性有何影响？
• RQ5 随着核素网络复杂度与宇宙学参数数量的增加，该代码的性能表现如何？

主要发现

• 代码以高精度计算标准 BBN 产额：当 failsafe=12 时，Yp = (2.472 ± 0.003)×10⁻¹，²H/H = (2.463 ± 0.038)×10⁻⁵，⁷Li/H = (5.370 ± 0.352)×10⁻¹⁰。
• 误差与相关性分析的引入使得与观测数据的 χ² 比较更加稳健，使用最新观测数据：Yp = 0.2449 ± 0.0040，²H/H = (2.527 ± 0.030)×10⁻⁵。
• 计算时间随 failsafe 等级显著变化：使用 GCC 与 OpenMP 时，failsafe=12 的计算耗时为 72.1262 秒，凸显其在详细参数扫描中的可扩展性。
• OpenMP 并行化可将复杂计算（如 failsafe=12）的运行时间相比单线程执行减少约 20%，尤其在激活 REACLIB 时效果更明显。
• 代码成功复现了已知的 BBN 约束，包括 Cooke 等人（2017）与 Aver 等人（2015）的保守参数集，并实现了统一的 χ² 评估。
• 对新情景（如衰变标量场与 WIMP 相互作用）的实现，展示了该代码在超越标准模型宇宙学研究中的灵活性。