[论文解读] Redshifts for galaxies in radio continuum surveys from Bayesian model fitting of HI 21-cm lines

I. Harrison, Michelle Lochner+1 more|arXiv (Cornell University)|Apr 26, 2017

Radio Astronomy Observations and Technology参考文献 1被引用 4

一句话总结

本文提出一种贝叶斯模型拟合技术，利用H I 21-cm线发射线在射电连续谱巡天中推导星系的光谱红移，将连续源位置作为先验信息。该方法在SKA1-MID模拟中对低红移星系实现了高精度红移测量（Δz/z < 10⁻⁵），覆盖率达10.14%，优于传统信噪比（SNR）基选择方法，可提供完整的后验概率分布和检测显著性。

ABSTRACT

We introduce a new Bayesian HI spectral line fitting technique capable of obtaining spectroscopic redshifts for millions of galaxies in radio surveys with the Square Kilometere Array (SKA). This technique is especially well-suited to the low signal-to-noise regime that the redshifted 21-cm HI emission line is expected to be observed in, especially with SKA Phase 1, allowing for robust source detection. After selecting a set of continuum objects relevant to large, cosmological-scale surveys with the first phase of the SKA dish array (SKA1-MID), we simulate data corresponding to their HI line emission as observed by the same telescope. We then use the MultiNest nested sampling code to find the best-fitting parametrised line profile, providing us with a full joint posterior probability distribution for the galaxy properties, including redshift. This provides high quality redshifts, with redshift errors $Δz / z &lt;10^{-5}$, from radio data alone for some 1.8 million galaxies in a representative 5000 square degree survey with the SKA1-MID instrument with up-to-date sensitivity profiles. Interestingly, we find that the SNR definition commonly used in forecast papers does not correlate well with the actual detectability of an HI line using our method. We further detail how our method could be improved with per-object priors and how it may be also used to give robust constraints on other observables such as the HI mass function. We also make our line fitting code publicly available for application to other data sets.

研究动机与目标

• 开发一种稳健的方法，从低信噪比射电巡天中的H I 21-cm线发射中获取光谱红移，尤其针对平方公里阵列（SKA）。
• 通过利用连续源位置作为先验信息，克服特征不明显射电谱中红移与流量之间的退化问题。
• 为红移及其他参数提供完整的后验概率分布，支持在宇宙学分析中的一致使用。
• 通过整合完整的数据级建模和贝叶斯证据以提高检测显著性，超越传统SNR基源选择方法。
• 仅基于射电数据实现红移估计，减少对多波段交叉匹配和光度红移的依赖。

提出的方法

• 对一维H I 21-cm线谱拟合一个六参数模型，包括线轮廓振幅、中心频率、线宽及基线参数。
• 使用MultiNest嵌套采样算法进行贝叶斯推断，以计算所有模型参数的联合后验分布。
• 以信号模型与噪声模型之间的贝叶斯证据比（ln(B)）作为检测标准，其中ln(B) > 0表示具有信心的检测。
• 利用SKA1-MID巡天中的连续源位置作为空间先验，将红移估计问题限制在单一视线上。
• 考虑仪器灵敏度响应函数，并通过S3-HI模拟和SKA1-MID灵敏度曲线模拟真实巡天条件。
• 该方法支持完整的不确定性量化，并可扩展用于从线轮廓形状中估计衍生量，如H I质量。

实验结果

研究问题

• RQ1 贝叶斯模型拟合H I 21-cm线是否能为低信噪比射电巡天中的星系提供可靠的红移，特别是在SKA1-MID条件下？
• RQ2 与传统SNR基源选择相比，该贝叶斯方法在源数和红移分布方面的性能如何？
• RQ3 在不依赖光学或近红外交叉匹配的情况下，能从射电数据中恢复多大程度的红移？
• RQ4 贝叶斯证据比（ln(B)）与传统SNR定义相比，与实际可检测性的相关性如何？
• RQ5 该方法能否扩展以对H I质量函数等额外可观测量提供稳健约束？

主要发现

• 该贝叶斯方法实现了高精度红移恢复，使模拟巡天中相当一部分星系的Δz/z < 10⁻⁵。
• 在SKA1-MID Band 2巡天中，高达10.14%的低红移星系被成功赋予高置信度红移（ln(B) > 0），同时Band 1中0.18%的高红移星系也被恢复。
• 检测源的数量和红红移分布与传统SNR > 5标准选择的结果相当，但实际源集存在显著差异。
• 常用的SNR定义与基于贝叶斯方法的实际可检测性相关性较差，表明标准预测实践存在局限性。
• 该方法为每个源提供了完整的后验概率分布P(z)和检测显著性B，支持在宇宙学参数估计中的一致集成。
• 该技术已在http://github.com/MichelleLochner/radio-z公开发布，可应用于其他数据集或模拟以进行方法比较。