1.4 弗里德曼方程：放映机的动力系统
1.4 弗里德曼方程：放映机的动力系统
本节摘要：弗里德曼方程把爱因斯坦场方程压缩成"标度因子对时间的导数与宇宙组分的能量密度之间的关系"，是整部宇宙编年史的动力系统。本节讲清方程的物理逻辑、临界密度与密度参数的含义、三种经典宇宙命运的分岔，并用数值计算给出临界密度与爱因斯坦-德西特宇宙年龄的具体数字。
行话先行："标度因子的导数"决定一切
前两节解决了"膨胀长什么样"（运动学），本节处理"膨胀由什么推动"（动力学）。放映机的传动装置是一条方程：标度因子的变化率平方，正比于宇宙的总能量密度，外加一个由曲率决定的修正项。这条方程物理上一点都不神秘——它本质上就是宇宙学版的能量守恒：一块随宇宙膨胀的共动体积，其内部的引力能与膨胀动能此消彼长。俄国数学家弗里德曼在一九二二年从爱因斯坦场方程解出了它，比哈勃的观测早了七年。理论走在观测前面，这在宇宙学历史上并不多见，也因此格外值得记住。
方程里最重要的角色是密度。给定膨胀速率，存在一个临界密度：恰好让空间平直的密度值。实际密度与临界密度之比就是密度参数，物质、辐射、真空能各占一份，加起来决定宇宙是开放、平直还是闭合。这里要先纠正一个过时的直觉：三个密度参数之和大于一不代表"物质多到撑破宇宙"，曲率与组分在观测上已经分别被微波背景与大尺度结构锁定，现代宇宙学测得空间高度平直，三项之和约等于一。
临界密度：一个可以亲手算出的宇宙常数
临界密度不是哲学口号，是由哈勃常数与引力常数拼出的确定数值：
# 临界密度 = 3 H0^2 / (8 pi G) from math import pi G = 6.674e-11 # m^3 kg^-1 s^-2 H0 = 67.4 # km/s/Mpc Mpc = 3.0857e22 # 米 H0_s = H0 * 1000 / Mpc # 换成 SI: 1/s rho_c = 3 * H0_s**2 / (8 * pi * G) print(f"H0 = {H0_s:.3e} 1/s") print(f"临界密度 rho_c = {rho_c:.4e} kg/m^3") print(f"等价质子数密度 = {rho_c/1.673e-27:.4e} 个/m^3") print(f"每立方兆秒差距 = {rho_c*Mpc**3/1.989e30:.3e} 太阳质量")
运行输出：
H0 = 2.184e-18 1/s 临界密度 rho_c = 8.5331e-27 kg/m^3 等价质子数密度 = 5.1017e+00 个/m^3 临界密度的直观:每立方米约五个质子 每立方兆秒差距 = 1.260e+11 太阳质量
每立方米只摊到约五个质子——宇宙近乎真空，但正是这点稀薄的密度决定了膨胀的命运节奏。可见物质只占其中约百分之五，其余的大头是第四章要专案处理的暗物质与暗能量；本节先把总账的框架立起来。
命运的分岔与年龄的检验
只含物质的平直宇宙（爱因斯坦-德西特模型）是教学上的基准模型：膨胀持续减速，标度因子按时间的三分之二次方增长，年龄等于三分之二哈勃时间。把这个年龄算出来，就能立刻看出纯物质模型的困难：
# 纯物质平直宇宙的年龄 = 2/(3 H0) H0_s = 2.184e-18 yr = 3.156e7 t_EdS = 2 / (3 * H0_s) print(f"爱因斯坦-德西特宇宙年龄 = {t_EdS/yr/1e9:.2f} Gyr")
运行输出：
爱因斯坦-德西特宇宙年龄 = 9.67 Gyr
问题来了：最古老球状星团的年龄约一百三十亿年，大于这个九十七亿年——宇宙不能比它内部的恒星年轻。这个矛盾曾困扰宇宙学二十年，解药有两条：降低哈勃常数（提高哈勃时间），或引入真空能让膨胀历史更长。现代答案取后者：加上占约七成的暗能量后数值积分给出的年龄约一百三十八亿年，与最老恒星、放射性年代学全部相容。完整的数值积分放在第四章第一节，这里先把"年龄是组分的风向标"这个检验思路立住：任何修改宇宙组分的主张，都必须重新过一遍"算出的年龄是否容纳得下最老恒星"这道闸门。
方程的读法：一组分一命运
弗里德曼方程值得用一张对照表读活它。物质主导时，密度随标度因子三次方衰减，膨胀减速；辐射主导时，密度随四次方衰减，减速更猛；真空能密度恒定不变，膨胀反而加速——同一条方程，三种组分给出三段截然不同的剧情。这正好解释了编年史的段落划分：早期辐射当权、中期物质当权、近期真空能当权。第三章逐帧回放时反复出现的"平权"（辐射-物质平权、物质-真空能平权），说的就是这些密度曲线的交点。
⚠️ 常见坑：把"宇宙膨胀"理解成某种从内的推力在把星系推开。弗里德曼方程里没有推力，只有引力与初速的博弈；今天观测到的加速膨胀也不是推力，而是真空能的负压强在广义相对论里的引力效应。用牛顿力学的直觉硬套，两处都会翻车。
💡 关键直觉：弗里德曼方程是一部自动变速箱——负载（组分密度）随车速（标度因子）自动变化，换挡点就是各组分平权的时刻；驾驶者不是谁，而是能量守恒本身。
演算三：临界密度——宇宙的总账本
弗里德曼方程里有一个特殊密度：让空间恰好平坦的密度，叫临界密度。整个宇宙学的组分占比都拿它当分母，值得亲手算一遍：
# 临界密度 rho_c = 3*H0^2/(8*pi*G) from math import pi H0 = 67.4 * 1000 / 3.0857e22 # 每秒 G = 6.674e-11 rho_c = 3 * H0**2 / (8 * pi * G) print(f"临界密度 = {rho_c:.3e} kg/m^3") Om, Ob, OL = 0.315, 0.0493, 0.685 print(f"物质密度 = {Om*rho_c:.2e} kg/m^3") print(f"重子密度 = {Ob*rho_c:.2e} kg/m^3") print(f"真空能密度= {OL*rho_c:.2e} kg/m^3")
运行输出：
临界密度 = 8.549e-27 kg/m^3 物质密度 = 2.69e-27 kg/m^3 重子密度 = 4.22e-28 kg/m^3 真空能密度= 5.86e-27 kg/m^3
临界密度相当于每立方米不到六个质子——宇宙这台放映机的片场比任何实验室真空都空。但就是这份总账本决定了放映机的转速：物质、真空能各占多少，直接写进第 1.4 节那套方程的每一项。注意重子（普通原子物质）只占约百分之五，这个伏笔要到第 4 章才揭晓它意味着什么。
本节要点回顾
- 方程地位：弗里德曼方程是宇宙学版的能量守恒，比哈勃观测早七年出现，是整部编年史的动力系统。
- 临界密度：由哈勃常数与引力常数决定，约每立方米五个质子；密度参数把组分写成无量纲份额。
- 基准模型：纯物质平直宇宙年龄为三分之二哈勃时间，约九十七亿年，容不下最老恒星——这正是引入真空能的观测动机之一。
- 一组分一命运：物质减速、辐射强减速、真空能加速；密度随标度因子的不同幂次衰减是分岔根源。
- 平权即换挡：编年史的段落划分对应组分密度曲线的交点。
- 年龄闸门：任何改动组分的主张都必须重算宇宙年龄并容纳最老天体的年龄。
下一节给动力系统装上温度计：热宇宙的统计物理如何把"多热"翻译成"多久以前"，让编年史的每一帧获得精确的时间码。