[论文解读] Comparison of physical and observational galaxy cluster modelling

Kamran Javid, Y. C. Perrott+4 more|arXiv (Cornell University)|May 5, 2018

Galaxies: Formation, Evolution, Phenomena参考文献 20被引用 4

一句话总结

本研究利用AMI与Planck SZ数据，对比了物理模型与观测模型在54个星系团中的表现，评估了Y(r500)与θ500的估计值。物理模型得出的Y值普遍低于观测模型，其中PM与OM I的后验分布差异最大；贝叶斯证据整体倾向于OM I，尽管PM对最低红移星系团表现出中等偏好。

ABSTRACT

We present a comparison between three cluster models applied to data obtained by the Arcminute Microkelvin Imager radio interferometer system. The physical model (PM) parameterises a cluster in terms of its physical quantities to model the dark matter and baryonic components of the cluster using NFW and GNFW profiles respectively. The observational models (OM I and OM II) model only the gas content of the cluster. The two OMs vary only in the priors they use in Bayesian inference: OM I has a joint prior on angular radius $ heta$ and integrated Comptonisation $Y$, derived from simulations, while OM II uses separable priors on $ heta$ and $Y$ which are based on calculations of the physical model. For the comparison we consider a sample of $54$ clusters which are a subsample of the second Planck catalogue of Sunyaev-Zel'dovich sources. We first compare the $Y$ estimates of the three models, and find that the PM generally yields lower estimates relative to the OMs. We then compute the Earth Mover's Distance between the $ heta$ - $Y$ posterior distributions obtained from each model for each cluster, and find that the two models which are most discrepant are PM and OM I. Finally, we compare the Bayesian evidence values obtained from each model for each cluster. OM I generally provides the best fit to the data but not at a statistically significant level, according to the Jeffreys scale. The highest evidence ratio obtained is actually in favour of the PM over OM I.

研究动机与目标

• 利用太阳yaev-Zel'dovich效应数据，比较星系团的物理建模与观测建模方法。
• 评估不同建模框架下Y(r500)与θ500估计值的一致性。
• 利用贝叶斯证据及后验分布之间的Earth Mover’s Distance（EMD）评估模型性能。
• 研究红移信息对模型差异的影响，尤其在低红移区域。
• 确定在54个星系团样本中，物理模型或观测模型哪一个更符合AMI数据。

提出的方法

• 应用物理模型（PM），通过NFW形式描述暗物质（NFW）与重子物质（广义-NFW）组分，推断星系团属性。
• 实施两种观测模型（OM I与OM II），仅对气体成分建模，二者在先验结构上不同：OM I采用联合先验，OM II采用可分先验。
• 使用贝叶斯推断与马尔可夫链蒙特卡洛方法，从AMI数据中推导Y(r500)与θ500的后验分布。
• 计算Y(r500)与θ500二维后验分布之间的Earth Mover’s Distance（EMD），以量化分布差异。
• 为每个星系团计算各模型的贝叶斯证据值，并应用Jeffreys尺度评估模型偏好。
• 使用Planck PSZ2星表数据进行星系团选择与红移信息获取，利用AMI在15 GHz的观测测量SZ效应。

实验结果

研究问题

• RQ1 在星系团样本中，物理模型得出的Y(r500)估计值与观测模型相比如何？
• RQ2 通过Earth Mover’s Distance衡量，三种模型之间Y(r500)与θ500后验分布的差异程度如何？
• RQ3 基于贝叶斯证据，哪一个模型——物理模型或观测模型——对AMI数据的拟合效果最佳？
• RQ4 红移如何影响模型预测之间的差异，尤其是在低红移区域？
• RQ5 OM I与OM II中的先验设置在多大程度上影响模型表现与后验一致性？

主要发现

• 物理模型（PM）通常得出低于观测模型（OM I与OM II）的Y(r500)估计值，但在低红移区域情况相反。
• PM与OM I模型后验分布之间的Earth Mover’s Distance最大，表明二者在Y(r500)–θ500分布预测上差异最显著。
• 在样本的三分之二中，OM I与OM II对Y(r500)的估计值在彼此的一倍联合标准差范围内。
• PM与OM II后验分布之间的差异随红移增加而减小，表明红移信息可降低模型差异。
• 贝叶斯证据在50个星系团中更倾向于OM I而非PM，但其中仅14个表现出“弱”或“中等”偏好；最高证据比实际上对最低红移星系团（z = 0.0894）更支持PM（中等偏好）。
• OM I在53个星系团中优于OM II，其中39个表现出“弱”或“中等”偏好，表明OM I整体表现优于OM II，但对大多数星系团而言证据尚不明确。