[论文解读] BGM FASt: Besan\c{c}on Galaxy Model for Big Data. Simultaneous inference of the IMF, SFH and density in the Solar Neighbourhood

R. Mor, A. C. Robin+2 more|arXiv (Cornell University)|Sep 10, 2018

Stellar, planetary, and galactic studies参考文献 57被引用 11

一句话总结

本文提出 BGM FASt，一种基于 Besançon 星系模型的快速、可扩展框架，利用近似贝叶斯计算（ABC）实现对太阳邻域初始质量函数（IMF）、恒星形成历史（SFH）和恒星密度的同步多参数推断。通过预采样模拟与重加权方案，BGM FASt 相较于标准 BGM 实现了约 10⁴ 的加速，同时准确再现了颜色-星等图等关键观测量，得到薄盘 SFH 随时间递减，当前恒星形成速率为 1.2±0.2 M⊙/yr，总恒星质量密度为 0.051+0.002−0.005 M⊙/pc³，以及 0.5 至 1.53 M⊙ 范围内复合 IMF 倾斜指数 α₂=2.1+0.1−0.3。

ABSTRACT

We develop a new theoretical framework to generate Besan\c{c}on Galaxy Model fast approximate simulations (BGM FASt) to address fundamental questions of the Galactic structure and evolution performing multi-parameter inference. As a first application of our strategy we simultaneously infer the IMF, the star formation history and the stellar mass density in the Solar Neighbourhood. The BGM FASt strategy is based on a reweighing scheme, that uses a specific pre-sampled simulation, and on the assumption that the distribution function of the generated stars in the Galaxy can be described by an analytical expression. To validate BGM FASt we execute a set of tests. Finally, we use BGM FASt with an approximate Bayesian computation algorithm to obtain the posterior PDF of the inferred parameters, by comparing synthetic versus Tycho-2 colour-magnitude diagrams. Results: The validation shows a very good agreement between BGM FASt and the standard BGM, with BGM FASt being $\approx 10^4$ times faster. By analysing Tycho-2 data we obtain a thin disc star formation history decreasing in time and a present rate of $1.2 \pm 0.2 M_\odot/yr$. The resulting total stellar mass density in the Solar Neighbourhood is $0.051_{-0.005}^{+0.002} M_\odot/pc^3$ and the local dark matter density is $0.012 \pm 0.001 M_\odot/pc^3$. For the composite IMF we obtain a slope of $\alpha_2={2.1}_{-0.3}^{+0.1}$ in the mass range between $0.5 M_\odot$ and $1.53M_\odot$. The results of the slope at the high mass range are trustable up to $4M_\odot$ and highly depend on the choice of the extinction map (obtaining $\alpha_3={2.9}_{-0.2}^{+0.2}$ and $\alpha_3={3.7}_{-0.2}^{+0.2}$ respectively, for two different extinction maps). Systematic uncertainties are not included. Conclusions: The good performance of BGM FASt demonstrates that it is a very valuable tool to perform multi-parameter inference using Gaia data releases.

研究动机与目标

• 解决传统银河系建模在使用现代贝叶斯推断探索大参数空间时面临的计算瓶颈。
• 开发一种与 Apache Spark 和 Hadoop 等大数据平台兼容的快速、可扩展模拟框架，以处理如 Gaia 等大型巡天数据。
• 利用 Tycho-2 的观测数据，实现对初始质量函数（IMF）、SFH 和恒星质量密度等基本银河参数的同步推断。
• 在保持动力学一致性的同时，将新框架的性能与标准 BGM 模拟进行验证。
• 展示该框架在利用真实测光和天体测量数据约束本地 IMF、SFH 和密度方面的实用性。

提出的方法

• BGM FASt 框架基于标准 Besançon 星系模型（BGM Std）的预采样模拟（母模拟）构建基础。
• 通过重加权方案根据新参数值（如 IMF、SFH、密度律）调整模拟的恒星族，而无需重新运行完整模拟。
• 该方法假设银河系中恒星的分布函数可由解析表达式描述，从而实现高效的参数空间探索。
• 与近似贝叶斯计算（ABC）算法集成，用于将合成颜色-星等图（CMD）与 Tycho-2 实测数据进行比较，实现参数推断的自动化。
• 该框架在 Apache Spark 和 Hadoop 上实现，支持在 Gaia 数据分析框架（GDAF）等大数据基础设施上进行分布式高性能计算。
• 通过施加解析约束并验证与完整 BGM Std 模拟的一致性，确保局部动力学统计平衡（LDSE）得以维持，从而保持动力学一致性。

实验结果

研究问题

• RQ1 能否开发一种快速、可扩展的模拟框架，以实现利用大型巡天数据对银河参数进行高效多参数推断？
• RQ2 与标准 BGM 相比，BGM FASt 在再现关键观测量（如 CMD 和恒星分布）方面有多准确？
• RQ3 利用 Tycho-2 数据，初始质量函数（IMF）、恒星形成历史（SFH）和太阳邻域恒星质量密度的推断值是什么？
• RQ4 系统性不确定度（特别是消光图）对高质质量区 IMF 倾斜度的影响如何？
• RQ5 该框架在多大程度上可扩展以包含非参数化 SFH、金属量依赖 IMF 或随时间演化的盘结构？

主要发现

• BGM FASt 在保持模拟颜色-星等图和恒星分布高度一致的前提下，相较标准 BGM 代码实现了约 10⁴ 的加速。
• 太阳邻域推断的恒星形成历史呈现随时间递减趋势，当前恒星形成速率为 1.2±0.2 M⊙/yr。
• 太阳邻域总恒星质量密度为 0.051+0.002−0.005 M⊙/pc³，本地暗物质密度为 0.012±0.001 M⊙/pc³。
• 对于复合 IMF，0.5–1.53 M⊙ 范围内的斜率为 α₂=2.1+0.1−0.3，与 Salpeter 倾斜度一致。
• 在高质质量区（4 M⊙），IMF 倾斜度对消光图高度敏感：使用 Drimmel 图时 α₃=3.7+0.2−0.2，使用 Stilism 图时 α₃=2.9+0.2−0.2，两者均在约 3σ 水平上拒绝 Salpeter 倾斜度。
• 薄盘的径向标长受限较弱（例如，使用 Stilism 图时为 2151+421−247 pc），且 SFH 形状与消光图选择存在退化，但当前恒星形成速率保持稳健。