Unbiased analysis of primordial non-Gaussianity: the multipoles of the full relativistic power spectrum

技术摘要：全相对论功率谱多极子的无偏分析：针对原始非高斯性的研究

问题陈述
下一代光谱巡天（如 Euclid、MegaMapper、DESI）旨在通过星系功率谱中的尺度相关晕偏置（scale-dependent halo bias）来约束局部型原始非高斯性（PNG），其参数为 $f_{\text{NL}}$。然而，随着这些巡天探测的体积增大及进入超大尺度，标准的牛顿力学红移空间畸变（RSD）处理方法已不再适用。广义相对论投影效应（特别是积分贡献，如引力透镜收敛、时间延迟和集成萨克斯-沃尔夫效应/ISW效应）以及宽间距修正（破坏局部平移不变性）会与局部 PNG 信号产生简并。忽略这些效应会导致对 $f_{\text{NL}}$ 的推断产生偏差。此外，巡天光度函数的不确定性（演化和放大偏差）引入了额外的简并，从而降低了约束精度。

方法论
作者开发了一个一致的理论框架，用于计算包含所有相关光锥效应的完整相对论 3D 功率谱 $P(k)$ 的多极子。

• 形式化方法： 基于用于功率谱多极子的 Yamamoto 估计量，作者定义了针对单个视线（LOS）$d$ 的“局部”功率谱 $P_{\text{loc}}(k, \mu; d)$。他们明确地将星系数密度建模为标准 Kaiser 项（$S$）、局部相对论投影项（$NI$）和积分项（$I$）之和。
• 积分效应： 一个核心重点是推导功率谱中积分-积分交叉项（$I \times I$）和积分-局部交叉项（$I \times S$）的贡献。不同于以往的角向分析，这些项直接在笛卡尔傅里叶空间中进行计算。作者采用了宽间距展开（利用 $\epsilon = x_{12}/d$）来处理由积分效应导致的破坏局部平移不变性的问题，这些效应取决于视线选择（由参数 $t$ 表征，其中 $t=0$ 为端点，$t=1/2$ 为中点）。
• 宽间距修正： 该框架还包括了用于多示踪器分析的协方差矩阵的宽角和径向演化修正，这是一项新颖的计算。
• 预报： 作者采用 Fisher 矩阵和马尔可夫链蒙特卡洛（MCMC）分析，为各种巡天配置（DESI BGS、Euclid H$\alpha$、Roman、MegaMapper LBG、SKAO2）预报 $f_{\text{NL}}$ 的约束。他们在线性偏置（$b_1$）、演化偏置（$b_e$）和放大偏差（$Q$）代表不确定性的参数化参数上进行了边缘化处理。
• 多示踪器方法： 为了减轻光度函数不确定性（演化和放大偏差）与 $f_{\text{NL}}$ 之间的简并，作者实施了一种“亮-暗”划分方案，根据亮度将星系样本分为两类，以利用交叉相关和奇宇称多极子。
• 实现： 所有计算均通过 Python 代码 CosmoWAP 及其配套的 Mathematica notebook MathWAP 实现。

主要贡献

1. 在 3D 功率谱中推导 $I \times I$ 项： 本文首次在 3D 傅里叶功率谱多极子中直接计算了 $I \times I$ 贡献（引力透镜-引力透镜、时间延迟-时间延迟等），扩展了 Limber 近似的应用范围。
2. 视线依赖性与奇多极子： 作者证明了积分贡献会破坏局部平移不变性，导致对视线选择（$t$）的依赖。至关重要的是，对于端点视线选择（$t \neq 1/2$），这些项会在自功率谱中产生非零的、虚数形式的奇宇称多极子（例如偶极子），这是 3D 背景下的一个新结果。
3. 偏差量化： 本研究量化了当忽略相对论和宽间距效应时，对 $f_{\text{NL}}$ 推断造成的偏差。
4. 新颖的协方差计算： 作者提出了一个新的解析计算方法，用于计算完整的多示踪器功率谱多极子协方差，其中包括宽间距修正和积分效应，这对于大尺度上的精确误差估计至关重要。
5. 多示踪器缓解作用： 研究表明，通过“亮-暗”划分的多示踪器分析可以部分打破 $f_{\text{NL}}$ 与光度函数不确定性之间的简并，从而提高约束精度。

结果

• 对 $f_{\text{NL}}$ 的偏差： 忽略积分和相对论效应会导致对 $f_{\text{NL}}$ 最佳拟合值的显著偏差。对于类 Euclid 的 H$\alpha$ 巡天，偏差约为 $3\sigma$；而对于类 MegaMapper 的莱曼断裂星系（LBG）巡天，偏差可达 $\sim 20\sigma$。偏差的主要驱动力是 $I \times I$ 项巨大的正向贡献，它增强了大尺度的功率，抵消了来自局部相对论项的负向抑制。
• 非线性效应： $I \times I$ 贡献对非线性功率谱建模具有敏感性，即使在轻微非线性尺度（$k \sim 0.05\, h/\text{Mpc}$）下也会表现出 $\sim 10\%$ 的功率提升，若不加以考虑，会增加对 $f_{\text{NL}}$ 的偏差。
• 光度函数不确定性： 在包含相对论效应时，光度函数的不确定性（通过对 $b_e$ 和 $Q$ 进行边缘化）会导致 $f_{\text{NL}}$ 的误差比固定偏置情景增加：Euclid 增加约 $30\%$，MegaMapper 增加约 $20\%$。
• 多示踪器改进： 利用“亮-暗”划分可以将 $f_{\text{NL}}$ 的不确定性降低 $15\text{--}20\%$（相比于具有相同相对论建模的单示踪器分析）。这种改进源于交叉相关中的奇宇称多极子有助于约束演化和放大偏差，从而打破了与 $f_{\text{NL}}$ 的简并。
• 可探测性： 据预测，引力透镜对功率谱的贡献在 Euclid、SKAO2 和 MegaMapper 巡天中是可探测的。然而，时间延迟和 ISW 贡献在 3D 傅里叶功率谱多极子中不太可能被探测到。

意义
本文确立了对于第四阶段及未来的星系巡天，采用“非牛顿”方法对于进行无偏的原始非高斯性分析是强制性的。作者认为，忽略积分相对论效应和宽间距修正会导致错误的探测结论或严重的 $f_{\text{NL}}$ 偏差。此外，他们展示了虽然这些效应引入了新的与光度函数之间的简并，但可以通过多示踪器技术进行有效管理。这项工作提供了一个完整的、公开的计算框架（CosmoWAP），用于将这些复杂的轻锥效应纳入标准的功率谱分析中，确保未来对早期宇宙物理的约束能够不受相对论系统误差的影响。