Improving constraints on primordial non-Gaussianity from Quaia with a new cosmological observable: angular redshift fluctuations

技术摘要：通过具有角红移涨落的 Quaia 改进对原初非高斯性的约束

问题与动机
对原初非高斯性（Primordial Non-Gaussianity, PNG）的探测——以 $f_{\text{NL}}$ 参数化——是检验暴胀模型的关键测试，其中非零值将表明偏离高斯初始条件的特征。虽然目前最严格的约束来自于 CMB 双谱（例如 Planck 2018），但大尺度结构（LSS）通过尺度相关偏差效应提供了一种互补的探测手段。该效应由 Dalal 等人（2008）提出，在局部型 PNG 下，会在星系偏差中引入特征性的 $1/k^2$ 信号。然而，由于需要对最大尺度上的观测系统误差进行精确控制，从 LSS 中测量这一效应具有挑战性。

近期的研究利用二维角功率谱以及与 CMB 透镜效应的互相关来约束 $f_{\text{NL}}$。本研究旨在通过引入一种新的宇宙学可观测量：角红移涨落（Angular Redshift Fluctuations, ARFs），来探讨提升这些约束的可能性。ARFs 测量给定壳层内平均红移的二维角偏差。理论上，ARFs 对星系偏差敏感，因此对 PNG 敏感，但此前尚未有将 ARFs 应用于实际数据以进行 $f_{\text{NL}}$ 约束的论证。

方法论
作者分析了 Quaia 目录，这是一个源自 Gaia DR3 和 unWISE 数据、包含约 130 万个类星体（QSOs）的全天样本，覆盖了广阔的宇宙体积。分析流程包括以下步骤：

1. 数据准备： 将 Quaia 样本分为两个红移分量（$\bar{z} \approx 0.97$ 和 $\bar{z} \approx 2.10$）。为两种密度场（数量计数）和 ARFs 生成了 HEALPix 地图（$N_{\text{side}}=128$）。ARF 场被定义为每个像素内加权平均红移的涨落，利用了一个与调查固有红移分布等效的宽高斯窗口。
2. 互相关： 本研究将 Quaia 密度图和 ARF 图与 Planck PR4 CMB 透镜收敛图相结合。这产生了四个可观测量：密度自功率谱（$C_{gg}^\ell$）、ARF 自功率谱（$C_{zz}^\ell$）以及它们与 CMB 透镜的互相关（$C_{\kappa g}^\ell$ 和 $C_{\kappa z}^\ell$）。
3. 理论建模： 使用修改版的 CAMB Boltzmann 代码（ARFCAMB）来计算理论功率谱。该代码整合了 ARF 核函数，并考虑了光度红移不确定性（$\sigma_z \approx 0.06(1+z)$）。模型包含了在普适性关系（假设 $p=1$）下由 $f_{\text{NL}}$ 引起的尺度相关偏差，并对星系偏差的红移演化（$b(z)$）进行了灵活参数化。
4. 系统误差与协方差： 为了处理观测系统误差，作者采用了针对 Quaia 模板进行线性回归的系统误差去投影法。通过对比数据偏移与模拟实现，确定了尺度截断：对于密度场应用 $\ell_{\text{min}}=14$，对于 ARF 自谱应用 $\ell_{\text{min}}=10$，以移除与统计涨落不符的低 $\ell$ 偏移。协方差矩阵通过 1,000 次相关模拟得出，并对模拟自谱进行了修正，以匹配在真实数据中观察到的低多极矩处的过剩功率。
5. 参数推断： 构建了基于测得功率谱和协方差矩阵的高斯似然函数。使用 MCMC 采样（emcee）来约束 $f_{\text{NL}}$ 和偏差参数（$b_0, \alpha$）。

关键结果
分析得出了对局部型非高斯性参数 $f_{\text{NL}}$（假设 $p=1$ 的普适性关系）的约束：

• 仅使用 ARFs： 结合 ARF 自相关和 ARF-CMB 透镜互相关，得到 $f_{\text{NL}} = -43^{+50}_{-44}$。
• 仅使用密度： 结合密度自相关和密度-CMB 透镜互相关，得到 $f_{\text{NL}} = -25^{+19}_{-17}$，这与之前的 Quaia 分析（Fabbian et al. 2026）一致。
• 联合分析： 结合所有四个可观测量（$C_{gg}^\ell, C_{\kappa g}^\ell, C_{zz}^\ell, C_{\kappa z}^\ell$），在 68% 置信度下得到 $f_{\text{NL}} = -3 \pm 14$。
• 拟合优度优化： 作者注意到，在包含 ARF-CMB 互相关（$C_{\kappa z}^\ell$）时，$\chi^2$ 出现了轻微张力。通过排除这一特定的互相关以改善模型的拟合优度，得到 $f_{\text{NL}} = -5^{+16}_{-15}$。

联合约束相比于之前仅使用密度和 CMB 透镜的 Quaia 测量结果，实现了约 25% 的改进。这一结果是迄今为止利用 LSS 二点统计量得到的第二紧的 $f_{\text{NL}}$ 约束（仅次于 DESI DR1 光谱数据），也是使用投影 2D 汇总统计量所能达到的最精确结果。

显著性与稳健性
论文声称，将 ARFs 作为一种宇宙学可观测量，为利用真实数据提升 $f_{\text{NL}}$ 约束提供了可衡量的改进。作者通过以下测试验证了结果的稳健性：

• 偏差建模： 当改变参考星系偏差演化模型（测试了 Piccirilli et al. 2024 和 Laurent et al. 2017 模型）以及允许密度场和 ARF 场具有独立的偏差参数时，结果保持稳定。
• 尺度截断： 约束结果对最小多极矩（$\ell_{\text{min}}$）截断的变化保持稳定，尽管在偏差参数上观察到了一些张力，这归因于 ARF 核函数相对于密度核函数对非线性尺度（$k$ 模式）具有更高的敏感性。
• 普适性关系： 当将 $p=1$ 的假设放宽至 $p=1.6$（某些文献中建议的 QSO 取值）时，约束变为 $f_{\text{NL}} = 2^{+20}_{-21}$，或者在模型无关的方法下为 $b_\phi f_{\text{NL}} = -63^{+79}_{-89}$。

作者总结道，尽管不同可观测量之间存在内部张力——可能由于需要改进对 ARF 非线性尺度的建模——但将 ARFs 加入 2D 分析是一个充满前景的方向。他们提倡在未来的大型巡天（如 DESI, Euclid 和 LSST）分析中纳入 ARFs，并指出 ARFs 可以作为标准 3D 聚类的稳健 2D 替代方案，特别是在红移不确定性不可忽略的光度巡天中。