技术摘要：应用场级推断重建宇宙最大尺度——以 Quaia 类星体目录为例

1. 问题陈述

绘制宇宙的大尺度结构（LSS）对于提取基础物理学至关重要，然而依赖于统计量（如功率谱、双谱）的传统分析方法存在数据压缩次优的问题，未能充分利用高维信息。类星体是能够探测高达红移 $z \sim 7$ 并覆盖超大尺度的独特示踪体，但其分析过程受到低数密度（高泊松噪声）、前景污染以及红移不确定性（特别是针对光度样本）的复杂影响。

最近发布的 Quaia 类星体目录（基于 Gaia 和 unWISE 数据）提供了全天覆盖范围（$f_{sky} \approx 0.71$）和广泛的红移范围（$0.084 < z \lesssim 4.0$）。然而，要利用该数据集进行完整的三维物质分布重建，需要一种能够同时处理巡天选择效应、红移空间畸变、类星体偏置以及示踪体固有稀疏性的方法。现有方法往往难以在不丢失信息的情况下，将这些系统误差整合进一个连贯的物理模型中。

本研究应用了 BORG（基于伽利略起源的星系贝叶斯重建）算法，这是一种场级推断（FLI）框架，将其应用于 Quaia 类星体目录。与通过压缩数据得到统计量的传统方法不同，BORG 执行的是对整个观测星系场的贝叶斯前向建模。

2. 方法论

本工作将 BORG（基于伽利略起源的星系重建）算法应用于 Quaia 类星体目录。与通过将数据压缩为统计量的方法不同，BORG 对观测到的整个星系场进行物理驱动的前向建模。

2.1 数据样本

分析使用了来自 Quaia 目录的两个不同样本（Storey-Fisher 等人，2024）：

• Quaia Clean：一个星等受限样本（$G < 20.0$），包含约 $755,850$ 个源。
• Quaia Deep Cut：一个扩展样本（$G < 20.5$），包含约 $1,295,502$ 个源。该样本采用了更严格的角向选择，排除了大、小麦伦尼星际云（LMC/SMC）附近的区域，并施加了大于 $0.5$ 的完备性阈值以减轻前景污染。

两个样本都被划分为八个径向分箱，以探测不同的宇宙时期，同时保持统计效力。红移被转换到 CMB 静止坐标系中，以消除由太阳特有运动引起的偶极子效应。

2.2 前向模型

BORG 算法采用基于物理的前向模型，将初始条件演化为现今的可观测量：

1. 初始条件：通过由原始功率谱（固定为 Planck 2018 宇宙学参数）调制的高斯白噪声场来生成初始引力势。
2. 引力演化：使用拉格朗日摄动理论（LPT）（一阶）来演化密度场，模拟从初始条件到现今的结构形成。考虑到当前研究的大尺度和粗分辨率（$39.1 \, h^{-1}\text{Mpc}$），这种方法在计算上是非常高效的。
3. 观测效应：模型纳入了红移空间畸变、光锥效应以及巡天窗口函数（包括径向选择函数和角向完备性掩模）。
4. 偏置建模：使用泊松幂律偏置模型来建立底层暗物质密度与观测到的类星体分布之间的联系。

2.3 推断引擎

该算法利用**哈密顿蒙特卡洛（HMC）**结合切片采样器，探索初始条件和晚期场的高维后验分布。这生成了一组物理上合理的实现，包括：

• 原始初始条件。
• 现今的暗物质密度场。
• 大尺度速度场（径向和绝对速度）。
• 预测的类星体分布。

2.4 验证策略

为了验证重建结果，作者执行了以下操作：

• 内部一致性检查：将恢复的功率谱和双谱与理论预测（CLASS）进行比较，并检查初始条件中的非高斯性。
• 外部互相关：一个关键的验证步骤是将 BORG 推断的物质密度场与 Planck 2018 CMB 引力透镜收敛图（$\kappa_{Planck}$）进行互相关。这测试了重建的大尺度结构是否与通过 CMB 引力透镜观测到的集成物质分布相一致。

3. 主要贡献与结果

3.1 前所未有的重建体积

该研究产生了迄今为止在共动体积上最大的可观测宇宙场级重建。

• 体积：$(10 \, h^{-1}\text{Gpc})^3$。
• 分辨率：最大空间分辨率为 $39.1 \, h^{-1}\text{Mpc}$（体素大小）。
• 范围：重建涵盖了整个 Quaia 巡天体积，提供了初始条件、现今暗物质和速度场的后验图谱。

3.2 推断场

• 密度场：重建的暗物质密度场在定性上与观测到的类星体分布表现出强一致性。它捕捉到了大尺度相干结构（纤维状结构和空洞），并通过引力耦合延伸至巡天边界之外，且在示踪体密度较低的区域，不确定性自然增加。
• 初始条件：推断的原始涨落与高斯统计一致（均值 $\approx 0$，偏度和峰度与零一致），证实了前向模型的有效性，并排除了由推断过程引入的人为非高斯性。
• 速度场：重建的速度场揭示了相干的大尺度体流，其中流入和流出过程追踪着底层的密度分布。信噪比遵循巡天窗口函数，突显了算法处理观测限制的能力。

3.3 统计验证

• 功率谱与双谱：推断的物质功率谱在广泛的尺度范围内与来自 CLASS 的理论预测一致。归一化双谱与零一致，表明不存在伪非线性贡献。
• 子体积一致性：在重建体积的子区域中计算的功率谱显示出统计一致性，排除了显著的空间模型失配或局部污染。

3.4 CMB 引力透镜互相关

最重要的外部验证是检测到了 BORG 重建密度场与 Planck CMB 引力透镜图之间的互相关信号。

• 显著性：对于分辨率为 $N=256$ 的 Quaia Deep Cut 样本，检测显著性达到 $\sim 4\sigma$（具体而言，通过振幅拟合为 $3.96\sigma$，通过相位对齐为 $3.0\sigma$）。
• 振幅：最佳拟合振幅参数 $A$ 为 $0.70 \pm 0.17$，在不确定性范围内与理论预期（$A=1$）一致。
• 意义：这一检测证实了 BORG 框架成功恢复了由类星体追踪的真实大尺度物质分布，通过独立宇宙学探针验证了重建结果。

4. 意义与主张

本文声称建立了一个利用类星体巡天进行场级宇宙学的稳健框架。其主要意义在于：

1. 首次应用于类星体：这是首次将场级推断应用于类星体数据集，证明了稀疏、多噪的示踪体可以用于重建三维物质分布和初始条件。
2. 处理系统误差：该工作表明，完全前向建模的层次贝叶斯方法可以有效地处理复杂的观测系统误差（选择效应、红移误差、偏置），而不会造成信息损失。
3. 验证大尺度物理：$\sim 4\sigma$ 的 CMB 引力透镜信号检测为重建场准确反映大尺度结构提供了独立确认。
4. 为未来宇宙学奠定基础：作者认为，该框架能够实现对超大尺度物理学的精密研究，包括对原始非高斯性（$f_{NL}$）的约束、对广义相对论的检验以及对集成萨克斯-沃尔夫（ISW）效应的研究。重建的场作为局部宇宙的“数字孪生”，为约束模拟和进一步的宇宙学分析提供了数据产品（后验图谱）。

作者指出，目前的重建分辨率受限于泊松噪声和光度光谱红移不确定性，这阻碍了在较小尺度上的可靠推断。然而，结果证实，场级推断是一种可行且强大的工具，可以将稀疏的类星体观测转化为统计稳健的三维宇宙大尺度结构图谱。