跳到主要内容
Sandbox Physics

引力把小起伏放大

宇宙蛛网与壳层穿越

几乎均匀的物质,为什么会长成宇宙蛛网? 引力会把更多物质拉向稍密的区域。不同方向收缩得快慢不一,形成片、丝与结点,周围留下更空的区域。无碰撞物质还能穿过彼此;同一位置可以有几股速度不同的流。

拖动初始起伏,观察物质怎样挤成片、聚成丝。对照位移近似与三维网格引力,看穿越以后同一位置如何容纳不同速度的物质。

交互模型宇宙蛛网与壳层穿越
宇宙膨胀尺度1.091.09
密度起伏的均方根6.366.36
网格峰值与均密度之比318318
一个长波模式的功率2790 Mpc32790\,\mathrm{Mpc}^3
网格密度位移近似参考同一排物质
穿过中心的网格截面

盒子与截面使用相同密度色标,最高一档颜色会饱和。

一条长波怎样增长

P100=L3∣1Np∑je−2πiyj,x∣2P_{100}=L^3\left|\frac1{N_p}\sum_j e^{-2\pi i y_{j,x}}\right|^2

两条曲线来自同一组初始波形。虚线把位移近似延续到穿越以后。

物理教程

物质可以折叠,也可以穿过彼此

背景引力会放大物质分布里的一点不均。某个方向先挤紧,另外两个方向还展开着,就出现一张片。继续收缩,片、丝和密集结点逐渐显出来。

为什么值得理解几股物质走到同一位置后,只看位置就说不清它们怎样运动。跟着同一排粒子走,才能看见密度图藏起的经历。

先抓重点

核心问题
粒子第一次穿过彼此时,哪里变了?
一句话直觉
同一位置可以有几个速度,密度要把各股流的正质量加起来。

核心数学模型

随膨胀增长的位移

y(q,a)=q+a s(q),r=aLy\boldsymbol y(\boldsymbol q,a)=\boldsymbol q+a\,\boldsymbol s(\boldsymbol q),\qquad \boldsymbol r=aL\boldsymbol y

第一式是本物质主导宇宙里的Zel’dovich位移参考。标签始终认出同一份物质;第二式把画面的共动坐标换成物理位置。

第一次折叠

J=det⁡ ⁣(I+a ∂s∂q),ρρˉ=1Jbefore crossingJ=\det\!\left(I+a\,\frac{\partial\boldsymbol s}{\partial\boldsymbol q}\right),\qquad \frac{\rho}{\bar\rho}=\frac1J\quad\text{before crossing}

穿越之前,同一份质量挤进更小的体积,密度便升高。第一次折叠时,单流公式出现奇点。穿越后每股流都贡献正质量;行列式为负表示取向翻转。

引力跟着变化的密度走

∇y2U=32δ,dyda=a−3/2u,duda=a−1/2(−∇yU)\nabla_y^2U=\frac32\delta,\qquad\frac{d\boldsymbol y}{da}=a^{-3/2}\boldsymbol u,\qquad\frac{d\boldsymbol u}{da}=a^{-1/2}(-\nabla_yU)

网格引力在每个积分步重新计算密度和受力。这里的无量纲正则动量已包含膨胀因子;有限网格会平滑较小结构。

常见难点

穿越以后,流还在走

容易误解挤成密片以后,粒子就粘在一起了。

正确理解这里追踪无碰撞物质。各股流穿过彼此,继续受引力影响;气体激波需要另外的物理描述。

密集结点与暗物质晕

容易误解画面上每个亮结点,都已经给出了一个晕质量。

正确理解颜色表示粗网格密度。识别晕还要选择边界定义,并有足够的力与质量分辨率。

动手实验

  1. 01

    先看一张片

    选“一张片的折叠”,拖动时间滑块,越过挤紧的时刻。

    你会看到: 选中的物质线会折回去。位置—速度图把位置重叠、速度不同的流分开。
  2. 02

    改变最初那点起伏

    水平拖动金色手柄。对照新的位移参考,再看网格引力计算完成后的结构。

    你会看到: 初始起伏更大,强烈收缩就来得更早。相机可以单独转动。
  3. 03

    比较两种演化

    选“斜织的蛛网”,叠加位移参考,再看盒子下方的功率曲线。

    你会看到: 两者从同一组波出发。网格引力会回应后来的密度,位置和长波功率便逐渐分开。

来源与数值范围

位移映射依据Zel’dovich在1970年提出的近似。宇宙学变量、CIC网格方法与时间因子参照GADGET-2方法论文;这里使用固定低通网格力。

独立对照包括物质主导的线性增长、实际半格粒子响应、正质量沉积、跨越穿越时刻的周期片层参考,以及另写的自适应积分。有限网格仍用于大尺度教学。

相关实验与文章

从小起伏到宇宙蛛网 · 阅读路径

膨胀、穿越与独立检验的引力

这条阅读路径里用到本实验的文章,后面跟着它们推荐的其他实验。

  1. 文章引力怎样把物质织成宇宙网
  2. 文章用一把随宇宙伸长的尺量运动
  3. 文章同一位置,为什么会有几种速度
  4. 文章网格怎样跟着折叠的物质求引力
  5. 文章功率记住了增长,也丢掉了什么
  6. 文章把导出的初值,交给另一个积分器
  7. 文章先说清什么叫晕,再读反馈怎样改变聚类
反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

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