1.2 星系盘的引力账本：从亮度分布推算转速

1.2 星系盘的引力账本：从亮度分布推算转速

本节摘要：面对旋转曲线异常，第一个辩解是"恒星质量被低估了"。本节把旋涡星系的指数盘亮度分布完整过账：表面亮度随半径指数衰减，恒星质量的内含分数有解析式，配上质光比即可算出盘的引力贡献与预期转速曲线。演算表明，无论质光比取多慷慨，盘贡献的曲线在中外区必然下弯，与实测平曲线的形状从根上不合。

证词交叉：亮度是能看到的一切

上一节钉死了"缺口存在"，本节处理第一条辩解。辩方的逻辑不难理解：质光比（质量与光度的比值）本来就有弹性——矮星里塞满暗弱红矮星，质光比能到 10 以上；年轻蓝星团里只有 2 上下。把银河系盘的质光比往高了估，恒星质量翻一两倍不是不可能，缺口不就小了吗？

这条辩解的软肋在形状而不在幅度。盘的贡献有个硬约束：指数盘。旋涡星系的表面亮度沿半径按 I(r) = I₀ exp(−r/R_d) 衰减，标长 R_d 典型 2 至 5 kpc。盘是薄盘子不是球，质量分布向中心集中的程度远不如"全部装在中心"，但十倍标长之外，账上的恒星也确实数得干干净净——光度测到哪儿，账本就到哪儿。

勘验记录：指数盘的完整对账

指数盘内含质量的分数有解析解。把面密度对半径积分，得到半径 r 之内占总质量的比例：

f(x) = 1 − (1 + x)·e^(−x)，其中 x = r/R_d。

这个式子是本节的主角：盘的"聚财速度"随半径飞快衰减。到 4 倍标长，账面已收进 98.6%；发光盘的光学边缘大致就在 3.2 至 4.4 倍标长处。

import math
G, MS, KPC = 6.674e-11, 1.989e30, 3.0857e19
Rd, Mstar = 3.5, 5e10          # 标长 kpc；恒星总质量 Msun（已按高质光比从宽）
print("指数盘内含质量分数 f = 1-(1+x)e^-x：")
for r in [3.5, 8, 14, 21, 35]:
    x = r / Rd
    f = 1 - (1 + x) * math.exp(-x)
    M = Mstar * f
    v = math.sqrt(G * M * MS / (r * KPC)) / 1e3   # 球对称近似下的转速
    print("  r=%4.1f kpc  f=%.3f  M(<r)=%.2e Msun  v_pred=%.0f km/s（实测~220）"
          % (r, f, M, v))
​

输出：

指数盘内含质量分数 f = 1-(1+x)e^-x：
  r= 3.5 kpc  f=0.264  M(<r)=1.32e+10 Msun  v_pred=127 km/s（实测~220）
  r= 8.0 kpc  f=0.666  M(<r)=3.33e+10 Msun  v_pred=134 km/s（实测~220）
  r=14.0 kpc  f=0.908  M(<r)=4.54e+10 Msun  v_pred=118 km/s（实测~220）
  r=21.0 kpc  f=0.983  M(<r)=4.91e+10 Msun  v_pred=100 km/s（实测~220）
  r=35.0 kpc  f=1.000  M(<r)=5.00e+10 Msun  v_pred=78 km/s（实测~220）
​

两条信息直接读出来。其一，8 kpc（太阳圈）处盘的贡献撑到头也就 130 到 150 km/s——球对称近似还会略低估薄盘的贡献，即便按所谓"最大盘"修正拉到 180 km/s，离 220 仍差一大截。其二，更致命的是形状：14 kpc 之后 f 几乎不再增长，预言转速单调下弯，78 对 220，越往外缺口越大。

数量级复核：质光比需要多离谱

再做一道反算题，让辩方自己报数。要单靠盘在 35 kpc 处维持 220 km/s，需要 M(<35 kpc) ≈ 4×10¹¹ Msun。盘的光度约 1×10¹⁰ L_sun（银河系口径），质光比得开到 40——是正常旋涡星系值（1 到 3）的十几倍。恒星族群模型里塞不进这么多暗物质星，褐矮星与行星的巡天也把这类"暗重子"封了顶（微透镜排查详见第4章 MACHO 一节）。

顺带把气体入账。中性氢加分子气体大约再添 10% 到 20% 重子，而且气体盘比恒星盘更延展，对中外区曲线略有帮助——但量级摆在那里，填不平九倍缺口。账本三栏（恒星、气体、残差）列清楚：

⚠️ 常见坑：把"质光比不确定"当成万能补丁。质光比是幅度旋钮，救不了形状错位——指数盘预言的中外区下弯是结构性的，与平曲线的分歧随半径拉大。

变式：拿这把尺子去量别的星系

同样的账法可以套到任何有测光的旋涡星系。低表面亮度星系（LSB）是最锋利的试金石：它们标长更长、面亮度更低，恒星盘更轻，平曲线却照样平——按盘算质量的预期曲线跌得更惨，暗物质占比反而更高，矮星系里重子连一成都占不到（第 3 章的矮星系解剖有演算）。反过来，近红外测光能穿透尘埃把老恒星数全，J 波段质光比稳定在 0.8 上下，把"藏起来的亮物质"这条退路也堵得差不多。

三问三答

问：质光比 40 不行，那 10 呢？10 已经够离谱了吧。
答：10 也不够，而且方向不对。质光比抬到 10，盘的质量翻倍，太阳圈转速也就从 134 涨到 190 km/s 上下；可 35 kpc 处的账目封顶在 5×10¹⁰ 乘以 10 也只有 5×10¹¹ 的一半不到——最要命的是形状照样下弯。审这个辩解要盯形状，别被幅度扯走。

问：核球与厚盘呢？会不会藏了没数进来的成分？
答：核球在 3 kpc 内就收完了，对外区曲线零贡献；厚盘质量约为主盘的一到两成、标长只长一点，改不了 20 kpc 之外的结论。每一笔"藏质量"的提案都得回答同一个问题：它能让 M(<r) 随半径线性增长吗？目前没有任何已知恒星成分能做到。

问：会不会远盘外有大量极暗的气体云，测光看不见？
答：气体不是测光看不见，是射电看得见——21 厘米氢线把中性氢数到了比恒星盘更远的地方，账上只有恒星盘的一到两成。这条退路在 1.3 节会正面审过。

本节要点回顾

• 盘有解析账本：f(x) = 1−(1+x)e^(−x)，四倍标长收进 98.6%，之后恒星账目封顶。
• 幅度不够：即便最大盘假设，太阳圈处盘贡献约 180 km/s，与 220 的差额就是暗物质的份额。
• 形状更不够：中外区预言单调下弯（35 kpc 处 78 km/s），与实测 220 km/s 的平直南辕北辙。
• 质光比救场失败：补平 35 kpc 曲线需质光比 40，恒星族群与微透镜观测双重排除。

恒星的账审完了，干净。下一节把观测手段升级到射电波段，看看发光边缘之外的世界——曲线还平吗？