跳到主要内容
Sandbox Physics
暗能量在变,还是引力在变

距离怎样记下宇宙的膨胀

Sandbox Physics ·

宇宙学里的每一种距离,都是膨胀速率的倒数对红移的积分。超新星亮度也好,重子声学振荡这把尺子也好,测到的都只是一个函数:每个红移处的膨胀速率,别的什么都没有。这一章先为随时间变化的暗能量写出这个函数,再把它变成距离,最后说明引力强弱为什么从头到尾不出场。

建议先读同样的膨胀,能配上不同的引力吗。需要的数学是指数、对数,以及把积分理解成许多小段相加。读完应该能推出暗能量随时间变化时平直宇宙的膨胀速率,亲手算两种距离,并看懂 Lab 里的距离图。

尺度因子、红移和膨胀速率

用尺度因子 aa 描述膨胀,规定今天 aa 等于一。任意两点的共动间隔都按 aa 的比例变大。尺度因子为 aa 时发出的光,到达时被拉长了 1/a1/a 倍,这就是红移的定义:

1+z=1a.1+z=\frac{1}{a}.

膨胀速率是单位时间里尺度因子增长的比例,H=a˙/aH=\dot a/a,今天的值记作 H0H_0,写成 H0=100 h km s−1 Mpc−1H_0=100\,h\,\mathrm{km\,s^{-1}\,Mpc^{-1}}。光在一个哈勃时间里走过的距离 c/H0c/H_0 是天然的长度单位,h=0.68h=0.68 时等于 4408.7 Mpc4408.7\,\mathrm{Mpc}。

空间平直时,弗里德曼方程把膨胀速率和总能量密度联系起来。各种密度都除以今天的临界密度,记今天的物质比例为 Ωm\Omega_{\mathrm m}。物质随体积稀释,密度按 a−3a^{-3} 下降。辐射只在比这里早得多的时候才重要,先略去。于是

E2(a)≡H2(a)H02=Ωma−3+(1−Ωm)ρDE(a)ρDE(1).E^2(a)\equiv\frac{H^2(a)}{H_0^2}=\Omega_{\mathrm m}a^{-3}+(1-\Omega_{\mathrm m})\frac{\rho_{\mathrm{DE}}(a)}{\rho_{\mathrm{DE}}(1)}.

第二项是暗能量,平直条件把它今天的份额定为 1−Ωm1-\Omega_{\mathrm m}。Weinberg 和 White,PDG 粒子物理评述 2025,§28.1

暗能量密度怎样变化

压强为 pp、密度为 ρ\rho 的流体,在体积膨胀时对外做功,自身能量随之减少。膨胀宇宙里的能量守恒给出

dln⁡ρdln⁡a=−3(1+w),w≡pρc2.\frac{\mathrm d\ln\rho}{\mathrm d\ln a}=-3(1+w),\qquad w\equiv\frac{p}{\rho c^2}.

物质的 w=0w=0,回到 ρ∝a−3\rho\propto a^{-3}。宇宙学常数的 w=−1w=-1,密度完全不变。DESI 和许多分析把变化的暗能量写成尺度因子的一条直线,w(a)=w0+wa(1−a)w(a)=w_0+w_a(1-a):w0w_0 是今天的值,waw_a 表示过去和今天差多少。DESI DR2 BAO 结果,2025

从今天 a=1a=1 往回积分到 aa:

ln⁡ρDE(a)ρDE(1)=−3∫1a1+w0+wa(1−a′)a′ da′=−3(1+w0+wa)ln⁡a−3wa(1−a).\ln\frac{\rho_{\mathrm{DE}}(a)}{\rho_{\mathrm{DE}}(1)}=-3\int_1^a\frac{1+w_0+w_a(1-a')}{a'}\,\mathrm da'=-3(1+w_0+w_a)\ln a-3w_a(1-a).

膨胀速率里的暗能量项因此是

ρDE(a)ρDE(1)=a−3(1+w0+wa) e−3wa(1−a).\frac{\rho_{\mathrm{DE}}(a)}{\rho_{\mathrm{DE}}(1)}=a^{-3(1+w_0+w_a)}\,\mathrm e^{-3w_a(1-a)}.

过去若 w>−1w>-1,往回走时它被压缩,密度比今天大。若 w<−1w<-1,密度反而随膨胀增大,这叫幻影行为。Lab 用的正是这个背景,不改的话 Ωm=0.31\Omega_{\mathrm m}=0.31。

从膨胀速率到距离

朝我们走来的光,在一小段时间 dt\mathrm dt 里走过共动距离 c dt/ac\,\mathrm dt/a。由 dz=−da/a2\mathrm dz=-\mathrm da/a^2 和 da=aH dt\mathrm da=aH\,\mathrm dt,这一小段等于 c dz/Hc\,\mathrm dz/H。从光源一路加到我们这里,得到共动距离

DM(z)=cH0∫0zdz′E(z′).D_{\mathrm M}(z)=\frac{c}{H_0}\int_0^z\frac{\mathrm dz'}{E(z')}.

平直宇宙里它也是横向共动距离:红移 zz 处共动尺寸为 rr 的东西,张开的角度是 r/DMr/D_{\mathrm M}。另一种距离管视线方向上的间隔:一小段红移 Δz\Delta z 对应共动长度 DH ΔzD_{\mathrm H}\,\Delta z,其中

DH(z)=cH(z).D_{\mathrm H}(z)=\frac{c}{H(z)}.

星系巡天在天空方向和视线方向上量重子声学尺度 rdr_{\mathrm d},得到 DM/rdD_{\mathrm M}/r_{\mathrm d} 和 DH/rdD_{\mathrm H}/r_{\mathrm d}。超新星量的是光度距离 (1+z)DM(1+z)D_{\mathrm M},差一个整体定标。这些数无一例外都由 E(z)E(z) 算出。

一个能精确核对的例子

只有物质的宇宙里,Ωm=1\Omega_{\mathrm m}=1,E=(1+z)3/2E=(1+z)^{3/2},积分有解析解:

DM(z)=2cH0(1−11+z).D_{\mathrm M}(z)=\frac{2c}{H_0}\left(1-\frac{1}{\sqrt{1+z}}\right).

z=1z=1 时等于 0.585786 c/H00.585786\,c/H_0。Lab 背后的数值积分算出的值,与它在小数点后十位一致。

例题:同一个红移,三个宇宙

取 Ωm=0.31\Omega_{\mathrm m}=0.31、h=0.68h=0.68,暗能量是宇宙学常数。z=0.5z=0.5 时尺度因子为 2/32/3,

E2=0.31×1.53+0.69=1.73625,E=1.3177.E^2=0.31\times1.5^3+0.69=1.73625,\qquad E=1.3177.

把 1/E1/E 从零数值积分到一半,得到 DM=0.4395 c/H0≃1938 MpcD_{\mathrm M}=0.4395\,c/H_0\simeq1938\,\mathrm{Mpc}。视线方向距离 DH=c/(H0E)=3346 MpcD_{\mathrm H}=c/(H_0E)=3346\,\mathrm{Mpc}。例题里规定声学尺度为 147 Mpc147\,\mathrm{Mpc},它张开 147/1938=0.0759147/1938=0.0759 弧度,天空上约 4.3 度。

改成 w0=−0.9w_0=-0.9、wa=0w_a=0。这时 z=0.5z=0.5 的暗能量密度是今天的 1.50.3=1.1291.5^{0.3}=1.129 倍,E=1.3511E=1.3511 变大,积分里每一小段都变小。距离缩到 0.4326 c/H00.4326\,c/H_0,短了 1.6%,同一把尺子看起来大 1.6%。

最后换成接近 DESI 偏好的穿越模型,w0=−0.752w_0=-0.752、wa=−0.86w_a=-0.86。今天 ww 在负一之上,到 z=0.5z=0.5 已经是 −1.039-1.039。与宇宙学常数相比,距离在 z=0.5z=0.5 短 2.3%,z=1z=1 短 1.8%,z=2z=2 短 1.0%。这种偏移只有百分之几,却铺在一整段红移上,现在的巡天正好测得出来。

引力强弱为什么不出场

再看两条距离公式,里面只有膨胀速率。物质对引力的响应,也就是后面几章的 μ\mu,不在里面;光线偏折也不在里面。

这句话可以说得很精确。对共动距离求导,dDM/dz=c/H(z)\mathrm dD_{\mathrm M}/\mathrm dz=c/H(z)。两个平直模型若在每个红移上给出相同的 DM(z)D_{\mathrm M}(z),它们的 H(z)H(z) 就相同,DHD_{\mathrm H} 和光度距离也相同。距离把膨胀历史完全定下来,也只定下膨胀历史。

修改引力理论同样会给出某段膨胀历史。不管它是什么,都可以写成上面公式里的一个等效 w(a)w(a)。hi_class 这样的宇宙学程序就用了这份自由:背景和引力响应是两项分开的输入。Zumalacárregui 等 2017,hi_class

Lab 任务

  1. 选“暗能量在变,引力照旧”,读探测壳处与宇宙学常数模型的距离之比,再把壳拖到 z=1z=1 和 z=2z=2。这个比值在 z=0.5z=0.5 约为 0.984,过了 z=1z=1 就停在 0.980 附近。说说为什么暗能量不再占主导以后,它就不怎么变了。
  2. 选“膨胀相同,引力更强”。两条增长曲线分开了,距离之比却在每个红移上都是一。用 DMD_{\mathrm M} 的公式说明原因。
  3. 拖动 waw_a,盯着光锥上的刻度。你这个宇宙的刻度会相对淡色的宇宙学常数刻度移动。找一组设置,让 z=0.5z=0.5 的距离缩短约 2%。

习题

  1. 用守恒律说明 w=−1w=-1 时密度不变,w=−1/3w=-1/3 时密度按 a−2a^{-2} 下降。
  2. 在只有物质的宇宙里,算 z=3z=3 的 DMD_{\mathrm M},分别以 c/H0c/H_0 和 h=0.68h=0.68 时的百万秒差距为单位。
  3. 不做积分,解释为什么 w0=−0.9w_0=-0.9 让到 z=0.5z=0.5 的距离比宇宙学常数情形短。
  4. 同事说:超新星样本和 BAO 样本在膨胀历史上一致,合起来就检验了广义相对论。它们真正检验的是什么?
答案

第一题,w=−1w=-1 时 dln⁡ρ/dln⁡a=−3(1+w)\mathrm d\ln\rho/\mathrm d\ln a=-3(1+w) 为零,ρ\rho 不变。w=−1/3w=-1/3 时斜率是 −2-2,所以 ρ∝a−2\rho\propto a^{-2}。

第二题,1/4=1/21/\sqrt{4}=1/2,所以 DM=2(1−1/2)=1D_{\mathrm M}=2(1-1/2)=1,单位是 c/H0c/H_0;h=0.68h=0.68 时为 4409 Mpc4409\,\mathrm{Mpc}。

第三题,w>−1w>-1 意味着暗能量过去更密,常数 ww 时多出因子 (1+z)3(1+w)(1+z)^{3(1+w)}。密度大,中间每个红移上的膨胀速率就大,距离积分里每一小段 dz/E\mathrm dz/E 就小。

第四题,它们检验的是同一段膨胀历史能否同时拟合两种距离,包括尺子和超新星亮度的定标。检验的是背景,不是引力规律。膨胀历史相同的修改引力理论两关都能过。要检验引力,得看增长或透镜,下一章就讲这两样。

来源

  1. Weinberg 和 White,PDG 粒子物理评述 2025,§28.1
  2. DESI DR2 BAO 结果,2025
  3. Zumalacárregui 等 2017,hi_class

继续阅读

这一页直接连着哪些内容

暗能量在变,还是引力在变 · 阅读路径

从距离、增长、透镜到稳定性

  1. 文章同样的膨胀,能配上不同的引力吗
  2. 当前页面距离怎样记下宇宙的膨胀
  3. 文章增长和透镜怎样感受引力
  4. 文章四个函数描述修改引力,以及它什么时候稳定
  5. 文章复算膨胀与增长网格
  6. 文章读 DESI 2024 修改引力论文
这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

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