Constraining Dark Acoustic Oscillations with the High-Redshift UV Luminosity Function

技术摘要：利用高红移紫外亮度函数约束暗声学振荡

问题陈述
虽然标准 $\Lambda$CDM 模型假定暗物质（DM）是冷且无碰撞的，但各种具有良好理论依据的暗部门扩展模型预测了“暗声学振荡”（Dark Acoustic Oscillations, DAOs）。这些振荡源于在早期宇宙中一部分暗物质与一种暗辐射（DR）物种发生相互作用，导致在小尺度上产生功率谱 $P(k)$ 的抑制和阻尼振荡。现有的来自宇宙微波背景（CMB）的约束在长尺度上非常严格，但对于强 DAO 而言，对相互作用暗物质的比例限制在 $\lesssim 5\%$。相反，来自莱曼-$\alpha$ 森林（Lyman-$\alpha$ forest）的约束对于弱 DAO 或完全相互作用的模型非常有效，但在区分强 DAO 特征与温暗物质（WDM）效应方面面临挑战。因此，需要对高红移（$z \sim 3-9$）处的小尺度结构进行模型无关的约束，因为在这些红移处，非线性结构演化和 DAO 对晕质量函数（HMF）的具体印记可以被探测到。

方法论
作者开发了一个流程，利用高红移紫外亮度函数（UVLF）来约束 DAO 参数。该方法由四个部分组成：

1. 唯象传递函数模型：
   为了避免为每一个马尔可夫链蒙特卡罗（MCMC）扫描点求解玻尔兹曼方程带来的计算开销，作者引入了一个新的四参数唯象传递函数 $T(k)$，该函数相对于 $\Lambda$CDM。该模型捕捉了 DAOs 的抑制和阻尼振荡，并允许存在小于 1 的相互作用暗物质比例（$f$）。其参数为：

   • $f$：经历 DAOs 的暗物质比例。
   • $k_{\text{peak}}$：第一个 DAO 峰值的波数。
   • $k_{\text{damp}}$：振荡的阻尼尺度。
   • $A$：振荡的振幅。
     该模型经过修改后的玻尔兹曼代码 CLASS-aDM（用于原子暗物质）和 CLASS（用于相互作用暗物质）生成的功率谱校准，证明了这些参数可以直接映射到物理量，如暗声速视界和扩散阻尼尺度。

2. 扩展压强-谢克雷斯塔（EPS）形式体系与 N 体模拟校准：
   为了从线性功率谱预测 HMF，作者采用了扩展压强-谢克谢斯塔（EPS）形式体系。考虑到标准窗口函数（例如实空间顶帽函数）对于截断功率谱失效，他们使用了一种“平滑-$k$”窗口函数。至关重要的是，他们使用 GIZMO 代码配合 MUSIC 初始条件进行了 13 组暗物质唯 N 体模拟。这些模拟涵盖了不同的相互作用暗物质比例（$f=0.1, 0.5, 1$）和 DAO 尺度（$k_{\text{DAO}} = 5, 20, 66, 190 \, h/\text{Mpc}$）。这些模拟被用于校准 EPS 参数（$q, \beta, c$），以确保半解析 HMF 在 $z=12$ 到 $4$ 的红移范围内能准确重现模拟结果。

3. 晕-星系连接：
   利用 GALLUMI 代码对 HMF 与可观测的 UVLF 之间的联系进行建模。该流程假设紫外亮度是由恒星形成率（SFR）驱动的，而 SFR 平均而言遵循与晕吸积率相关的双幂律关系。该模型包含八个天体物理学上的疑难参数（描述 SFR-晕质量关系的归一化、转折标度、演化以及量级中的随机性），通过对这些参数进行边缘化处理，以解释高红移星系形成物理过程中的不确定性。

4. 数据分析：
   作者结合了覆盖红移 $z=3$ 到 $9$ 的哈勃空间望远镜（HST，包括空白场和透镜场）、詹姆斯·韦伯空间望远镜（JWST）、苏鲁比望远镜（Subaru）以及加拿大-法国-希马拉雅望远镜（CFHT）的 UVLF 测量数据。他们执行了贝叶斯 MCMC 扫描和频率派剖面似然分析来约束 DAO 参数，并将 $\Lambda$CDM 参数固定为 Planck 2018 的最佳拟合值。

核心贡献

• 新型唯象模型： 引入了一种灵活的传递函数模型，能够容纳次级相互作用暗物质比例，弥补了特定粒子物理模型（如原子暗物质）与通用约束之间的鸿沟。
• 模拟校准： 首次利用专门的 N 体模拟对 DAO 模型的 EPS 形式体系进行校准，验证了在存在阻尼振荡的情况下，半解析方法预测 HMF 的有效性。
• 模型无关的约束： 推导出了关于 DAO 参数空间（$f$ 对比 $k_{\text{peak}}$）的约束，这些约束独立于特定的暗部门粒子物理学，允许在各种模型（如 IDM、原子暗物质）中进行轻松的重新解释。

结果

• 对 $f$ 和 $k_{\text{peak}}$ 的约束： 分析显示相互作用暗物质的比例与 DAO 的尺度之间存在强相关性。
  • 如果相互作用暗物质的比例显著（$f \gtrsim 0.07$），则第一个 DAO 峰值的波数必须满足 $k_{\text{peak}} \gtrsim 50 \, h/\text{Mpc}$ 才能与 UVLF 数据保持一致。
  • 反之，如果 $k_{\text{peak}}$ 很小（表示更早的解耦或更强的相互作用），则相互作用暗物质的比例被限制为 $f \lesssim 0.07$。
• 与其他探测手段的比较：
  • 对于高比例的相互作用暗物质（$f \sim 1$），UVLF 约束在 $k_{\text{peak}}$ 上比 CMB 衍生的界限提高了约一个数量级。
  • 对于极小尺度（$k_{\text{peak}} \lesssim 1 \, h/\text{Mpc}$），CMB 约束在相互作用暗物质的丰度（$f \lesssim 5\%$）方面仍然更为严格。
  • 在相互作用暗物质（IDM）模型中，对暗物质-重子散射截面的约束比 CMB 约束强五个数量级以上，并且与来自银河系卫星丰度的界限具有竞争力。
• 未受约束的参数： 分析发现，在对其他参数进行边缘化处理时，无法对 DAO 振幅（$A$）或阻尼尺度比（$k_{\text{damp}}/k_{\text{peak}}$）给出显著约束，因为 UVLF 的敏感性主要由功率谱的初始抑制而非振荡特征所驱动。

意义
本文表明，高红移紫外亮度函数测量是约束小尺度暗物质物理的一种强大且互补的手段，可与 CMB 和莱曼-$\alpha$ 森林数据相媲美。通过成功地利用 N 体模拟校准 EPS 形式体系，并采用灵活的晕-星系连接模型，作者为相互作用暗物质经历声学振荡的比例建立了稳健的限制。这些结果将相互作用暗物质的解耦时期推向了更早的时期（更高的 $k_{\text{peak}}$），或者要求在特定的相互作用尺度下，相互作用比例必须是次级的（$f < 7\%$）。这项工作强调了高红移可观测物理在区分不同暗部门范式方面的效用，并为未来使用薇拉·鲁宾天文台、南希·格雷斯·罗曼空间望远镜及平方公里阵列等数据的约束奠定了基础。