[论文解读] The Bayesian Global Sky Model (B-GSM): Validation of a Data Driven Bayesian Simultaneous Component Separation and Calibration Algorithm for EoR Foreground Modelling

Gregory Carter, Will Handley+2 more|arXiv (Cornell University)|Jan 2, 2025

Geophysics and Gravity Measurements被引用 3

一句话总结

B-GSM 引入一个数据驱动的贝叶斯框架，在合并建模低频前景与标定扩展天空地图的同时，使用嵌套采样选择分量数量及其谱，且在合成数据上得到验证。

ABSTRACT

We introduce the Bayesian Global Sky Model (B-GSM), a novel data-driven Bayesian approach to modelling radio foregrounds at frequencies <400~MHz. B-GSM aims to address the limitations of previous models by incorporating robust error quantification and calibration. Using nested sampling, we compute Bayesian evidence and posterior distributions for the spectral behaviour and spatial amplitudes of diffuse emission components. Bayesian model comparison is used to determine the optimal number of emission components and their spectral parametrisation. Posterior sky predictions are conditioned on both diffuse emission and absolute temperature datasets, enabling simultaneous component separation and calibration. B-GSM is validated against a synthetic dataset designed to mimic the partial sky coverage, thermal noise, and calibration uncertainties present in real observations of the diffuse sky at low frequencies. B-GSM correctly identifies a model parametrisation with two emission components featuring curved power-law spectra. The posterior sky predictions agree with the true synthetic sky within statistical uncertainty. We find that the root-mean-square (RMS) residuals between the true and posterior predictions for the sky temperature as a function of LST are significantly reduced, when compared to the uncalibrated dataset. This indicates that B-GSM is able to correctly calibrate its posterior sky prediction to the independent absolute temperature dataset. We find that while the spectral parameters and component amplitudes exhibit some sensitivity to prior assumptions, the posterior sky predictions remain robust across a selection of different priors. This is the first of two papers, and is focused on validation of B-GSMs Bayesian framework, the second paper will present results of deployment on real data and introduce the low-frequency sky model which will be available for public download.

研究动机与目标

• 开发一个面向低频天空前景（<400 MHz）的数据驱动贝叶斯模型，提供鲁棒的误差量化与标定。
• 通过对扩展天空 surveys 和绝对温度数据进行条件化，实现分量的同时分离与标定。
• 使用贝叶斯证据来确定发射分量的最优数量及其光谱参量化。
• 提供经过验证的方法学，在存在部分天空覆盖和噪声时实现准确的天空预测与标定。

提出的方法

• 将观测天空表示为 k 个分量的和，空间图 M_c(Ω) 和谱 S^c(v)；
• 为扩展地图 D、标定参数 a_v, b_v 与绝对温度数据 E 构建联合似然 P(E,D|a,b,M,S)；
• 推导对分量图 M 的解析边际化以获得边际似然；对 M 采用共轭高斯先验 P(M|S)；
• 使用嵌套采样计算贝叶斯证据和 a、b 与谱 S 的后验分布，选择分量数量及其参数化；
• 采用分块矩阵形式以高效处理扩展前景和绝对温度似然，并讨论一种实际的近似边际化以降低计算成本；
• 引入对分量振幅的条件先验 c0(S) = (S^T C_Sky^{-1} S)^{-1}，以及具有幂律频率依赖的天空协方差先验 C_Sky。

实验结果

研究问题

• RQ1 B-GSM 是否能够从低频数据中恢复出正确数量及谱形的扩展发射分量？
• RQ2 对绝对温度测量进行条件化是否能够实现扩展地图的同时标定与可靠的分量分离？
• RQ3 后验天空预测对先验假设的变动是否鲁棒？
• RQ4 在现实数据问题（部分天空覆盖、热噪声和标定不确定性）存在时，B-GSM 的表现如何？

主要发现

• B-GSM 识别出一个包含两种发射分量且谱为曲线幂律的模型。
• 贝叶斯证据强烈拒绝错误模型，表明模型选择有效。
• 后验天空预测与合成天空的真实值在统计不确定性范围内一致。
• 在以绝对温度数据进行标定后，真实天空与后验预测之间的 RMS 残差显著降低，相较于未标定情况。
• 谱参数与分量振幅对先验有一定敏感性，但在不同先验下后验天空预测仍然稳健。