[论文解读] AMICO galaxy clusters in KiDS-DR3: the impact of estimator statistics on the luminosity-mass scaling relation

Merijn Smit, Andrej Dvornik+5 more|arXiv (Cornell University)|Sep 24, 2021

Galaxies: Formation, Evolution, Phenomena参考文献 93被引用 2

一句话总结

本研究比较了在 KiDS-DR3 的 6,925 个 AMICO 星系团中，使用 L1 回归（最小绝对偏差）与加权均值（L2）估计器测量弱引力透镜 shear 的表现。结果表明，L1 回归可降低偏差并提高统计效率，尤其在信噪比较高时的形状测量中表现更优，从而得到一条约束更紧密的光度-质量标定关系，其幂律斜率为 1.24 ± 0.08，且两种方法在质量估计上的 1σ 差异显著。

ABSTRACT

As modern-day precision cosmology aims for statistical uncertainties of the percent level or lower, it becomes increasingly important to reconsider estimator assumptions at each step of the process, and their consequences on the statistical variability of the scientific results. We compare $L^1$ regression statistics to the weighted mean, the canonical $L^2$ method based on Gaussian assumptions, for inference of the weak gravitational shear signal from a catalog of background ellipticity measurements around a sample of clusters, in many recent analyses a standard step in the process. We use the shape measurements of background sources around 6925 AMICO clusters detected in the KiDS 3rd data release. We investigate the robustness of our results and the dependence of uncertainties on the signal-to-noise ratios of the background source detections. Using a halo model approach, we derive lensing masses from the estimated excess surface density profiles. The highly significant shear signal allows us to study the scaling relation between the $r$-band cluster luminosity $L_{200}$, and the derived lensing mass $M_{200}$. We show the results of the scaling relations derived in 13 bins in $L_{200}$, with a tightly constrained power law slope of $\sim 1.24\pm 0.08$. We observe a small, but significant relative bias of a few percent in the recovered excess surface density profiles between the two regression methods, which translates to a $1\sigma$ difference in $M_{200}$. The efficiency of $L^1$ is at least that of the weighted mean, relatively increasing with higher signal-to-noise shape measurements. Our results indicate the relevance of optimizing the estimator for infering the gravitational shear from a distribution of background ellipticities. The interpretation of measured relative biases can be gauged by deeper observations, while increased computation times remain feasible.

研究动机与目标

• 评估不同统计估计器（特别是 L1 回归与加权均值）对从背景星系椭率中推断弱引力透镜 shear 的影响。
• 评估估计器选择对星系团中推导出的光度-质量标定关系的影响，特别是在现代精确宇宙学背景下的表现。
• 调查已知对非高斯噪声具有鲁棒性的 L1 回归是否在存在复杂噪声与几何效应的真实数据中优于标准的 L2 方法。
• 量化在背景源形状测量的不同信噪比范围内，不同估计器之间的相对偏差与不确定性差异。
• 通过将模拟结果（论文 I）与 KiDS-DR3 星系团的真实数据分析进行比较，验证结果的稳健性，确保在方法论与观测复杂性下的一致性。

提出的方法

• 将 L1 回归（最小绝对偏差）与加权均值（L2）估计器应用于 KiDS-DR3 中 6,925 个 AMICO 星系团周围背景星系椭率测量，以估算弱引力透镜 shear。
• 采用晕模型方法，从两种估计器获得的超额表面密度轮廓中推导出引力透镜质量（M200）。
• 将样本按 13 个光度（L200）区间分箱，以研究在宽动态范围内的光度-质量标定关系。
• 通过使用 Σ−2cr 加权噪声与径向平均的方法，对晕中心偏移与径向依赖性等系统误差进行建模。
• 比较 L1 与 L2 估计器所得的超额表面密度轮廓及推导出的 M200 值，以量化偏差与不确定性的差异。
• 利用背景源检测的信噪比评估估计器性能，衡量在不同信噪比区间下的效率增益。

实验结果

研究问题

• RQ1 在具有非高斯本征星系形状分布的真实观测数据中，L1 与 L2 估计器在恢复真实弱引力透镜 shear 信号方面表现如何？
• RQ2 估计器选择对推导出的光度-质量标定关系有何影响，特别是在 M200 的偏差与不确定性方面？
• RQ3 在存在真实噪声与几何效应的弱引力透镜分析中，L1 回归是否提供了比加权均值更高效的统计估计器？
• RQ4 L1 与 L2 估计器之间的相对偏差与不确定性如何随背景源形状测量的信噪比而变化？
• RQ5 像晕中心偏移与本征离散度等系统效应在多大程度上影响不同估计器之间结果的一致性？

主要发现

• 与加权均值相比，L1 回归在恢复的超额表面密度轮廓中表现出约 +1% 至 -3% 的相对偏差，而加权均值的偏差为 -4% 至 -4.5%，表明 L1 显著降低了绝对偏差。
• L1 估计器在效率方面表现出微小但具有统计显著性的改进，误差条缩小了 3% 至 18%，且在高信噪比下增益更明显。
• 由 L1 回归推导出的光度-质量标定关系具有高度约束的幂律斜率 1.24 ± 0.08，且在各光度区间内保持一致。
• L1 与 L2 估计器在 M200 估计上的 1σ 差异对应于数个百分点的相对差异，表明估计器选择带来的系统效应是可测量且不可忽视的。
• 结果证实了论文 I（Smit & Kuijken 2018）的发现，在真实数据中 L1 回归对非高斯本征星系形状分布与噪声更具鲁棒性。
• 尽管计算时间更长，L1 回归仍具可行性，并为弱引力透镜分析中对有偏估计器进行事后校正提供了一种更合理、更客观的替代方案。