[论文解读] Simulating Cosmic Reionization

Ilian T. Iliev, Paul R. Shapiro+3 more|ArXiv.org|Jun 17, 2008

Galaxies: Formation, Evolution, Phenomena参考文献 21被引用 5

一句话总结

本论文展示了迄今最大规模且最详细的宇宙再电离模拟，采用 $3072^3$ 粒子的N体模拟与一种快速射线追踪辐射转移代码（C²-Ray），以模拟矮星系的电离辐射。主要成果为对电离结构拓扑与晕质量函数的自洽、高分辨率模拟，与Sheth-Tormen质量函数高度一致，并揭示了电离泡在星际介质中复杂且动态演化的结构。

ABSTRACT

The Cosmic Dark Ages and the Epoch of Reionization constitute a crucial missing link in our understanding of the evolution of the intergalactic medium and the formation and evolution of galaxies. Due to the complex nature of this global process it is best studied through large-scale numerical simulations. This presents considerable computational challenges. The dominant contributors of ionizing radiation were dwarf galaxies. These tiny galaxies must be resolved in very large cosmological volumes in order to derive their clustering properties and the corresponding observational signatures correctly, which makes this one of the most challenging problems of numerical cosmology. We have recently performed the largest and most detailed simulations of the formation of early cosmological large-scale structures and their radiative feedback leading to cosmic reionization. This was achieved by running extremely large (up to 29 billion-particle) N-body simulations of the formation of the Cosmic Web, with enough particles and sufficient force resolution to resolve all the galactic halos with total masses larger than 10^8 Solar masses in computational volumes of up to (163 Mpc)^3. These results were then post-processed by propagating the ionizing radiation from all sources by using fast and accurate ray-tracing radiative transfer method. Both of our codes are parallelized using a combination of MPI and OpenMP and to this date have been run efficiently on up to 2048 cores (N-body) and up to 10000 cores (radiative transfer) on the newly-deployed Sun Constellation Linux Cluster at the Texas Advanced Computing Center. In this paper we describe our codes, parallelization strategies, scaling and some preliminary scientific results. (abridged)

研究动机与目标

• 模拟宇宙再电离期间早期大尺度结构的形成及其辐射反馈。
• 在高达 $114/h\,\text{Mpc}^3$ 的宇宙体积中，解析所有质量超过1亿倍太阳质量的星系晕。
• 利用快速并行化的辐射转移方法，精确模拟数百万个源的电离辐射传播。
• 评估源的聚类特性与数密度对再电离可观测特征的影响。
• 开发并扩展先进的N体与辐射转移代码，以支持百亿亿级（exascale）宇宙学模拟。

提出的方法

• 使用CubeP³M N体代码，通过 $1728^3$ 至 $3072^3$ 粒子，在 $64/h\,\text{Mpc}^3$ 至 $114/h\,\text{Mpc}^3$ 的体积中模拟暗物质结构的形成。
• 采用快速并行化的射线追踪方法（C²-Ray），在三维空间中计算来自最多数百万个源的电离前沿。
• 将N体模拟与后处理辐射转移结合，以模拟星际介质中非平衡电离化学过程。
• 使用MPI与OpenMP并行化技术，将CubeP³M扩展至2048个核心，C²-Ray扩展至10,000个核心，在Ranger超级计算机上运行。
• 以Sheth-Tormen质量函数作为基准，验证模拟所得的晕质量函数。
• 实现基于网格的晕探测算法，用于识别含40至50个粒子或以上的晕，并计划发展对更小晕的网格无关方法。

实验结果

研究问题

• RQ1 质量超过1亿倍太阳质量的矮星系的聚类特性与数密度在多大程度上影响再电离期间电离泡的拓扑结构？
• RQ2 模拟体积与宇宙方差在多大程度上影响稀有、大质量晕的形成红移与丰度？
• RQ3 与Press-Schechter模型相比，Sheth-Tormen质量函数在大规模再电离模拟中对晕质量函数的描述有多准确？
• RQ4 力的分辨率与软化长度对小晕动力学及其电离反馈的影响如何？
• RQ5 如何高效耦合N体与辐射转移代码，以模拟反馈对形成结构中气体的影响？

主要发现

• 模拟在 $z=6$ 时的晕质量函数与解析预测的Sheth-Tormen函数高度一致，且显著优于Press-Schechter模型。
• 在两个最大规模模拟中，坍缩分数与活跃源数量的差异在百分之几以内，尽管模拟体积不同。
• 在 $163\,\text{Mpc}$ 箱体中，最早晕形成于 $z \sim 30$；在 $91\,\text{Mpc}$ 箱体中，形成于 $z \sim 27$，差异归因于宇宙方差与箱体尺寸效应。
• 电离结构拓扑呈现出复杂且动态演化的电离泡与中性区域网络，表明再电离过程高度不均匀。
• 模拟实现了 $3072^3$ 粒子的分辨率（290亿粒子），仅使用Ranger超级计算机约3%的资源，表明仍有潜力开展更大规模模拟。
• 作者预测，基于现有硬件，实现高达1650亿粒子的模拟是可行的，可实现可观测宇宙尺度（$h^{-1}\text{Gpc}^3$）的全体积模拟。