宇宙学模拟已经成为现代天体物理学不可或缺的研究手段。在望远镜无法触及的时空尺度上，数值模拟承担着检验理论、解释观测和预测未知的三重使命。从1941年瑞典天文学家Erik Holmberg用74个灯泡模拟星系交汇，到2026年“千衍”模拟追踪4.2万亿个暗物质粒子跨越120亿光年的演化历程，宇宙学模拟走过了八十余年的技术演进。本文将从方法原理、技术前沿与科学产出三个维度，系统梳理宇宙学模拟的现状与走向。
一、宇宙学模拟的基本逻辑：为什么需要“虚拟宇宙”
现代宇宙学建立在以大爆炸理论为基础的标准宇宙学模型之上。暗物质通过引力驱动物质聚集，塑造宇宙大尺度结构；暗能量则占据宇宙总能量约三分之二，推动宇宙加速膨胀。然而，宇宙大尺度结构的形成以及星系形成过程具有高度非线性特征，简单的解析模型无法精确描述这些过程。计算机数值模拟的核心价值，正在于能够覆盖观测无法触及的时间和空间尺度，将原初密度涨落如何在百亿年间演化为今日所见宇宙网的过程逐帧还原。
宇宙学模拟得以成立的逻辑前提是宇宙学原理：宇宙在大尺度上均匀且各向同性，没有任何一点是特殊的。这意味着模拟虽然无法精确复现银河系或本星系群，却能够产生统计意义上可比的“类比体”。研究者可以在模拟中筛选出与真实观测对象具有相似质量、环境和并合历史的暗晕，进而检验物理模型在极端条件下的表现。
在技术实现层面，宇宙学模拟将暗物质和重子物质视为无压流体，用N个虚拟粒子对相空间分布进行离散采样，在膨胀宇宙背景下求解多粒子体系的引力动力学演化。纯引力模拟（N体模拟）处理暗物质聚集的主导过程，而流体动力学模拟则在此基础上加入气体动力学、恒星形成和反馈等亚网格物理。当前主流模拟项目通常采用“大体积N体+星系形成半解析模型”或“流体动力学模拟”两条技术路线并行推进。
二、规模竞赛与技术突破：当前前沿格局
2026年是宇宙学模拟领域标志性的一年。4月，中国科学家领衔的国际团队发布“千衍”模拟首批成果，在边长约120亿光年的立方体空间内追踪4.2万亿个暗物质粒子，输出100个演化阶段的快照，首次以超高精度描绘宇宙百亿年来的结构演化全景。该模拟使用国家天文台团队自主开发的PhotoNs软件，面向国产超算平台深度优化，实现了在上万张加速卡上长时间稳定计算，消耗超过一亿CPU核心小时和千万加速卡小时，产出原始及处理数据约13PB。
仅半年后，美国阿贡国家实验室借助“极光”百亿亿次超级计算机完成了两次规模更大的模拟，每次追踪13.4万亿个暗物质粒子，覆盖286亿光年宽度的区域。这一规模跃升的目标直指宇宙暴胀物理：标准暴胀模型预测早期密度涨落服从高斯分布，而任何可测量的“原始非高斯性”偏差都可能揭示暴胀时期存在其他场或相互作用。由于这类信号极其微弱且可能被宇宙方差掩盖，探测需要覆盖最大宇宙尺度的模拟体积。
在欧洲，FLAMINGO项目走出了另一条路径。该项目首次在单一框架下实现普通物质、暗物质与暗能量复杂相互作用的全流体力学模拟，数据量达2.5PB，在兼顾大尺度宇宙网的同时呈现星系内部的气体湍流与恒星形成过程。与纯引力模拟相比，全流体力学模拟的算力代价高出一个量级以上，但其优势在于能够自洽地处理重子物理对暗物质聚集的反作用——这正是当前宇宙学最棘手的前沿问题之一。
值得注意的还有“缩放模拟”策略的成熟。对于星系形成等小尺度问题，研究者在大体积低分辨率模拟中选取感兴趣的暗晕区域，以高分辨率重新模拟其形成历史。IllustrisTNG项目中的TNG50和TNG100系列即采用这一思路，在保证大尺度环境引力场正确的前提下，将星系内部结构的解析度推进到千秒差距量级。AGORA合作项目则更进一步，用六种不同的星系形成代码模拟同一批高红移星系，系统比较各代码在恒星质量、金属丰度和紫外光度上的收敛性，揭示亚网格反馈模型对模拟结果的敏感程度。
三、从模拟到观测：验证链条如何闭合
宇宙学模拟的可信度最终取决于与真实观测的比对。2026年“千衍”团队选择了一个极端的检验对象：距离地球约40亿光年的并合星系团阿贝尔2744。该星系团由多个子星系团碰撞并合而成，内部动力学状态活跃，结构极端复杂。研究团队从模拟中筛选出相似类比体，生成类詹姆斯·韦布空间望远镜的观测图像和引力透镜质量图，基于像素统计的比较显示，模拟与观测在星系数密度、光度密度和投影质量密度分布上高度一致。这一结果的意义在于：即便在高度非线性、结构极端复杂的星系团环境中，ΛCDM标准宇宙学模型仍能给出与观测相符的描述。
在更高红移处，检验则呈现出更复杂的图景。韦布望远镜在宇宙黎明期（红移10至12）发现了比理论预期更明亮的星系，AGORA合作项目的对比研究表明，质量较大的暗晕（在红移10时质量超过5×10¹⁰太阳质量）能够在无需额外物理的条件下重现观测到的恒星质量和紫外光度，但在更高红移处模拟倾向于低估这些性质。不同代码之间的差异主要体现在金属丰度上，说明高红移星系性质对反馈模型的实现方式高度敏感。
暗能量性质的限制则暴露出模拟方法的固有局限。一项基于IllustrisTNG的局部哈勃流分析表明，虽然星系群和星系团环境中的速度—距离关系可以稳健地约束暗晕质量和哈勃常数，但由于环境方差的本征散布，当前数据尚不足以区分不同的暗能量状态方程。这提示了一个重要的认识边界：宇宙学模拟能够排除某些极端模型，但在暗能量本质这类核心问题上，模拟的判别力受限于统计精度而非计算能力。
四、未解问题与技术走向
当前宇宙学模拟面临的核心张力在于分辨率与体积的矛盾。模拟一个星系内部的气体冷却和恒星形成需要秒差距量级的空间分辨率，而为了获得宇宙大尺度结构的统计代表性，模拟盒子的边长又必须达到数百兆秒差距以上。二者之间的尺度跨度接近十个数量级，没有任何单一模拟能够同时覆盖两端。解决路径无非两条：一是继续提升算力以同时增大体积和分辨率，二是发展更智能的亚网格模型和嵌套策略，在关键区域投入高分辨率而其余部分保持粗粒度。
另一个深层问题是重子反馈模型的校准困境。当前主流流体模拟通过调节活动星系核反馈和恒星反馈的参数来匹配星系 stellar mass function 等观测量，但这些亚网格模型的形式和强度缺乏第一性原理约束。eROSITA X射线观测最新揭示，星系群尺度暗晕中的热气体比例仅为宇宙重子份额的20%至40%，远低于早期模拟预测，暗示反馈过程将重子物质驱逐到了远超出维里半径的区域。这对模拟中反馈能量的注入方式和传播机制提出了更苛刻的约束。
超算架构的变迁也在重塑模拟软件的设计。PhotoNs针对国产加速卡的移植优化、GADGET系列向异构计算的适配、以及机器学习代理模型在亚网格物理中的探索性应用，都在推动宇宙学模拟从“跑得更大”向“跑得更准”转向。在巡天观测进入大规模统计时代的背景下，模拟数据的角色正在从“理论家的玩具”转变为“观测家的必备工具”——中国空间站巡天望远镜、欧几里得望远镜和LSST等项目的科学产出，都将高度依赖配套模拟对选择效应、系统误差和统计不确定性的精细刻画。
相关问题与回答
问1：宇宙学模拟中的“粒子”是真实存在的物质粒子吗？
不是。模拟中的“粒子”是相空间分布的离散采样点，每个粒子代表一大块体积内暗物质或重子物质的平均密度和速度。例如“千衍”模拟的4.2万亿个粒子，每个粒子代表的暗物质质量远大于任何基本粒子。粒子数越多，对引力聚集过程中精细结构的解析能力就越强。
问2：纯引力模拟与流体动力学模拟有什么区别？
纯引力模拟（N体模拟）只计算暗物质粒子在引力作用下的运动，忽略气体、恒星和反馈过程。流体动力学模拟则额外求解气体的激波、冷却、恒星形成和超新星/活动星系核反馈。纯引力模拟的算力效率高，适合大体积统计研究；流体模拟物理更完整，但计算代价通常高出一到两个数量级。
问3：“千衍”模拟为什么选择阿贝尔2744作为检验对象？
阿贝尔2744是一个正在经历剧烈并合的星系团，多个子团碰撞导致动力学状态高度复杂。选择极端环境进行检验的逻辑是：如果一个模型能在最复杂的条件下与观测吻合，它在相对平静的系统中更有可能可靠。反之，若在极端条件下失效，则说明模型存在需要修正的物理缺失。
问4：宇宙学模拟能“预测”尚未观测到的天体吗？
可以，但预测是统计性的而非精确的。模拟能够告诉研究者某一类天体（如特定质量和红移范围内的高温气体晕、罕见大质量星系团）在给定宇宙学模型下的预期数量密度和性质分布。大型巡天项目常利用模拟数据来优化观测策略和筛选标准，但模拟无法预测某个特定天体的精确坐标。
问5：为什么暗能量状态方程难以用模拟来约束？
模拟可以生成大量虚拟星系群，测量其局部速度场与距离的关系，从而拟合哈勃常数和暗晕质量。但局部环境的本征差异（如附近大质量结构的引力扰动）引入了与暗能量信号量级相当的统计噪声。当前模拟研究表明，仅凭局部运动学数据，不同暗能量模型之间的差异被环境方差淹没，需要更大体积的模拟和更多独立观测探针才能突破。
问6：韦布望远镜的高红移星系观测对模拟提出了什么挑战？
韦布在红移10至12发现的星系比许多模型预测的更明亮、恒星质量更大。AGORA项目的多代码对比显示，大质量暗晕在红移10附近能够重现观测，但更高红移处模拟倾向于低估光度。可能的原因包括：尘埃衰减的处理方式不同、金属丰度对恒星光谱的影响、以及早期恒星形成历史假设的差异。这一问题尚未完全解决。
问7：中国在宇宙学模拟领域的软件自主性如何？
“千衍”模拟使用的PhotoNs软件由国家天文台团队自主开发，针对国产超算的加速卡架构进行了多年深度优化，已实现在上万张加速卡上的高效稳定运行。这标志着中国在宇宙学模拟的核心工具链上不再完全依赖GADGET、RAMSES等国外开源代码，具备了根据自身硬件条件定制模拟软件的能力。