Balancing bias, baryons, and scale cuts in LSST 3x2pt analysis

技术摘要：在 LSST $3 \times 2\text{pt}$ 分析中平衡偏差、重子效应与尺度截断

问题陈述
第四阶段宇宙学巡天（如鲁宾天文台 LSST）将探测结构形成的非线性机制，在此机制下，来自星系偏差（galaxy bias）和重子反馈（baryonic feedback）的系统效应将占据主导地位。在 $3 \times 2\text{pt}$ 分析（结合星系聚类、宇宙剪切及其互相关）中，这些效应引入了通常难以约束且可能与宇宙学参数表现出简并性的参数化扰动（nuisance parameters）。具体而言，非线性星系偏差与物质谱中重子抑制之间的相互作用提出了挑战：纯线性偏差模型在小尺度（$k \gtrsim 0.1 \, h/\text{Mpc}$）上是不够的，然而更复杂的微扰模型会引入额外的自由度，而这些自由度可能与重子物理学产生简并。核心问题在于如何平衡模型复杂度与尺度截断（scale cuts），以在不引入显著偏差的前提下，最大化对宇宙学参数（$\Omega_m, \sigma_8$）及扩展参数（中微子质量 $M_\nu$）的约束能力。

方法论
作者使用一种名为 MGL（Modified Gravity Lensing）的新型开源贝叶斯推断流水线进行了全面的分析。该研究生成了模拟 LSST 第一年（Y1）和第十年（Y10）规范的数据向量，包括宇宙剪切和星系聚类谱，跨越 20 个对数间隔的 $\ell \in [20, 5000]$ 统计区间。

关键方法组成部分包括：

• 建模框架： 分析采用 BACCO 模拟器 来模拟非线性物质谱和通过“重子化”（baryonification）机制实现的重子反馈。对于星系偏差，对比了三种模型：
  1. 线性偏差（Linear Bias）： 仅在较大尺度上有效。
  2. 混合有效场论（Hybrid Effective Field Theory, HEFT）： 一种结合了拉格朗日偏差展开与 N 体模拟演化的微扰方法，其有效范围高达 $k_{\text{max}} \approx 0.7 \, h/\text{Mpc}$。
  3. 极简偏差（Minimal Bias）： 一种简化的 HEFT 变体，其中高阶偏差项（$b_{s^2}$ 和 $b_{\nabla^2}$）被固定为零。
• 重子处理： 重子反馈通过应用于物质谱的抑制因子 $S_{mm}$ 进行建模。对于星系-星系透镜（$P_{gm}$），作者采用了 Zennaro 等人 (2024) 的近似方法，即 $P_{gm}$ 由 $\sqrt{S_{mm}}$ 进行修正，而 $P_{gg}$ 则依赖于偏差展开的灵活性。
• 分析策略： 作者在由完整 HEFT 模型和标定参数生成的模拟数据上进行似然分析。他们测试了各种分析假设，包括模型失配的情况（例如，使用线性或极简偏差模型来分析由 HEFT 生成的数据）以及包含重子参数的不同情况。
• 指标： 通过优值（Figure of Merit, FoM）来量化约束能力，并通过偏差图（Figure of Bias, FoB）来衡量与标定值的偏差。

主要贡献

1. 开发 MGL 流水线： 本文介绍了一种独立的、开源的 MGL 流水线，旨在服务于第四阶段巡天。该流水线利用模拟器实现功率谱的高效计算，并支持使用 nautilus 采样器的 MCMC 分析。
2. 偏差模型的系统比较： 该工作在 LSST Y1 和 Y10 的背景下，对线性、HEFT 和极简偏差模型进行了直接比较，特别评估了它们针对重子反馈和尺度截断的稳健性。
3. 简并性分析： 研究明确调查了高阶偏差参数与重子反馈之间的简并性，以及偏差建模与质量中微子探测之间的相互作用。

结果

• 尺度截断与偏差模型： 线性偏差模型仅在 $k_{\text{max}} \approx 0.1 \, h/\text{Mpc}$ 以下能提供对 $\Omega_m$ 和 $\sigma_8$ 的无偏约束（FoB $< 0.5\sigma$）。若要扩展到更小尺度，则需要微扰方法。
• 极简偏差的效能： “极简偏差”模型（将 $b_{s^2}$ 和 $b_{\nabla^2}$ 固定为零）的表现与完整的 HEFT 模型相当。只要标定偏差参数处于标准范围内，即使在 $k_{\text{max}} = 0.7 \, h/\text{Mpc}$ 时，它对 $\Omega_m$ 和 $\sigma_8$ 的约束仍保持无偏。这表明，对于标准宇宙学约束，可能不需要 HEFT 的全部复杂度，从而为减少参数空间提供了一条路径。
• 重子-偏差相互作用： 高阶偏差参数可以模仿星系-星系功率谱中的重子反馈抑制效应。然而，重子反馈无法完全复制由高阶偏差项产生的行为范围。忽略星系-星系透镜部分的重子修正不会显著使 $\Omega_m$ 或 $\sigma_8$ 产生偏差，但可能会改变推断出的重子参数（$\log_{10} M_c$）。
• 中微子质量探测： 当扩展到非线性尺度（$k \geq 0.3 \, h/\text{Mpc}$）时，对于 Y1 和 Y10 均可实现总中微子质量（$M_\nu$）的探测。然而，特定的推断值在所有可能的场景下并不稳健。当使用极简偏差模型分析具有较高阶偏差参数的模拟数据时，会对推断的 $M_\nu$ 产生显著偏差（对于 Y10 高达 $\sim 2\sigma$）。这表明中微子质量、高阶偏差项以及投影效应之间存在强烈的简并性。
• Y1 与 Y10 对比： 从 Y1 转向 Y10，在 $k_{\text{max}} = 0.7 \, h/\text{Mpc}$ 时，FoM 提升了约 2.3 倍。虽然 $\Omega_m$ 和 $\sigma_8$ 保持稳健，但简化偏差模型下 $M_\nu$ 约束的脆弱性依然存在。

意义与主张
本文声称，对于标准 $\Lambda$CDM 参数（$\Omega_m, \sigma_8$），极简偏差参数化足以从 LSST 数据中提取无偏约束，直至 HEFT 的有效限度（$k_{\text{max}} = 0.7 \, h/\text{Mpc}$），即使在存在重子反馈的情况下也是如此。这允许纳入准非线性尺度，而无需承担完整 15 项偏差展开带来的计算和统计惩罚。

然而，作者谨慎地指出，这一结论取决于特定的标定偏差参数值。虽然 LSST 应该能够对非零中微子质量进行明确探测，但精确的显著性和探测值高度依赖于偏差建模的准确性。研究强调，简化的建模方法（如极简偏差）可能对标准宇宙学研究结果稳健，但如果潜在的星系偏差物理特性偏离了假设的标定值，则可能对中微子质量等扩展参数引入显著偏差。因此，作者主张在采用简化模型之前，应针对特定的科学案例进行仔细测试。