宇宙学 · 宇宙学模拟
宇宙学模拟
宇宙结构的模拟，以及驱动这些模拟的数值方法。
简介
宇宙学中的许多问题无法单凭纸笔解答。引力把早期近乎均匀的宇宙中的物质聚拢起来，形成由片状结构、纤维状结构和致密团块组成的网络（见图）；要详细追踪这一过程，就需要计算机模拟（cosmological simulation）。最致密的团块里是星系团（galaxy cluster）——靠自身引力束缚在一起的天体中，质量最大的就是它们。
模拟连接着理论与观测，但前提是两者要在同等条件下比较。模拟出的星系团必须先像真实的星系团那样“观测”一遍，才能与之比较：例如把它转换成仿真的X射线图像，因为星系团中星系之间的热气体会发出X射线。
模拟内部的数值方法同样重要。模拟要分成许多小的时间步来推进，每一步都是一种近似。采用某些方法时，能量等本应守恒的量会出现微小的误差，并在数百万步中不断累积；辛方法（symplectic method）则能让这些误差始终保持在较小的水平。因此，数值方法决定了模拟结果在长时间尺度上有多可信。
为什么这很重要？像南希·格蕾丝·罗曼空间望远镜这样的巡天将精确测量宇宙结构；而在不同的宇宙学模型下宇宙应该是什么样子，要靠模拟来预测。我的工作涉及这两个方面：我的第一个项目比较了模拟星系团的仿真观测，之后我又发展了新的数值方法（见下文）。
对这一领域还不熟悉？术语表用通俗的语言解释了相关术语；带下划线的术语都链接到术语表。
我的工作
我的科研工作是从宇宙学模拟起步的：2019年，我作为访问本科生在麻省理工学院做研究。我们利用Mock-X框架，研究了如何判定星系团处于动力学弛豫状态。我们考察了九种常用的形态学指标，所用的是5.4万幅星系团合成X射线图像，来自IllustrisTNG、BAHAMAS和MACSIS模拟。结果发现，由不同指标定义的弛豫子样本彼此吻合得很差，因此挑选“弛豫”星系团会引入显著的选择效应（Cao et al. 2021）。
读博期间，我还发展了ContEvol——一族隐式数值方法，只需求解线性方程，并且近乎保辛。这些方法给出的是连续解，能处理从轨道运动到薛定谔方程的各种问题（Cao 2025；代码）。
论文
第一作者论文
- ContEvol formalism: numerical methods based on Hermite spline optimization
· DOI · arXiv:2405.05188 · ADS · 代码 - Studying galaxy cluster morphological metrics with Mock-X
· DOI · arXiv:2006.10752 · ADS · 报告与墙报：SMURF，麻省理工学院，2019年
合作论文
- Characterizing hydrostatic mass bias with Mock-X
· DOI · arXiv:2001.11508 · ADS
我对Mock-X框架的早期数据产品提供了反馈，参与了关于弛豫判据的讨论，并对稿件提出了意见。
下一步
ContCosmos：我一直希望重返宇宙学模拟领域，包括N体模拟。如今我对大尺度结构的统计量了解得深入多了，也更清楚我们能从模拟中学到什么；我计划探索能否用ContEvol来构建这类模拟。
这些工作在我未来五年计划中的位置，见研究方向。
本页术语（4）
- 宇宙学模拟（cosmological simulation）
- 用计算机从早期宇宙一直算到今天，追踪物质（有时还包括气体和恒星）在引力作用下的运动。N体模拟用大量粒子来代表物质。
- 星系团（galaxy cluster）
- 由引力束缚在一起的成百上千个星系，其间充满数千万度的气体，这些气体在X射线波段发光。星系团是引力束缚天体中质量最大的。如果一个星系团自上次并合以来已经稳定下来，就称它是弛豫的（relaxed）。
- 辛方法（symplectic method）
- 一种数值方法，用来追踪物体在引力等力作用下的运动。它的构造使每个时间步都遵守这类运动的几何结构（辛结构）；由此带来的一个结果是“相空间体积”守恒（相空间体积综合度量了位置和速度的弥散程度）。因此，即使计算很长时间，能量等守恒量的误差也始终很小，不会稳步增长。我提出的一族数值方法ContEvol近乎保辛。
- 泽尔多维奇近似（Zel'dovich approximation）
- 描述结构形成的一个简单模型：每一小团物质都沿着早期涟漪所决定的方向作直线运动，于是物质在路径交汇处堆积起来。它刻画了宇宙网最初的形成过程。