玻尔百科我们宇宙中绝大多数的物质是不可见的。这种“暗物质”塑造了宇宙，形成了星系赖以建立的引力支架，但其基本性质仍然是现代科学中最巨大的谜团之一。为了解开这个谜团，科学家们开发了一种强大的工具：宇宙学暗物质模拟。这些数字宇宙使我们能够根据观测到的宇宙结构来检验我们的宇宙学和粒子物理学理论。但是，要有效地运用如此强大的工具，我们必须首先了解它是如何构建的，其固有的局限性，以及它与几乎所有物理学领域的深刻联系。
本文将带领读者深入这些宇宙实验室的核心。在第一部分“原理与机制”中，我们将揭开在计算机内部构建宇宙的幕后过程。我们将探讨其基本概念，从创造一个有限但无边界的宇宙，到设定反映时间之初的初始条件，并审视驱动宇宙演化的复杂算法。我们还将直面“机器中的幽灵”——每个模拟研究者都必须学会驯服的数值赝象。随后的“应用与跨学科联系”部分将展示这些模拟令人难以置信的科学影响力。我们将看到它们如何超越了创造宇宙网的精美图像，成为解读星系巡天、检验暗物质本质以及指导我们在地球上最灵敏的实验中寻找这种难以捉摸的物质的不可或缺的工具。
模拟宇宙是一项相当大胆的尝试。我们试图在一个盒子里构建一个宇宙，一个我们自己宇宙的数字回响，并由相同的物理定律支配。这样的壮举怎么可能实现呢？这并非通过蛮力，而是通过一系列巧妙的原理和机制，每一项都证明了我们对宇宙宏伟设计的理解。让我们踏上创造这样一个数字世界的旅程，从它的最初时刻到其成熟、结构化的状态。
第一个、也是最突出的问题是，宇宙实际上是无限的。然而，我们的计算机是有限的。我们如何可能捕捉一个无限的广袤空间呢？我们从宇宙本身得到启示。在最大尺度上，宇宙在任何地方、任何方向上看起来都是相同的。这就是宇宙学原理，一个深刻的宇宙谦逊宣言：我们在宇宙中并不占据特殊的位置。
为了模拟这一点，我们不模拟一个有硬边界的盒子。相反，我们使用一个巧妙的技巧，称为周期性边界条件（PBC）。想象宇宙被我们模拟盒子的无限个相同副本所平铺。一个从盒子右侧离开的粒子不会撞到墙上；它会立即从左侧重新进入，就像经典街机游戏中的角色环绕屏幕一样。这种数学上的巧妙手法创造了一个没有边缘、没有中心的空间，一个有限的区域，完美地模仿了真正均匀宇宙的平移不变性。
当然，这个漂亮的技巧也有其局限性。通过将自己限制在一个边长为的盒子里，我们对任何大于盒子本身的宇宙结构波都是盲目的。我们的模拟有一个对应于波长的“基本模式”；所有更长的模式都根本不存在。这意味着我们模拟的这片宇宙，根据其构造，其平均密度将恰好等于宇宙的平均密度。然而，一个同样大小的真实区域会有其自身的轻微过密或欠密，这是由这些缺失的长波模式引起的涨落。这种“有限盒子效应”是不确定性的一个根本来源，提醒我们我们总是在观察一个无限拼图的有限碎片。[@problem-id:3494821]
随着我们的宇宙舞台搭建完毕，我们必须决定如何安排演员。我们不能简单地随机撒播暗物质粒子。我们今天看到的宇宙是从其婴儿期存在的无限小的密度涨落中成长起来的，这些结构种子的印记在宇宙微波背景中仍然可见。这些初始涨落的“配方”是一个称为功率谱的统计描述，记为。它告诉我们不同物理尺度上密度变化的幅度，从微小的涟漪到广阔的宇宙波。
但是，我们如何从一个抽象的统计配方转变为在我们的盒子中放置具体的粒子呢？在这里，我们使用优雅的Zel'dovich近似。在宇宙早期的线性区域，结构的增长很简单。粒子们还没有开始它们那复杂、旋转着汇入星系的舞蹈；它们大多只是随着哈勃膨胀而相互远离。Zel'dovich近似为每个粒子提供了“行动指令”，计算出其偏离一个完美均匀网格的微小初始位移。这些位移被精心设计，使得最终的粒子分布在统计上平均恰好具有功率谱所描述的正确属性。这是连接原始宇宙的理论描述与我们模拟的第一个快照的惊人直接的联系。
这个位移场告诉粒子从哪里开始。但它们移动的速度如何呢？对于标准的冷暗物质（CDM），其热运动被假定为可以忽略不计，这就足够了。但如果暗物质并非完全冷呢？像温暗物质（WDM）这样的候选者会拥有其从原始等离子体中退耦时残留的热速度。这不仅仅是一个小细节；这些初始的“踢动”可能产生深远的影响，因为快速移动的粒子可以从小的密度涨落中流出，抹去最小星系的种子。
为了模拟这一点，我们必须给我们的粒子一个从适当的统计分布中采样的初始速度——在这种情况下，对于费米子遗迹，是费米-狄拉克分布。通过从该分布计算预期的均方根速度，我们可以为我们的模拟粒子赋予正确的热运动量，确保我们的初始设置忠实地代表了我们希望测试的特定暗物质模型。 这突显了一个关键主题：模拟的设置本身就是一个物理假设。
场景已经布置好，粒子各就各位。现在，我们说“要有引力！”，然后观察我们的宇宙演化。但这是在一个相当奇特的环境中的引力：一个膨胀的空间。为了描述运动，我们不使用熟悉的物理坐标。相反，我们使用共动坐标，这就像在一个气球上画一个网格，然后看着网格点随着气球膨胀而分开。一个粒子的运动随后被分成两部分：总体的哈勃膨胀（网格的拉伸）和它相对于网格的“奇特”运动。
当我们在该共动参考系中从牛顿定律推导运动方程时，出现了一个引人入胜的新项：一个“哈勃阻力”。这是一个摩擦力，与粒子的奇特速度成正比，作用是减缓其相对于宇宙网格的速度。它在传统意义上不是一个“真实”的力；它是我们共动视角的结果，一个像科里奥利力一样的虚拟力，代表了随着宇宙膨胀动量的红移。完整的运动方程优美地将这种宇宙摩擦与密度涨落产生的我们所熟悉的引力结合在一起。
求解这些方程是最大的计算障碍。我们数十亿个粒子中的每一个都感受到来自其他每一个粒子的引力拖拽。直接逐个粒子地求和力将是一个的计算，一个计算上的噩梦，完成它需要比宇宙年龄还长的时间。物理学家们，作为聪明的生物，已经开发出两种主要的变通方法。
粒子-网格（PM）方法：这种方法效率极高。我们不计算单个粒子之间的力，而是首先将质量分布的“模糊图像”描绘到一个规则的网格上。一旦质量在网格上，我们就可以使用一种称为快速傅里叶变换（FFT）的数学工具，几乎瞬间解出泊松方程以得到引力势。这种方法在计算温和、长程的引力时表现出色，但就其本质而言，它在描述邻近粒子之间尖锐、细致的力时是模糊且不准确的。
树方法：这种方法采用分层的方式。为了计算作用在某个粒子上的力，它将一个遥远的粒子团簇不视为数千个个体，而是一个单一的大物体。它将所有粒子组织成一个树状数据结构，其中邻近的粒子被分组到小的“叶”节点中，而更大的群体被捆绑到“枝”节点中。这个近似，由一个决定一个群体是否足够远以便被视为一个整体的“张角”控制，是高度准确的，并且捕捉到了PM方法遗漏的关键的短程力。
在实践中，许多最先进的模拟使用混合的Tree-PM方法，利用快速的PM方法处理长程力，利用准确的树方法处理短程力，从而兼得两者的优点。
模拟是一种近似，一个好的科学家必须了解他们工具的局限性。在构建我们的袖珍宇宙时，我们被迫引入一些“非物理”的成分，以使计算易于处理。其艺术在于确保这些成分不会毒害最终结果。
其中一个成分是引力软化。在真实世界中，暗物质是一种光滑的流体。在我们的模拟中，它由离散的粒子代表。如果两个这样的粒子发生非常近的相遇，它们之间的引力将急剧飙升，将它们以荒谬的高速抛出，并使模拟陷入停顿。为了防止这种情况，我们在非常小的尺度上修改引力，使相互作用“软化”。力不再发散，而是在零分离时平滑至零。然而，这个必要的修正有其后果。它阻止了暗物质晕预测中应有的尖锐中央密度“尖峰”的形成，反而创造了人为的恒定密度“核心”。当我们分析我们模拟的暗晕时，我们必须小心区分这些数值核心和可能存在的真实物理核心，并在我们的测量中考虑软化引入的偏差。
一个相关的赝象是双体弛豫。我们大质量的模拟粒子可以通过引力相互散射，而真实暗物质那种几乎无碰撞的、细粒度的流体则不会。这种虚假的散射可以人为地“加热”系统，使暗晕的核心膨胀，并且最有害的是，扰乱较小子暗晕的轨道，可能摧毁它们。我们可以估计这种数值加热变得显著的时间尺度，即弛豫时间，。一个稳健的模拟必须被设计——通过仔细选择粒子质量和软化长度——以确保这个弛豫时间远远长于宇宙的年龄，从而保证我们的结果不被这种数字摩擦所污染。
最后，我们必须选择我们的时间步长，。在两次计算之间，我们可以跳跃多长的时间？规则很简单：在单一步长内，任何粒子的移动都不能太远。我们必须选择一个足够小的来精确追踪在最紧密束缚区域内移动最快的粒子的路径。这在模拟WDM时尤其苛刻，因为粒子开始时就有很大的热速度，迫使模拟采取极小的初始步长来跟上。
虽然我们的“盒子里的宇宙”是一个强大的工具，但它在各处的分辨率都是均匀的。如果我们只对单个星系感兴趣，但希望以极高的细节观察它呢？用如此高的精度模拟整个宇宙体积将是一种浪费。
这就是放大模拟的用武之地。我们首先运行一个大宇宙体积的低分辨率模拟，以确定一个感兴趣的暗晕将形成的区域。然后，我们只对该区域重新进行模拟，但使用高得多的分辨率。高分辨率区域中的粒子质量要小得多，我们使用更精细的网格或更小的软化长度来捕捉更多细节。这种技术，通常用自适应网格加密（AMR）实现，就像在宇宙的某一部分放置了一个计算显微镜。为了确保高低分辨率边界处的物理是一致的，我们必须遵循严格的标度律。当我们将分辨率提高一个因子时，粒子质量必须减少，软化长度必须减少，这是一个优美的结果，直接源于我们对采样质量和解析力的要求一致性。[@problem_targ_id:3503453]
模拟不仅仅是为了证实我们的标准模型；它们是我们测试新思想的主要实验室。如果暗物质不仅仅是一个引力上的旁观者呢？在自相互作用暗物质（SIDM）模型中，粒子可以像台球一样碰撞和散射。我们可以将这种物理添加到我们的模拟中。在每个时间步长的引力计算之外，我们搜索附近的粒子对。根据相互作用截面决定的概率，我们让它们发生散射。要正确地做到这一点，我们必须执行一个经典的教科书物理计算：跳入粒子对的质心系，根据期望的散射角旋转它们的相对速度矢量，然后再跳回模拟参考系。这确保了碰撞中动量和能量都完美守恒，使我们能够准确预测这种相互作用会产生的独特特征——比如在暗晕中形成真实的物理核心。
在我们的模拟运行了模拟的138亿年后，我们得到了一张包含数十亿个点的快照。我们如何从这个数字星云中提取科学信息呢？第一步是找到已经形成的引力束缚结构——承载星系的暗物质晕。
这是暗晕寻找器的工作。最常见的类型之一是球形过密度（SO）寻找器。该算法在概念上很简单：它找到一个局部密度峰值，并在其周围生长一个球体，计算所包含的质量。它继续增长球体，直到内部的平均密度下降到一个特定的阈值，例如，宇宙临界密度的200倍，。达到这个条件的半径，，和所包含的质量，，定义了这个暗晕。这个阈值并非任意的；它受到优美的球形塌缩模型的启发，这是一个简单的解析理论，表明一个均匀的球形过密度将会塌缩并在大约这个密度值上达到维里化。 通过在我们的最终粒子分布上运行这些寻找器，我们可以生成庞大的暗晕目录，记录它们的质量、大小、形状和子结构。这个目录是我们模拟的最终产品，是连接我们的“盒子里的宇宙”和我们用望远镜观测到的真实星系宇宙的桥梁。
在深入了解了驱动暗物质模拟的原理和机制之后，我们可能会倾向于将它们仅仅视为奇观——一个沉默、不可见的宇宙的宏伟数字立体模型。但这将是一个深远的错误。这些模拟不仅仅是为了展示；它们是有史以来构想出的最强大、用途最广泛的科学仪器之一。它们是计算实验室，我们可以在其中预演宇宙，对我们的理论进行压力测试，并学习如何解读我们从真实宇宙中收集到的微弱信息。它们是连接宇宙学宏伟架构与地面实验所寻求的精微信号的必要桥梁。让我们踏上一次旅程，看看这些虚拟宇宙如何延伸并触及现代物理学的几乎每一个角落。
想象一下飞越一个暗物质模拟。你不会看到均匀的迷雾，而是一个令人惊叹、错综复杂的结构网络。巨大的空旷区域，称为空洞，被纤细的壁和长而粗的物质绳索所勾勒。这些绳索，或称丝状结构，构成了宇宙的超级高速公路，将物质输送到星系团诞生的巨大十字路口。这整个结构被称为宇宙网。
模拟不仅仅是产生这幅美丽的图像；它们使我们能够对其进行量化。例如，如果我们测量一个模拟丝状结构的质量作为其长度的函数，我们可能会发现像 这样的关系。这非常有趣！一条简单的线其质量会与其长度成正比，即 。模拟的结果告诉我们，这些宇宙丝状结构不是简单的线。它们具有更复杂的、“褶皱”的结构，这是被称为分形的数学对象的标志。这里的指数 是分形维度，是对丝状结构复杂几何形状的精确测量。通过模拟，我们可以超越诗意的描述，为宇宙的纹理赋予具体的数字。
从宇宙网放大到其节点，我们找到了基本的构建模块：暗物质晕。在这里，模拟也打破了我们最简单的直觉。为了数学上的简便，我们可能倾向于将这些暗晕建模为完美的球体。然而，模拟一致地向我们表明，暗晕不是球形的。它们通常是三轴的——像稍微压扁和拉长的橄榄球。这不仅仅是一个美学细节；它对我们如何解读天文观测有实际后果。
考虑一下Tully-Fisher关系，这是天文学家用来根据旋涡星系的旋转速度估算其质量的一个经验法则。这个速度是由星系的暗物质晕的引力决定的。如果一位天文学家假设暗晕是球形的，但实际上它是一个垂直延伸的扁长形状，那么他们的质量估计将出现系统性错误。模拟通过揭示暗晕真实的、非球形的性质，迫使我们改进我们的观测模型。它们起到了关键的修正作用，提醒我们宇宙往往比我们最简单的假设要复杂得多。
模拟中的暗物质晕是宇宙的无形骨架。为了使它们真正有用，我们必须学会如何用构成我们实际用望远镜看到的星系的恒星和气体来充实它们。这种将星系“描绘”到暗物质分布上的过程是宇宙学模拟最活跃和最重要的应用之一。
关于如何做到这一点，主要有两种学派。一种方法，称为暗晕占据分布（HOD），是统计性的。它就像一个宇宙食谱：对于一个特定质量的暗晕，它提供了一个关于它应该包含多少个星系的概率分布。一个大质量暗晕几乎肯定会拥有一个大的中心星系，并且有一定概率也拥有一小撮较小的卫星星系。
另一种方法，子暗晕丰度匹配（SHAM），更具确定性。它基于一个简单但强大的假设：最重的子暗晕应该承载最亮的星系。人们将模拟中的所有子暗晕按某个属性（如它们的质量）排序，并将巡天中的所有观测星系按其光度排序，然后将它们一一匹配。模拟揭示的一个微妙之处是，子暗晕的当前质量并不是最好的使用属性。卫星子暗晕可能会被其宿主暗晕的潮汐力严重剥离其暗物质，使其大大缩小。然而，它的星系，由于更紧凑，可能大部分保持完整。事实证明，一个更好的光度代理是子暗晕在其历史上曾经拥有的峰值质量（），在它被剥离之前。
这就把我们带到了一个更深层次的模拟分析：追踪暗晕的完整生命史。通过将模拟的一个快照中的暗晕与下一个快照中的暗晕联系起来，我们构建了并合树，追踪它们的生长、合并和演化。在这里，我们同样必须巧妙。一个在一个巨大暗晕内运行的小子暗晕可能会被潮汐力剥离得如此严重，以至于其粒子数低于模拟的分辨率极限。暗晕寻找器算法“丢失”了它。但它所承载的真实星系仍然在那里，在更大的暗晕内部运行。这些被称为“孤儿”星系。为了正确地追踪它们，我们必须用解析物理来增强我们的模拟输出。当一个子暗晕从我们的目录中消失时，我们不只是从故事中删除它。我们继续将其位置作为一个孤儿来追踪，使用公认的动力学摩擦理论来模拟其轨道衰减，直到它最终与中心星系合并。这种原始数值数据和解析物理之间优美的相互作用对于创造忠于现实的模拟宇宙至关重要。
也许这些模拟最深刻的应用在于检验暗物质本身的基本性质。我们的标准模型，冷暗物质（CDM），已经取得了显著的成功，但它在小尺度上仍面临一些潜在的挑战。模拟是我们探索替代方案的主要实验室。
一个著名的难题是“核-尖点问题”。简单的CDM模拟预测，暗物质的密度应该在暗晕的正中心急剧上升，形成一个“尖峰”。然而，对一些真实星系的观测似乎表明存在一个更平坦的“核心”。这是否意味着暗物质并非完全无碰撞？如果暗物质粒子能够相互作用，即使是微弱的，它们也可能在致密的中心区域交换能量，有效地“加热”中心，并将尖峰压平为一个核心。这就是自相互作用暗物质（SIDM）的理论。SIDM的模拟证实了这种行为，产生了有核心的密度剖面。这些模拟结果随后启发了解析模型，人们可以在其中平滑地将内部的核心剖面与外部的尖峰剖面相匹配，从而提供一个可以与观测数据进行检验的具体预测。
或者，答案可能更加奇异。如果暗物质根本不是一个粒子，而是一个极其轻的量子场，表现得像一个巨大的、相干的波呢？这就是超轻暗物质（ULDM）（有时称为“模糊暗物质”）背后的思想。在这种图景中，每个粒子都有一个德布罗意波长，。对于一个典型的重暗物质粒子，这个波长是微不足道的。但如果质量极其微小（比如说， eV），波长就会变得巨大。一个简单的计算，将德布罗意关系与星系暗晕中粒子的典型速度相结合，表明对于这些超轻质量，德布罗意波长可以变得和星系本身一样大！当这种情况发生时，量子干涉效应变得主导，阻止了小尺度结构的形成，并自然地在暗晕中心形成了一个核心。模拟这一点需要全新的技术——在宇宙尺度上求解薛定谔-泊松方程——它为核-尖点问题提供了完全不同的解决方案。
拥有如此强大的力量也伴随着巨大的责任。模拟是现实的近似，而非现实本身。一个好的科学家必须了解其局限性，并能够区分物理发现和数值赝象。
以温暗物质（WDM）为例。在这个模型中，暗物质粒子比CDM中的更轻，在早期宇宙中移动得更快。这种速度使它们能够从小的密度涨落中“自由流出”，抹去了小尺度上的结构。当我们模拟一个WDM宇宙时，初始条件缺少这些小尺度波。然而，模拟本身是由离散粒子构成的。这种粒子表示引入了其自身的人为噪声（“散粒噪声”）。在宇宙网的长而密的丝状结构中，没有真正的物理能量来形成小团块，这种数值噪声可以被人为地通过引力放大。结果是一种奇异的现象：光滑的丝状结构分裂成一系列完全均匀间隔的团块，就像线上串珠一样。一个天真的观察者可能会声称发现了一种新的碎裂丝状结构的物理机制。但一个谨慎的实践者知道他们看到的是一个赝象——一个机器中的幽灵。理解和识别这样的数值陷阱是模拟艺术的一个关键部分。
在我们的模拟中生长的结构，其影响波及整个宇宙学，并指导着我们最雄心勃勃的真实世界实验。
宇宙微波背景（CMB）是大爆炸的余晖，是宇宙在38万岁时的一张婴儿照片。但来自那古老火焰的光子已经向我们传播了138亿年。在它们漫长的旅途中，它们穿过了我们的模拟一直在忙于构建的宇宙网。当一个光子落入一个正在增长的暗物质晕的引力势阱时，它获得能量（蓝移）。当它爬出时，它失去能量（红移）。如果暗晕的势阱是静态的，这两种效应会完美抵消。但暗晕在增长；它们的势阱随着时间加深。这意味着光子在进入时获得的能量比它在离开时失去的要多一点。这种微小的净能量增益，被称为Rees-Sciama效应，在CMB上创造了新的、微弱的温度点。我们今天模拟的宏伟结构在宇宙最古老的光上留下了它们微妙的、随时间变化的指纹。
最后，模拟是我们直接寻找暗物质粒子的不可或缺的指南。这场搜寻主要在两条战线上进行。在间接探测中，我们寻找暗物质湮灭的产物，如伽马射线。湮灭率取决于暗物质密度的平方（），因此我们期望最亮的信号来自最密集的暗物质团块。我们的模拟预测了整个宇宙中暗物质的“团块性”，因此对于计算这种伽马射线背景信号的预期强度和空间涨落至关重要。
在直接探测中，我们建造超灵敏的探测器，将它们置于地下深处，等待我们自己星系暗晕中的暗物质粒子撞击我们探测器原子核的微弱“叮当”声。这类实验的预期事件率取决于三件事：未知的粒子属性（其质量和相互作用截面）、探测器的属性，以及当地的天体物理环境——就在我们太阳系这里的暗物质粒子的密度和速度分布。我们从哪里获得这些信息？从宇宙学模拟中。我们可以运行一个“放大”模拟来创建一个高分辨率的类银河系星系模型，然后在一个相当于我们太阳的位置测量暗物质“风”。得到的局域密度和速度分布，然后被直接代入预测我们实验中反冲率的基本方程中。这在最大的宇宙学模拟和地球上最小、最灵敏的实验之间形成了一个惊人的直接联系。
从宇宙网的分形性质到地下实验室的反冲率，暗物质模拟远不止是数字新奇事物。它们是我们锻造对宇宙理解的熔炉，是帮助我们翻译宇宙语言的罗塞塔石，也是我们用来设计寻找自然界最深刻秘密之一的蓝图。