[论文解读] The thesan project: public data release of radiation-hydrodynamic simulations matching reionization-era JWST observations

Enrico Garaldi, Rahul Kannan+14 more|arXiv (Cornell University)|Sep 12, 2023

Astronomy and Astrophysical ResearchPhysics and Astronomy被引用 3

一句话总结

theSAN 项目发布了涵盖 z ≥ 5.5 的全面辐射流体动力学模拟，结合 arepo-rt 辐射传输代码与 IllustrisTNG 星系形成模型及尘埃物理，用于模拟再电离时期。该研究提供高分辨率、公开可访问的数据产品——包括高时间采样率的笛卡尔输出和事后处理目录——与近期詹姆斯·韦伯太空望远镜对高红移星系的观测结果高度一致，涵盖主序列和质量-金属丰度关系等关键关系。

ABSTRACT

Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope, SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the THESAN simulation project at https://thesan-project.com. This project focuses on the $z \geq 5.5$ Universe, combining a radiation-hydrodynamics solver (AREPO-RT), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The THESAN suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope, affirming the accuracy and applicability of THESAN. The examples also serve as prototypes for how to utilise the released dataset to perform comparisons between predictions and observations.

研究动机与目标

• 为解决缺乏公开可用、高保真度的宇宙学模拟这一问题，这些模拟需在宇宙历史最初十亿年内同步模拟星系形成与宇宙再电离过程。
• 克服技术与计算障碍，提升科学界对先进辐射流体动力学模拟在高红移宇宙中的广泛可及性。
• 提供全面、公开可访问的数据集，使理论预测可与未来 JWST、罗马空间望远镜及其他新一代设施的多波段观测结果直接对比。
• 为科学界提供丰富且文档详尽的数据产品——包括高时间采样率输出和事后处理目录——并针对广泛使用的分析工具进行优化。
• 通过支持未来在体积、分辨率和物理解析方面的扩展（包括聚焦模拟与再电离历史解析模拟），为持续的、社区驱动的研究奠定基础。

提出的方法

• 项目采用 arepo-rt 代码自洽模拟宇宙体积中的辐射传输，捕捉早期光源产生的非均匀电离辐射场。
• 整合经过充分验证的 IllustrisTNG 星系形成模型，包含恒星与黑洞的反馈机制，并引入宇宙尘埃物理以提升真实性。
• 在 16 种不同配置下运行模拟，涵盖体积（100 至 1000 Mpc/h）、分辨率（高达 20 pc）及物理亚网格模型，以系统研究物理不确定性。
• 生成高时间采样率的笛卡尔输出，以精细时间分辨率追踪关键物理量的时间演化，支持对再电离过程的详细时序分析。
• 产出并公开发布大量事后处理产品，包括星系群目录、恒星形成历史及合成测光数据，可通过项目官网获取。
• 将数据与广泛使用的科学软件（如基于 Python 的 numpy、matplotlib 及自定义分析工具）集成，显著降低各类专业水平研究人员的使用门槛。

实验结果

研究问题

• RQ1 包含完整辐射传输的辐射流体动力学模拟在多大程度上能再现 z ≥ 6 时高红移星系的观测特性？
• RQ2 反馈方案、空间分辨率和模拟体积等物理模型的变化，对再电离时期星系预测属性的影响如何？
• RQ3 模拟星系群体在多大程度上能匹配主序、质量-金属丰度关系及紫外谱指数-星星质量关系等关键观测标度关系？
• RQ4 在 z ≥ 6 时，尘埃遮蔽的恒星形成比例如何随星系恒星质量变化？这对早期宇宙巡天中的观测选择效应有何启示？
• RQ5 局部再电离历史与环境密度偏移对高红移模拟中单个星系的结构与恒星形成特性有何影响？

主要发现

• theSAN 模拟成功再现了 z ≥ 6 时的关键观测标度关系，包括星系主序、质量-金属丰度关系及紫外谱指数-星星质量关系，与近期 JWST 数据高度一致。
• 模拟预测存在不可忽视的尘埃遮蔽恒星形成比例，且该比例随星系恒星质量增加而上升，表明尘埃消光可能显著影响高红移时的紫外光度观测。
• 高时间采样率的笛卡尔输出揭示了电离前沿与恒星形成活动的详细时序演化，使再电离过程可被视作动态、非均匀的演化过程，而非全局转变。
• 模拟表明再电离过程对局部物理过程（如反馈与辐射传输）极为敏感，不同模拟配置间在电离结构与星系属性上存在显著差异。
• 本次数据发布包含 16 种不同模拟运行，涵盖体积、分辨率与物理模型的差异，支持对早期星系形成物理中不确定性的系统性探索。
• theSAN 数据与标准分析工具的集成及全面的文档支持，显著降低了无高性能计算经验的研究人员的使用门槛。