PAC in DESI. I. Galaxy Stellar Mass Function into the $10^{6}{\rm M}_{\odot}$ Frontier

这篇论文《PAC in DESI. I. Galaxy Stellar Mass Function into the 10^6 M⊙ Frontier》（DESI 中的 PAC 方法 I：进入 10^6 太阳质量前沿的星系恒星质量函数）详细阐述了作者团队如何利用改进的“宇宙网周围的光度学天体”（Photometric Objects Around Cosmic webs, PAC）方法，结合 DESI Y1 光谱巡天和 DECaLS 光度巡天数据，将星系恒星质量函数（GSMF）的测量范围推低至 $10^6 M_\odot$ 量级。

以下是该论文的详细技术总结：

1. 研究背景与问题 (Problem)

• 光谱巡天的局限性： 传统的星系恒星质量函数（GSMF）主要依赖光谱巡天（如 SDSS, GAMA）。然而，光谱巡天受限于深度，通常只能探测到较亮的星系（$M_* \gtrsim 10^8 M_\odot$）。在低质量端，由于选择效应和局部体积效应（如本星系群位于“本地空洞”中），光谱数据往往低估了矮星系的丰度。
• 光度巡天的挑战： 虽然光度巡天（如 DECaLS, LSST）深度更深、覆盖更广，但缺乏精确的红移信息，难以直接构建三维空间分布和物理性质。
• 现有方法的不足： 之前的 PAC 方法（Xu et al. 2022a, Paper I）虽然结合了光谱和光度数据，但依赖于红移分箱（redshift bins），且在低红移处由于共动距离变化显著，近似处理会引入误差。此外，之前的测量下限仍在 $10^8 M_\odot$ 左右。
• 核心目标： 需要一种更严谨的方法，消除红移分箱的依赖，利用 Stage-IV 巡天（DESI + DECaLS）将 GSMF 的测量下限推向 $10^6 M_\odot$ 甚至更低，以研究暗物质性质、再电离物理及矮星系形成。

2. 方法论与改进 (Methodology & Improvements)

作者对 PAC 方法进行了重大改进，使其无需红移分箱即可工作，并提高了低红移下的精度：

• 全红移范围的物理性质推断：
  • 不再假设光度样本（pop2）中的星系具有分箱内的平均红移。
  • 通过定义一组参考红移（44 个，覆盖 $0.001 < z < 1$），利用 SED 拟合代码 CIGALE 计算每个参考红移下的物理性质（如恒星质量）。
  • 利用插值法（假设 $\log M_*$ 与 $\log D_L$ 线性相关）获得任意红移下的物理性质，精度达到 $\Delta \log M_* < 0.02$。
• 新的 PAC 估计量（Estimator）：
  • 摒弃了传统的分箱平均 $\bar{n}_2 w_p$ 方法。
  • 引入了加权对计数方法，将 $\bar{S}_2$（面密度）和 $r_1^2$（距离平方）作为权重直接纳入 Landy-Szalay 估计量中。
  • 新估计量公式：$\bar{n}_2 w_p(r_p) = \frac{D_{1,w}D_{2,w} - D_{1,w}R_{2,w} - R_{1,w}D_{2,w} + R_{1,w}R_{2,w}}{R_1 R_2}$。这消除了对红移分箱的依赖，并更准确地处理了 $r_p$ 随 $r_1$ 的变化。
• 颜色切割（Color Cuts）：
  • 为了排除前景和背景噪声，根据 DESI BGS 样本的红移分布，在 $g-r$ vs $r-z$ 颜色空间定义了保守的颜色切割，剔除不相关的星系，提高信噪比。
• 恒星质量完备性极限的确定：
  • 提出了一种仅基于 DECaLS 数据本身的方法来确定质量极限：比较深场区域（$r$ 波段深度 $>23.8$）与全场的角面密度。当面密度下降至深场区域的 95% 时，定义为该红移下的质量极限。
• GSMF 推导策略：
  • 将光度样本分为蓝星和红星。
  • 假设星系偏差（Bias）主要由恒星质量和颜色决定。
  • 结合 PAC 测得的 $\bar{n}_2 w_p$ 和光谱样本（BGS）测得的 $w_p$（投影相关函数），推导 $\bar{n}_2$（数密度）。
  • 对于 BGS 中无法直接测量 $w_p$ 的极低质量星系，利用低质量蓝星 $w_p$ 几乎恒定的特性，以及红星 $w_p$ 在特定尺度（$8 h^{-1}$ Mpc）的幂律行为进行外推。

3. 关键贡献 (Key Contributions)

1. 方法学突破： 成功去除了 PAC 方法对红移分箱的依赖，显著提高了低红移（$z < 0.1$）下的测量精度，并验证了在 $M_* < 10^9 M_\odot$ 范围内误差小于 3%。
2. 数据深度挖掘： 首次利用 DESI Y1 BGS Bright 样本和 DECaLS 数据，将 GSMF 的测量下限从之前的 $10^8 M_\odot$ 大幅扩展至 $10^6 M_\odot$（蓝星）和 $10^6.3 M_\odot$（红星）。
3. 系统误差处理： 详细量化并修正了 DECaLS 中由于星系碎片化（如旋臂、恒星形成区被误判为独立星系）导致的低质量端高估问题。通过掩膜大质量星系周围区域，发现 $M_* < 10^5 M_\odot$ 时碎片化比例显著上升（最高达 40%），并据此对结果进行了修正。
4. 环境效应揭示： 对比了宇宙平均（PAC 结果）与局部体积（Local Volume, 如 GAMA, BGS Vmax）的 GSMF，发现局部体积严重低估了低质量红矮星系的丰度，揭示了低质量星系颜色对环境的强烈依赖性。

4. 主要结果 (Results)

• 测量范围：
  • 蓝星系： 完备测量至 $10^{5.3} M_\odot$，下限延伸至 $10^{3.8} M_\odot$。
  • 红星系： 完备测量至 $10^{6.3} M_\odot$，下限延伸至 $10^{5.6} M_\odot$。
• GSMF 拟合参数：
  • 蓝星系： 符合单 Schechter 函数，低质量端斜率 $\alpha_{blue} = -1.54 \pm 0.02$。
  • 红星系： 符合双 Schechter 函数，低质量端斜率 $\alpha_{red} = -2.50 \pm 0.08$。
  • 全样本： 符合三 Schechter 函数，低质量端斜率介于 $-1.85$到$-2.56$ 之间。
• 关键发现：
  • 红矮星系的主导地位： 在 $M_* < 10^{7.6} M_\odot$ 以下，红矮星系的数密度急剧上升并超过蓝星系，成为主导。这与之前的模型预测（如 Peng et al. 2010）不符，暗示红矮星系的淬灭机制可能不同。
  • 局部体积偏差： 传统的局部体积巡天（如 GAMA, BGS Vmax）由于位于本地空洞，严重低估了低质量红矮星系的丰度。经过宇宙空洞（Local Void）修正后，BGS 数据与 PAC 结果在 $10^{7.6} M_\odot$ 以上吻合良好。
  • 与高红移对比： 与 COSMOS-Web 在 $z \sim 0.5$ 的测量相比，PAC 结果在 $M_* > 10^{7.6} M_\odot$ 范围内一致，进一步支持了低质量星系丰度被低估的观点。

5. 科学意义 (Significance)

• 暗物质性质约束： 测量到的 $10^6 M_\odot$ 量级星系数密度对暗物质粒子质量（如温暗物质模型）非常敏感。观测到的陡峭上升趋势可能为限制暗物质粒子质量提供强有力的观测约束。
• 再电离与星系形成物理： 低质量星系的丰度直接关联于宇宙再电离时期的物理过程。目前的星系形成模型难以重现观测到的低质量端陡峭斜率，表明我们对矮星系形成和淬灭的物理机制（特别是环境淬灭）理解尚不完整。
• 未来巡天的潜力： 证明了 PAC 方法在 Stage-IV 巡天（如 LSST, Euclid, Roman, CSS-OS）中的巨大潜力。这些巡天比 DECaLS 深 2-3 个星等，结合 PAC 方法有望探测到宇宙中几乎所有的星系，彻底解决矮星系缺失问题（Missing Satellites Problem）并深化对宇宙结构形成的理解。

总结： 该论文通过方法学的革新和高质量数据的结合，成功将星系普查的边界推向了 $10^6 M_\odot$ 的矮星系领域，揭示了红矮星系在低质量端的意外主导性，并指出了当前星系形成模型和局部体积观测的局限性，为未来宇宙学研究开辟了新途径。