跳到主要内容
Sandbox Physics
从小起伏到宇宙蛛网

从小起伏到宇宙蛛网 · L1

引力怎样把物质织成宇宙网

引力会把更多物质拉向稍密的区域。周围的引力在各个方向并不相同,收缩便有快有慢;较空的体积之间,片状、丝状和结点状的密集结构逐渐显出来。

一开始,只需要一点不均匀

把一大片空间分成许多小格。若某几格里的物质稍多一点,周围物质受到的引力就有一点偏向。更多物质流进来,原先微弱的密度差便逐渐放大。这样的局部差别,叫密度起伏。

这个过程发生在一个整体膨胀的宇宙里。起伏刚开始增长时,两处物质可能仍在远离,只是比平均膨胀分开得慢。随着聚集加强,足够密的区域才会真正转向收缩。平均空间在变大,局地引力聚集仍然可以继续。

问题不只是“引力会不会把东西吸到一起”。还要问:哪几个方向上的物质更快地靠拢,周围哪些地方在把物质送过来。形状就藏在这些方向差别里。

对应阅读:Zel’dovich (1970), Gravitational instability: an approximate theory for large density perturbations

同一团物质,可以先扁下去

在一小团相邻物质外画一个盒子。若盒子每一处都得到完全相同的加速度,它只是一起移动;改变形状,需要盒子两端受到不同的拉动。引力在邻近位置的这种差别,称为潮汐作用。它可以沿一个方向挤压,同时让别的方向收缩得慢一些,甚至伸长。

想象盒子沿上下方向先变薄,左右和前后还展开着,物质便集中到一张片附近。若另外一个方向也挤紧,就能留下细长的结构。这里的“先”只指这一小块如何变形。整座宇宙并没有统一的施工时间表。

Zel’dovich在1970年提出一个简洁的位移近似:先用初始引力确定运动方向,再让位移随宇宙演化放大。它把这种不等速的压缩画得很清楚,也指出了地图第一次折叠的地方。

同一保存尺度因子下,自有三维片层、细丝与斜向网络PM状态
自有已验收PM,每幅40964096粒子,同在原生第116116帧,a=1.0904635200240378a=1.0904635200240378。坐标为20 Mpc20\,\mathrm{Mpc}周期盒的共动Mpc;颜色从各自CIC密度取样,共用标尺,4040以上封顶。不同受控初值展示不同几何,并非普遍的先后演化顺序。 查看完整图像

对应阅读:Zel’dovich (1970), Gravitational instability: an approximate theory for large density perturbations · Original finite-mesh calculations and independent small-wave Python replay

丝的走向,还记着更远处的联系

一张宇宙网并不是许多孤立小盒各自坍缩之后,再随意接起来。一个位置的引力来自周围很大范围的物质。邻近密集区怎样分布,会让不同地方的运动方向彼此关联。

Bond、Kofman和Pogosyan在1996年的工作强调了这层联系:后来醒目的丝状网络,在早期密度与潮汐场中已经有较宽、较弱的线索。相近的密集峰若有合适的形变方向,它们之间的物质桥会更连贯;距离变远或方向错开,联系会减弱。后来的引力演化把这些联系压得更窄、更鲜明。

这也解释了为什么只看局部的一次片状折叠,还不足以读懂整张网。片、丝和结点的相互位置,带着初始分布的空间联系。被物质流出的区域则越来越稀,形成网孔般的空洞;空洞里仍有物质。

对应阅读:Bond, Kofman & Pogosyan, How filaments are woven into the cosmic web (1996), astro-ph/9512141v1

几股物质到了一处,还可以继续穿过

在冷暗物质常用的无碰撞近似里,粒子之间的直接碰撞很少,物质流相遇时可以互相穿过。引力仍在改变它们的速度,却没有一层气体激波把这些流立即拦住。

给原来的物质层分别做上标签。开始时,一个位置只对应一股流;压缩到后来,不同标签会到达同一位置。这叫壳层穿越。名字里的“壳”不要求球形,一排平面物质层也能互相穿越。每个标签仍能跟着自己的物质走,只是从一个位置反查它来自哪一层,答案开始不止一个。

这时只画密度,会把几股流叠在一起。再加上速度,一股朝左、另一股朝右的物质就分开了。把位置和速度一起描述的空间,叫相空间。密度图里的一个结点,可以容纳仍在快速运动的不同流。

位移近似能继续画出折叠的线。要跟随穿越以后更长的运动,还得让引力回应已经改变的密度。粒子网格方法所做的,就是反复统计物质到了哪里、重新求引力,再推进下一小段运动。把重叠位置上速度不同的物质分开,我们就能接着追踪各股流之后的路。

对应阅读:Zel’dovich (1970), Gravitational instability: an approximate theory for large density perturbations · Springel (2005), The cosmological simulation code GADGET-2, astro-ph/0505010v1

练习与解释性答案

  1. 若一团物质的每一部分都受到完全相同的额外加速度,它会因此被拉成一条丝吗?

    展开解释

    这份共同加速度只会改变整团物质的移动。要改变形状,需要邻近位置的加速度有差别,让不同部分靠拢或分开。

  2. 密度图中有一个很亮的结点,怎样判断里面的物质是否已经停住?

    展开解释

    要继续查看速度。几股流可以在同一位置重叠,并沿不同方向运动;相空间图能把它们区分开。

动手检验预测

跟着同一排物质,看它怎样折叠

  1. 先猜一向和两向压缩会留下什么形状,再依次选择“一张片的折叠”和“两向挤出一条丝”。
  2. 在“斜织的蛛网”里转动盒子,看清深度;金线始终连接同一排物质标签。

来源与计算范围

  1. Zel’dovich (1970), Gravitational instability: an approximate theory for large density perturbations

    核读原文全部六页,区分原记号的物理位移、现代共动增长与正则动量。使用首穿越前的映射;原文气体激波讨论没有套到无碰撞物质上。论文保留外链。

  2. Bond, Kofman & Pogosyan, How filaments are woven into the cosmic web (1996), astro-ph/9512141v1

    核读正文、图注及原图1。初始密度与潮汐相关组织网络,非线性演化使它更鲜明。局部首焦散变平,不等于整张网必先经历普遍片层阶段。课程图使用自有PM输出。

  3. Springel (2005), The cosmological simulation code GADGET-2, astro-ph/0505010v1

    定向核读第2.1节式1至3、第3.2及3.2.1节式20至23、第4.1节式26至31,查看原PDF第3和11页。用于宇宙学动量、PM窗口与积分方法背景。教学模型另行规定截断和有限梯度,没有短程树力。

  4. Original finite-mesh calculations and independent small-wave Python replay

    自有已验收核/线性、非线性及五预设双精度计算用于课程图与校准。公开小ZIP仅含一份精确小波当前状态、独立NumPy/SciPy轨迹复算及预期结果。示例由Node执行冻结模型和规范序列化函数产生,不冒充浏览器下载。大型原始预设另行留存,包内没有Engine或原论文。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策