Probing Evolution of Long Gamma-Ray Burst Properties through Their Cosmic Formation History

技术摘要：通过长伽马射线暴的宇宙形成历史探测其性质演化

问题陈述
长伽马射线暴（LGRB）前身星的天体物理学及其性质潜在的宇宙学演化仍是重大争议课题。尽管先前的研究表明，LGRB 率密度（LGRB-RD）在高红移（$z$）处追踪宇宙恒星形成率密度（SFRD），并将其归因于金属丰度的演化，但在低红移处的关系却存在争议。一些研究指出，与 SFRD 相比，LGRB-RD 在低$z$处出现“上升”或过剩，这可能由一类独特的低光度伽马射线暴（LLGRBs）、样本污染或不同的前身星子类引起。此外，现有文献对于 LGRB 是无偏的恒星形成示踪体，还是受金属丰度和初始质量函数（IMF）影响的有偏示踪体，观点相互冲突。本研究旨在利用具有特定 X 射线平台特征的 LGRB 均匀样本，通过调查 LGRB-RD 及其与 SFRD 的映射关系来解决这些差异。

方法论
作者利用了一个由 242 个 LGRB 组成的样本，这些暴发生于 2005 年 3 月至 2023 年 12 月之间，其特征是在尼尔·盖雷尔斯·斯威夫特（Neil Gehrels Swift）天文台观测到的余辉光变曲线中存在 X 射线平台。选择这一特定子样本的动机在于数据的均匀性，以及相较于瞬时辐射研究，其相关性中的弥散度更低。

分析按以下步骤进行：

1. 光度与流量限制：计算平台结束时的观测光度（$L_{a,obs}$）和最小可观测光度（$L_{min,obs}$）。应用流量限制截断（$\log F_{lim} = -11.9$），剔除 17 个暴以减轻 Malmquist 偏差。
2. 非参数修正：为了解决内禀光度演化和选择偏差，研究采用了 Efron & Petrosian（EP）非参数方法。该技术利用 Kendall $\tau$统计量检验光度（$L$）与红移（$1+z$）的独立性。
3. 演化因子确定：分析确定了一个光度演化因子 $(1+z)^k$，其中 $k = 2.1^{+0.5}_{-0.7}$。该因子用于修正累积率密度（$\sigma$）和光度函数（$\Phi$），从而得出修正后的分布（$\sigma_{corr}$和 $\Phi_{corr}$）。
4. 率密度计算：微分 LGRB-RD（$\rho(z)$）由修正后的累积分布推导得出。作者应用修正因子（$F_1$和 $F_2$）以考虑被剔除的暴（由于流量限制或无关联集）以及仪器限制（Swift-BAT 视场、运行时间和运行年份），从而估算真实的 LGRB-RD。
5. 理论比较：将修正后的 LGRB-RD（$\rho_{corr}$）与基于 SFRD（Madau & Dickinson 2014）推导出的理论 LGRB-RD（$\rho_{theor}$）进行比较，采用类似“德雷克方程”的方法。该理论模型纳入了大质量恒星比例、核心坍缩超新星、Ic 型超新星、喷流束流以及金属丰度等因素。

主要贡献

• 更新的光度函数：本研究利用比先前工作（如 Dainotti 等人 2013）更新且更大的样本，更新了具有 X 射线平台的 LGRB 的光度函数（LF）。
• X 射线平台样本的首次 RD 分析：这是第一项专门针对 X 射线平台样本推导 LGRB-RD 并将其与 SFRD 模型直接比较的研究。
• 演化约束：本文提供了对光度演化因子（$k \approx 2.1$）的定量约束，该因子用于使光度和红移分布相互独立，从而修正 Malmquist 偏差。
• 比较分析：该工作将推导出的 $\rho_{corr}$ 与理论模型及先前的观测结果（包括 Dainotti 等人 2024 年的光学平台样本）进行比较，以分离选择效应和演化趋势。

结果

• 光度演化：分析确认了显著的光度演化，最佳拟合演化因子为 $(1+z)^{2.1}$。若无此修正，LGRB-RD 的内禀性质将被掩盖。
• 率密度行为：修正后的 LGRB-RD（$\rho_{corr}$）在大部分红移范围内呈现相对平坦的分布，并在 $(1+z) \approx 2.5$ 附近有一个微弱的峰值。
• 与 SFRD 的不匹配：
  • 低$z$过剩：研究确认了在低红移（$z < 1$）处 LGRB 相对于 SFRD 存在过剩，这与 Yu 等人（2015）、Petrosian 等人（2015）及其他人的发现一致。这种过剩在 $z \sim 0.1$ 处减弱，可能是由于小样本统计所致。
  • 高$z$差异：在高红移处，$\rho_{corr}$ 与大多数 SFRD 模型显示出不匹配，尽管它与 Hopkins & Beacom（2006）的模型部分吻合。SFRD 峰值出现在略高于 LGRB-RD 峰值的红移处。
• 检验演化假设：作者测试了三种将 $\rho_{theor}$ 映射到 $\rho_{corr}$ 的情景：
  1. 无演化：被排除，因为分布不兼容。
  2. 束流角演化：$(1+z)^{1.6}$ 因子使理论曲线的高$z$斜率变平，但无法解释低$z$处的行为。
  3. 简单幂律演化：单一幂律因子 $(1+z)^\delta$ 无法同时匹配观测率密度的低$z$端和高$z$端。负$\delta$值拟合低$z$端，正$\delta$值拟合高$z$端，但两者均无法适用于整个范围。

意义与主张
本文结论认为，“无演化”假设不足以解释 LGRB 形成与恒星形成之间的关系。此外，仅涉及束流角或单一性质幂律演化的简单演化模型，不足以在整个红移范围内将 LGRB-RD 映射到 SFRD。

作者声称，要完全匹配 LGRB-RD 和 SFRD，可能需要同时纳入 LGRB 的多种演化性质。观测到的差异表明，必须考虑选择偏差（如 LLGRB 污染或被错误分类为 LGRB 的并合驱动伽马射线暴）以及复杂的演化因素（如金属丰度依赖性、演化的 IMF 或前身星系统的演化）。研究强调，虽然推导出的 $\rho_{corr}(z)$ 是探测前身星系统和恒星质量黑洞群体的有力工具，但要完全理解 LGRB 的宇宙形成历史，需要更详尽的红移演化函数。作者谦逊地建议，未来的数据（可能借助机器学习进行红移推断）将需要用于进一步解决这些复杂问题。