POLAR-II: modeling star formation history of galaxies on the 21-cm signal from Epoch of Reionization

该研究利用 Jiutian-300 模拟与 L-Galaxies 2020 模型，揭示了星系恒星形成历史（特别是恒星年龄加权年龄）的变化会显著改变再电离时期的电离与加热过程，进而影响 21 厘米信号的全局特征和功率谱。

要点

这是一篇关于宇宙早期历史的科学论文，标题为《POLAR-II：基于 21 厘米信号模拟星系恒星形成历史》。为了让你轻松理解，我们可以把这篇论文的研究内容想象成**“侦探在调查宇宙婴儿期的‘成长日记’"**。

1. 背景：宇宙在“打嗝”和“发烧”

想象一下，宇宙在诞生后的几亿年里，就像是一个刚出生的婴儿，充满了中性氢气体（一种冷冰冰、不发光的气体）。

• 再电离时期（EoR）： 随着第一批恒星和星系诞生，它们开始向宇宙中喷射强烈的紫外线和 X 射线。这就像无数个小太阳同时点亮，把原本冷冰冰、中性（不带电）的氢气“加热”并“电离”（撕碎成带电粒子）。这个过程叫“再电离”。
• 21 厘米信号： 科学家想通过一种特殊的“无线电波”（21 厘米信号）来观察这个过程。这就像是用一种特殊的“热成像仪”去扫描宇宙，看看哪里热、哪里冷、哪里被电离了。

2. 问题：以前的“成长日记”太简单了

以前的科学家在模拟这段历史时，对星系的描述有点“偷懒”。他们假设星系里的恒星形成是匀速的，就像一个人每天吃同样多的饭，每天长同样的肉。

• 现实情况： 实际上，星系的成长充满了戏剧性。它们会经历“爆发期”（像青春期突然猛长，恒星爆发式形成），也会经历“停滞期”（像生病或绝食，恒星停止形成）。这种忽快忽慢的节奏，被称为恒星形成历史（SFH）。
• 之前的误区： 以前的模型忽略了这种“忽快忽慢”，只算总账（总共造了多少恒星），却不管这些恒星是“一口气造完”还是“细水长流”。这导致模拟出来的宇宙温度分布和电离状态不够准确。

3. 新方法：POLAR-II 的“精细账本”

这篇论文的作者开发了一个叫 POLAR-II 的新工具（就像给侦探配了一台更高级的显微镜）。

• 数据来源： 他们利用了一个超级计算机模拟（Jiutian-300）和一个复杂的星系模型（L-Galaxies 2020），这些模型能详细记录每个星系在历史上每一刻的“饮食记录”（恒星形成速率）。
• 核心创新： 他们不再假设恒星是匀速形成的，而是把每个星系真实的“成长曲线”（什么时候爆发、什么时候停顿）都输入到模拟中。
• 辐射来源： 他们不仅考虑了恒星发出的紫外线（负责电离），还考虑了 X 射线双星和超新星爆发产生的 X 射线（负责加热气体）。

4. 发现：节奏不同，结果大不同

通过对比“匀速生长”和“真实节奏”两种模拟，他们发现了一些有趣的现象：

• 年龄越老，气泡越小：
  想象星系是一个吹气球的人。如果一个人（星系）在很久以前就吹了很多气（恒星形成早，年龄大），那么他现在的“吹气能力”就变弱了，产生的电离气泡（被照亮的区域）反而比那些“年轻、正在猛吹”的星系要小一点。

  • 比喻： 就像两个烧水壶，一个是一直小火慢炖（老星系），一个是刚开大火猛烧（年轻星系）。虽然总热量可能差不多，但刚开大火的那个，周围的水（气体）瞬间沸腾得更剧烈，气泡扩散得更快。

• 真实节奏 vs. 匀速节奏：
  使用真实“成长曲线”的星系，产生的电离气泡比假设“匀速生长”的星系稍微大一点、热一点。

  • 比喻： 真实的星系像是一个有节奏的鼓手，时而重击（爆发），时而轻敲（停顿）。这种节奏产生的声波（辐射）在空气中传播的方式，比一直用同样力度敲鼓（匀速）要更复杂，能“震”开更大的范围。

• 对宇宙信号的影响：
  这些微小的差异，最终会改变我们接收到的"21 厘米信号”（宇宙的热成像图）。

  • 如果忽略恒星形成的真实节奏，我们可能会误判宇宙什么时候变热、什么时候完全电离。这就像如果你不知道一个人是“暴饮暴食”还是“少食多餐”，你就很难准确预测他的体温变化。

5. 结论：细节决定成败

这篇论文告诉我们，要真正看懂宇宙早期的“热成像图”（21 厘米信号），我们必须把星系当成有血有肉、有起伏历史的个体，而不是简单的“匀速机器”。

• 未来的意义： 随着像 SKA（平方公里阵列）这样强大的新一代射电望远镜即将投入使用，我们将能捕捉到更清晰的宇宙早期信号。POLAR-II 这种更精细的模型，将帮助科学家像拼图一样，把宇宙从“黑暗时代”走向“光明时代”的完整故事拼凑得更准确。

一句话总结：
这篇论文就像是在告诉天文学家：“别再把星系当成只会匀速跑步的机器人了，它们是有节奏、有爆发力的‘舞者’。只有听懂了它们的舞步（恒星形成历史），我们才能准确描绘出宇宙婴儿期是如何从寒冷黑暗走向温暖光明的。”

技术摘要

这是一篇关于宇宙再电离时期（Epoch of Reionization, EoR）中星系恒星形成历史（SFH）对 21 厘米信号影响的学术论文。以下是该论文的详细技术总结：

1. 研究背景与问题 (Problem)

• 背景：再电离时期是宇宙从完全中性、寒冷状态向高度电离、高温状态转变的关键阶段。21 厘米中性氢谱线是探测这一时期宇宙物理状态（如电离度、气体温度）的理想探针。
• 现有挑战：目前的半数值模拟（Semi-numerical simulations）在研究 EoR 时，通常假设星系具有恒定的恒星形成率（Constant SFR）或简化的恒星形成历史。然而，实际星系由于并合、反馈等机制，其恒星形成历史（SFH）往往包含星暴（starburst）和淬灭（quenched）阶段，导致复杂的演化过程。
• 核心问题：尽管具有相同恒星质量的星系在历史总电离光子预算上可能相似，但其不同的 SFH 如何影响星际介质（IGM）的电离和加热过程，进而如何改变 21 厘米信号（全球信号和功率谱）？现有的半数值方法通常未能准确包含这种演化的 SFH 效应。

2. 方法论 (Methodology)

作者对之前的代码 POLAR 进行了更新（称为 POLAR-II），采用了一套结合数值模拟与半解析模型的混合方法：

• 暗物质模拟：使用 Jiutian-300 N-body 暗物质模拟（300 $h^{-1}$ cMpc 盒子，$6144^3$ 粒子），提供大尺度结构密度场。
• 星系形成模型：结合半解析模型 L-Galaxies 2020 (LG20)。该模型基于分子氢（H2）的恒星形成定律，包含了气体冷却、星系并合、超新星/AGN 反馈等物理过程，能够生成详细的星系目录（包括恒星质量、金属丰度、恒星形成历史等）。
• 辐射源建模：
  • 恒星源：利用 BPASS 模型计算不同年龄和金属丰度下的光谱能量分布（SED）。
  • X 射线双星 (XRB)：基于 Madau & Fragos (2017) 的标度关系，结合恒星质量和金属丰度计算高/低质量 XRB 的光度。
  • 激波加热热 ISM (hot-ISM)：假设光度与恒星形成率（SFR）呈线性关系，产生软 X 射线辐射。
• 再电离与 21 厘米信号计算：
  • 使用一维辐射传输代码 Grizzly 进行后处理。
  • 关键改进：不再假设恒定的 SFR，而是将 LG20 生成的详细 SFH（特别是恒星质量加权年龄 $\tau_{age}$）作为输入。
  • 计算每个星系周围的 1D 电离和温度剖面，并通过 FFT 技术构建 3D 电离和温度图，最终计算 21 厘米微分亮度温度（DBT, $T_{21cm}$）。
• 对比实验：运行了六组不同的 Grizzly 模拟，对比了不同的 SFH 建模方式：
  1. 固定 $\tau_{age}$ 的恒定 SFR 模型（$\tau_{age} = 0.02, 0.1, 0.2$ Gyr）。
  2. 基于 LG20 真实 SFH 的模型（Fiducial）。
  3. 基于 LG20 星系但强制 SFR 恒定的模型（Const）。

3. 主要发现与结果 (Key Results)

A. 恒星形成历史 (SFH) 的特征

• $\tau_{age}$ 的敏感性：对于相同恒星质量（$M_*$）的星系，恒星质量加权年龄 $\tau_{age}$ 是描述 SFH 的关键参数。$\tau_{age}$ 越小，意味着近期形成的恒星比例越高（星暴）；$\tau_{age}$ 越大，意味着恒星主要在早期形成。
• 质量与红移依赖：在 $M_* > 10^6 M_\odot$ 时，平均 $\tau_{age}$ 随质量增加而减小（大质量星系持续形成恒星）；随红移降低，平均 $\tau_{age}$ 增加。

B. 电离与加热效应

• 电离气泡大小：
  • 对于相同恒星质量，$\tau_{age}$ 越大，电离气泡半径越小（约小 13%，体积小 34%）。这是因为长 $\tau_{age}$ 导致更多的复合发生（recombination），抵消了部分电离光子。
  • 对于相同的 $\tau_{age}$，使用 LG20 真实 SFH 的星系产生的电离气泡比使用 恒定 SFR 的略大（约 10%），因为 LG20 模型中早期的星暴阶段提供了更强的早期电离。
• 气体温度 ($T_k$)：
  • $\tau_{age}$ 越大，X 射线加热效率越低，导致气体温度略低。
  • 使用 LG20 SFH 的模拟中，气体温度高于恒定 SFR 的模拟，因为 LG20 模型中的星暴阶段增强了 X 射线辐射。
• 拓扑结构：不同的 SFH 建模改变了 IGM 电离和加热的拓扑结构及时序。

C. 21 厘米信号特征

• 全球信号 ($\bar{T}_{21cm}$)：
  • 小 $\tau_{age}$（强星暴）导致更强的加热和更大的电离区，使 EoR 结束得更早（$\Delta z \approx 0.5$）。
  • 使用 LG20 SFH 的模拟（如 simul_Fiducial）比恒定 SFR 模拟（simul_Fiducial_const）在 $z > 10$ 时表现出更高的 $T_{21cm}$（更少的吸收或更强的发射），因为加热更强。
• 功率谱 ($\Delta^2_{21cm}$)：
  • $z > 9$（加热主导期）：加热较弱的模拟（大 $\tau_{age}$ 或恒定 SFR）表现出更高的功率谱幅度，因为局部 21 厘米涨落更大。
  • $z < 8$（电离主导期）：电离气泡的不均匀性主导信号。使用 LG20 SFH 的模拟，其功率谱峰值出现的时间比恒定 SFR 模拟略早（因为电离完成得更快）。
  • 不同 SFH 模型在功率谱演化上留下了可观测的印记，特别是在峰值位置和演化斜率上。

4. 关键贡献 (Key Contributions)

1. 模型更新：首次在半数值框架下，将基于半解析模型（L-Galaxies 2020）生成的详细、演化的恒星形成历史（SFH）引入到 EoR 的辐射传输模拟中。
2. 多源辐射耦合：不仅考虑了恒星紫外辐射，还统一处理了 X 射线双星（XRB）和激波加热 ISM 对 IGM 的加热作用，并展示了它们与 SFH 的耦合关系。
3. 量化差异：定量证明了即使总电离光子预算相同，不同的 SFH 建模（特别是 $\tau_{age}$ 的变化）也会显著改变电离气泡的大小、气体温度分布以及 21 厘米信号的时序和功率谱特征。

5. 科学意义 (Significance)

• 解释观测数据：随着 JWST 对高红移星系性质的精确测量，以及 SKA、HERA 等射电望远镜对 21 厘米信号探测能力的提升，理解 SFH 对 21 厘米信号的影响对于正确解释观测数据至关重要。
• 打破简并：研究指出，忽略 SFH 的演化可能会导致对再电离时间、星系电离效率以及 IGM 加热历史的误判。准确的 SFH 建模是区分不同物理模型（如反馈机制、星系形成效率）的关键。
• 未来预测：该研究为下一代 21 厘米实验（如 SKA）提供了更精确的理论模板，表明未来的数据分析必须考虑星系形成历史的复杂性，而不仅仅是简单的平均参数。

总结：这篇论文通过改进 POLAR 代码，揭示了星系恒星形成历史的细节（特别是 $\tau_{age}$）对再电离时期 IGM 物理状态和 21 厘米信号的显著影响，强调了在解释未来高红移观测数据时，必须采用更物理、更演化的星系形成模型。