[论文解读] Fisher Matrix Stability

Naren Bhandari, C. Danielle Leonard+2 more|arXiv (Cornell University)|Jan 1, 2021

Fractional Differential Equations Solutions参考文献 7被引用 13

一句话总结

本文研究了弱引力透镜巡天中费雪矩阵预测的数值不稳定性，重点关注类LSST的宇宙学探针。利用CCL和CosmoSIS计算汇聚功率谱，发现费雪矩阵中微小的相对差异（最高达0.6%）由于条件数过高，可导致逆协方差矩阵出现数量级差异，尤其在同时变化多个宇宙学参数时更为显著。

ABSTRACT

Fisher forecasts are a common tool in cosmology with applications ranging from survey planning to the development of new cosmological probes. While frequently adopted, they are subject to numerical instabilities that need to be carefully investigated to ensure accurate and reproducible results. This research note discusses these challenges using the example of a weak lensing data vector and proposes procedures that can help in their solution.

研究动机与目标

• 研究类LSST巡天中弱引力透镜功率谱费雪矩阵预测的数值不稳定性。
• 评估导数计算中的数值误差对预测参数约束可靠性的影响。
• 评估使用不同宇宙学代码库（CCL与CosmoSIS）计算时费雪矩阵及其逆矩阵的稳定性。
• 提出方法论测试以确保费雪预测在高维参数空间中的数值准确性。
• 量化费雪矩阵中微小扰动在矩阵求逆过程中因高条件数而被放大的程度。

提出的方法

• 构建了一个与LSST科学需求一致的基准宇宙学模型和巡天设置，包括红移分布和天区覆盖。
• 使用CCL v2.0.1和CosmoSIS v1.6计算剪切-剪切角功率谱，作为弱引力透镜汇聚功率谱的数据向量。
• 通过改变步长计算数值导数以评估对数值精度的敏感性，从而构建费雪矩阵。
• 计算费雪矩阵的条件数，以量化矩阵求逆过程中的数值不稳定性。
• 比较CCL与CosmoSIS在所有参数对之间费雪矩阵及其逆矩阵（协方差矩阵）的相对差异。
• 评估DETF指标（FoM）以评估尽管存在数值差异，预测约束的稳定性。

实验结果

研究问题

• RQ1 导数计算中的数值不稳定性如何影响弱引力透镜巡天中费雪矩阵预测的可靠性？
• RQ2 费雪矩阵元素中微小的相对差异在多大程度上传播为协方差矩阵逆矩阵中的巨大差异？
• RQ3 费雪矩阵的条件数与参数预测精度的数值不稳定性之间有何关联？
• RQ4 重新参数化宇宙学参数（例如使用Ω_c h²而非Ω_c）对矩阵条件数和预测稳定性有何影响？
• RQ5 在计算类LSST巡天的费雪矩阵时，能否在不同宇宙学代码库（CCL与CosmoSIS）之间确保数值稳定性？

主要发现

• CCL与CosmoSIS数据向量之间的信噪比相对差异仅为-0.04%，表明两者在数据向量上高度一致。
• 尽管费雪矩阵元素的相对差异很小（最大约0.6%），但由此产生的逆协方差矩阵却出现数量级差异，原因在于条件数过高。
• 当使用Ω_c h²与Ω_b h²作为参数时，费雪矩阵的条件数达到3.33×10⁷，表明对数值误差极度敏感。
• CCL的条件数为1.07×10⁷，高于CosmoSIS的8.08×10⁶，但两者均足够大，足以放大微小的数值误差。
• 尽管FoM值相近（CCL为25.24，CosmoSIS为29.50），但由于条件数过高，其底层费雪矩阵仍存在数值不稳定性。
• 本研究证实，高条件数可放大费雪矩阵中的微小数值误差，导致参数协方差预测不可靠，与Blanchard等（2020）先前估算一致：为实现10%的稳定性，每元素误差需控制在0.01%以内。