Angular bispectrum of matter number counts in cosmic structures

技术摘要：宇宙结构中物质数计数数的角双谱

问题陈述
两点相关函数（功率谱）仍然是约束标准 $\Lambda$CDM 模型的主要工具。然而，随着大尺度结构（LSS）演化进入弱非线性机制，高阶统计量——特别是三点相关函数及其傅里叶对应物，即双谱（bispectrum）——对于探测由非线性引力演化和原初非高斯性（PNG）引起的傅里叶相位相关性和模式耦合变得至关重要。

目前对星系双谱的分析严重依赖于局部平坦天空近似、Limber 近似以及牛顿动力学。对于覆盖广阔区域并利用有限红移分箱（测距法/tomography）的下一代巡天项目（如 Euclid、DESI、SPHEREx），这些近似将会失效。此外，虽然相对论效应（投影、动力学广义相对论和辐射）在理论上已知存在，但在不使用 Limber 近似且针对有限红移分箱设置的情况下，其对角双谱的精确影响尚未得到充分计算。此外，由于数值伪影以及牛顿项的主导地位，从模拟中提取这些微妙的相对论信号具有挑战性。

方法论
作者利用二阶摄动理论，提出了物质数计数的全天角双谱计算方法。该方法的技术选择特征如下：

1. 框架： 计算采用了“测距球面谐波”（TSH），利用窗口函数来定义有限红移分箱（$1.75 \le z \le 2.25$ 以及 $0.55 \le z \le 0.65$）。这种方法避免了 Limber 近似和平坦天空假设。
2. 摄动展开： 对数计数摄动 $\Delta$ 进行二阶展开。作者根据弱场展开参数 $H/k$ 对各项进行组织：
   • 牛顿项 ($O((H/k)^0)$)： 包括线性密度、红移空间畸变（RSD）以及二次项（密度-密度、密度-速度耦合）。
   • 相对论投影效应 ($O(H/k)$)： 在弱场展开的一阶内包含领先阶的投影效应（多普勒、萨克斯-沃尔夫效应等）。
   • 动力学 GR 和辐射效应 ($O((H/k)^2)$ 及更高阶)： 包括非线性动力学 GR 校正和辐射校正（遵循 Tram 等人 2.2016 的研究）。
   • 耦合项： 包括投影与动力学效应之间高达 $O((H/k)^3)$ 的耦合。
3. 排除项： 为了简化并实现直接的一致性检查，作者忽略了示踪体偏置（tracer bias）、视线方向积分效应（引力透镜、集成萨克斯-沃尔夫效应）以及观测者位置项。由于复杂性和文献中缺乏共识，他们也省略了二阶中的次领先投影效应。
4. 数值实现： 作者扩展了 Assassi 等人 (2017) 高效的数学框架，利用 FFTLog 分解和超几何函数求值，将问题简化为一维积分。他们提供了公开代码 ang_bispec，用于计算 14 个不同的树级（tree-level）双谱贡献项。
5. 验证： 理论预测与相对论 $N$ 体模拟（Montandon 等人 2023）进行了对比。为了确保公平比较，作者从模拟数据中移除了积分投影效应，并对理论预测应用了分箱和高斯平滑，以匹配模拟估计器（Bucher 等人 2016）。

主要贡献

• 首次全天计算： 本工作提供了在二阶摄动理论下，针对有限红移分箱且不使用 Limber 近似的角双谱的全天计算。
• 新型牛顿项计算： 即便是针对该特定设置（有限分箱、无 Limber）下的牛顿部分双谱，也被报告为具有新颖性。
• 全面的相对论包含： 作者在统一框架内，最高达二阶，包含了所有的动力学 GR 效应、辐射校正以及领先的相对论投影效应。
• 代码发布： ang_bispec 代码已公开发布，扩展了以往的计算，涵盖了所有主导的牛顿项和相对论效应。

结果
作者评估了两个红移分箱的双谱，并对比了各项的贡献：

1. 牛顿项的主导地位： 正如预期，在整个谱范围内，对于两个红移度，牛顿项（密度、RSD、二次项）的贡献通常比相对论信号大一个或多个数量级。
2. $z \approx 2$ 处的相对论层级：
   • 投影效应和动力学 GR 效应在长波尺度上的振幅相当。
   • 意外发现： 在挤压极限（squeezed limit）下，辐射效应在相对论信号中占据主导，超过其他相对论项达一个数量级。
3. $z \approx 0.6$ 处的相对论层级： 恢复了预期的层级结构，即投影效应占主导。然而，动力学校正（来自非线性演化）仍然不可忽视，其规模仅比投影效应小 2–3 倍。
4. 模拟比较：
   • 总双谱： 理论结果与模拟测量值在总双谱方面吻合较好，除在模拟分辨率限制功率的小尺度处外。
   • 相对论信号： 在将“纯”相对论信号（GR 与辐射项之和）与模拟测量值（通过从相对论模拟中减去牛顿模拟得出）进行比较时，作者发现存在差异。在高红移处，理论总和比模拟测量值小约 5 倍。
   • 对差异的解释： 作者将这一差距归因于模拟中的数值伪影，特别是残余的规范效应（gauge effects，即 Poisson 规范与 N-body 规范之间的不匹配）和回标定效应（backscaling effects，即辐射如何在牛顿模拟中被包含）。他们得出结论，这些数值误差在振幅上与物理性的动力学相对论效应相当。
   • 平滑效应： 分析强调，平滑处理（对于检测微小的相对论信号是必要的）会引入系统偏差，将小尺度的数值误差泄漏到大尺度，这可能解释了在低红移比较中观察到的偏移。

意义与主张
本文声称建立了一个稳健的理论流程，用于分析光度巡天中物质数计数的角双谱，弥合了引力透镜与星系聚类观测之间的鸿隙。

• 振幅验证： 本工作证实了 Montandon 等人 (2023) 通过模拟得出的相对论效应振幅估计，同时提供了模拟难以轻易分离的逐项分解。
• 动力学效应的相关性： 与认为投影效应是唯一显著相对论贡献的假设相反，作者证明了在高质量红移处，动力学 GR 和辐射效应与投影效应处于同一数量级，并且在较低红移处依然显著（差距在 2 倍以内）。
• 对模拟提取的警示： 一个重要的发现是，当前相对论模拟中的数值伪影（规范效应和回标定效应）与正在寻找的物理信号本身一样大。这表明，使用目前的模拟减法方法提取“纯”相对论信号需要极其谨慎。
• 未来用途： 作者认为，他们的框架对于未来的下一代巡天（如 Euclid）至关重要，以便将原初非高斯性与内在非高斯性区分开来，并测试非线性部门的引力理论，前提是最终能够纳入星系偏置和视线方向积分等简化项。