Prospects for measuring the Doppler magnification dipole with LSST and DESI

技术摘要：利用 LSST 与 DESI 测量多普勒放大偶极矩的前景研究

问题与动机
多普勒放大效应源于星系的本征速度（peculiar velocities），这些速度导致其在红移空间中的表观位置发生偏移，从而引起其通量和角尺度的表观变化。在低红移（$z \lesssim 0.5$）时，这种由多普勒引起的收敛度（$\kappa_v$）主导了标准的引力透镜收敛度（$\kappa_g$）。这一现象在过密度周围产生了一个特征性的偶极放大模式：向观测者靠近的星系表现为退化（变小且变暗），而远离观测者的星系则表现为放大。虽然该效应的理论基础已经建立（例如 Bonvin 2008; Bacon 等人 2014），且先前的预测表明未来巡天具有高信噪比（SNR）探测的可能性，但该效应尚未被探测到。目前仍存在一个关键差距，即需要量化现实的观测不确定性，特别是关于星系尺寸测量的不确定性，这对于重建收敛场至关重要。本文旨在通过对暗能量光谱仪（DESI）与薇拉·鲁宾天文台（Vera C. Rubin Observatory）时空遗产巡天（LSST）的联合分析，探讨探测多普勒放大偶极矩的可行性。

方法论
作者提出了一种交叉相关策略，利用星系数量计数示踪子（源自 DESI 明亮星系巡天 BGS 的光谱红移）与从星系尺寸测量中重建的收敛场（源自 LSST 成像）进行关联。

1. 样本选择： 研究使用了基于 Outer Rim $N$ 体模拟的真实合成星系目录 cosmoDC2。通过星等截断（$r < 19.5$）和颜色标准选择了一个类似于 DESI BGS 的样本，重点关注“明亮”样本以确保高质量的光谱红移。
2. 本征尺寸方差： 星系尺寸的本征散射（$\sigma_{\kappa, \text{intrinsic}}$）作为收敛度重建的基础噪声底限，直接从 cosmoDC2 目录中的真实半光半径计算得出。该值按红移和表观星等进行分箱处理。
3. 图像模拟与测量误差： 为了评估测量不确定性（$\sigma_{\kappa, \text{measured}}$），作者使用 GalSim 模拟了 LSST 第一年（Y1）质量的图像。这些模拟纳入了真实的视波函数（PSF）、探测器噪声和天光背景。星系被建模为单组分 Sérsic 轮廓（固定指数 $n=1.5$）。
4. 尺寸恢复： 使用 Galight 轮廓拟合软件包对模拟图像进行分析，以恢复星系的半光半径。测量误差被量化为恢复半径与真实半径之间的均方差。
5. 预测： 通过组合本征误差和测量误差（以平方和开根号的方式）来计算总不确定性（$\sigma_{\kappa, \text{total}}$）。该不确定性被传播到偶极矩估计量的协方差矩阵中（遵循 Andrianomena 等人 2019 的方法），以预测在四个红移箱（$0.1 \le z < 0.5$）内，数量计数与收敛度之间交叉相关的信噪比（SNR）。

主要贡献

• 现实的尺寸噪声量化： 不同于以往可能假设理想条件的研究所采用的方法，本工作通过显式模拟 LSST 类的成像条件，并应用特定的轮廓拟合流程（Galight）来量化测量误差对收敛度不确定性的贡献。
• 验证了本征主导地位： 分析证实，正如预期，由于星系尺寸变化的本征散射导致的方差主导了尺寸测量误差。然而，作者指出，cosmoDC2 中发现的本征散射可能低于此前部分文献所建议的水平。
• 联合巡天预测： 本文提供了首个专门结合 DES_I 光谱数据与 LSST 成像尺寸测量来预测多普勒放大偶极矩的详细预测，并考虑了特定的重叠区域和噪声特性。

结果

• 测量精度： 模拟结果显示，在所有红移箱中，星系的测量不确定性（$\sigma_{\kappa, \text{measured}}$）都显著小于本征尺寸离散度（$\sigma_{\kappa, \text{intrinsic}}$）。因此，总不确定性由本征散射驱动。测量偏差（$b_\kappa$）被发现是微不足道的。
• 偶极振幅： 预测的偶极信号（$\xi_{\Delta\kappa}$）在所有尺度上均为负值，在小尺度（$d \sim 10\text{--}50$ Mpc/$h$）迅速下降，在 $100\text{--}125$ Mpc/$h$ 附近达到宽阔的极小值，并在较大尺度上变得不再那么负。振幅随红移增加而减小。
• 可探测性 (SNR)： 预测表明该偶极信号具有高度可探测性。每个红移箱的总 SNR 计算如下：
  • $0.1 \le z < 0.2$ 时为 $12.98$
  • $0.2 \le z < 0.3$ 时为 $12.97$
  • $0.3 \le z < 0.4$ 时为 $12.66$
  • $0.4 \le z < 0.5$ 时为 $10.41$
    每个分离尺度下的峰值 SNR 在小尺度（$d \sim 10\text{--}25$ Mpc/$h$）处达到 $\sim 8\text{--}9$。
• 方差分解： 协方差分析显示，在低红移和小尺度下，方差由尺寸噪声项（$\Xi_2$）主导。在较大尺度和较高红移下，宇宙方差（$\Xi_1$）成为主要的误差来源。散点噪声（$\Xi_3$）在所有尺度上均可以忽略不计。

意义与主张
本文主张，利用结合了 DESI 和 LSST 数据集的联合分析，多普勒放大偶极矩是一个稳健的可观测物理量，可以在 $0.1$到$0.5$之间的多个红移箱内以高显著性（$\text{SNR} \ge 10$）被探测到。作者强调，这种探测之所以可行，是因为来自类 LSST 成像的测量误差远小于星系尺寸的本征散射。

研究结论认为，这种联合分析提供了一种独立于传统红移空间畸变（RSD）测量的方法，用于探测结构增长率和本征速度。作者指出，偶极矩估计量依赖于跨不同实验（数量计数与尺寸）对两个不同场进行交叉相关，这为抵御不相关的系统误差（如 PSF 建模误差或视宁度变化）提供了内在保护。

局限性
作者明确指出了若干限制：

• 模拟依赖于 cosmoDC2 的尺寸分布，该分布基于半解析模型，可能无法完美匹配真实数据。
• 图像模拟假设星系是孤立的，并未考虑混合（blending）（即重叠源）现象，而这在深空 LSST 图像中预计会对很大比例的星系产生影响，并可能增加测量不确定性。
• 分析假设所有星系都具有简化的单组分 Sér_sic 轮廓，而现实中的星系具有复杂的结构。
• 预测假设了一个特定的重叠区域（$\sim 5,900$ deg$^2$），虽然实际情况可能更小，但未来的巡天（如 4MOST）可能会增加这一面积。

尽管存在这些局限性，作者仍断言，在所建立的模型假设下，本征尺寸散射是主要的限制因素，由此产生的 SNR 足以实现对多普勒放大偶极矩的确定性探测。