5.4 小尺度危机：尖点缺失与卫星不足
5.4 小尺度危机：尖点缺失与卫星不足
本节摘要：冷暗物质在星系团尺度的胜利记录没有争议，可把分辨率调到矮星系，模拟与观测就连吃三张黄牌：模拟的暗晕中心密度该有"尖点"而观测偏"平核"；模拟造的上万子晕远多于看得见的卫星星系；最大子晕的速度比最亮的矮星系高出太多——"太大而不能倒"。三道题各有辩护席（重子反馈、温暗物质、自相互作用暗物质），但没有任何一方拿到全票。这是悬案至今未结的核心战场。
三道黄牌从哪来
先复述控方立场。冷暗物质（CDM）模拟给出的暗晕密度剖面遵循 NFW 形式（第3章）：内区 ρ ∝ r⁻¹，中心是"尖点"。同时，一个银河系质量的暗晕按等级成团必然携带海量子晕——质量函数按幂律一路延伸到分辨率极限。控方由此提出三条指控：
- 尖点-核心问题（cusp-core）：观测到的矮星系旋转曲线（第1章的老朋友）在中心几百秒差距趋于平坦，对应密度核心，而 NFW 尖点给出的内区质量明显超出运动学允许值。
- 卫星缺失问题（missing satellites）：模拟预言银河系应有一万个以上可分辨子晕，观测到的明亮卫星只有几十个。
- 太大而不能倒（too-big-to-fail）：模拟中最大的子晕质量太大，不可能没形成恒星——却在你看得见的地方找不到对应物；它们的速度上限（Vmax 30–40 km/s）远超最亮矮星系的 12 km/s 左右。
三条指控里第三条最狠：卫星缺失还能辩解"小晕没点亮"，太大而不能倒堵死了这条路——按重子物理，这些大子晕没有理由不发光。
图：冷暗物质的三道小尺度黄牌——模拟预期与观测的对质
二、把差距算成数字
指控不能停在形容词。拿一条 M200 = 1e10 Msun、浓度 c = 12 的矮星系暗晕，分别按 NFW 尖点与 α = 1 kpc 的平核剖面重算内区质量，看看"尖"与"平"差出多少个太阳：
import math G = 4.302e-6 # kpc (km/s)^2 / Msun h = 0.7 rhoc = 2.775e11 * 0.25 * h * h * 1e-9 # 临界密度 Msun/kpc^3 M200, c = 1e10, 12.0 r200 = (3 * M200 / (4 * math.pi * 200 * rhoc)) ** (1.0 / 3.0) rs = r200 / c f = lambda x: math.log(1 + x) - x / (1 + x) rho_s = M200 / (4 * math.pi * rs ** 3 * f(c)) print("r200 = %.1f kpc, r_s = %.2f kpc, rho_s = %.4f Msun/pc^3" % (r200, rs, rho_s * 1e-9)) mnfw = lambda r: 4 * math.pi * rho_s * rs ** 3 * f(r / rs) vmax = max(math.sqrt(G * mnfw(rs * i / 500) / (rs * i / 500)) for i in range(1, 1001)) print("NFW: M(<1 kpc) = %.2e Msun, Vmax = %.0f km/s" % (mnfw(1.0), vmax)) # 平核剖面 rho = rho0 / (1 + r/rc)^2, rc = 1 kpc，同 M200 归一 rc = 1.0 def mcore(r, rho0): n = 20000 s = r / n m, rr = 0.0, s / 2 for _ in range(n): m += 4 * math.pi * rr ** 2 * rho0 / (1 + rr / rc) ** 2 * s rr += s return m rho0 = M200 / mcore(r200, 1.0) print("cored rho0 = %.4f Msun/pc^3" % (rho0 * 1e-9)) for r in (0.3, 1.0): print("r=%.1f kpc: NFW %.2e vs cored %.2e Msun" % (r, mnfw(r), mcore(r, rho0))) print("300 pc 内质量比 NFW/cored = %.1f 倍" % (mnfw(0.3) / mcore(0.3, rho0)))
输出：
r200 = 70.5 kpc, r_s = 5.88 kpc, rho_s = 0.0024 Msun/pc^3 NFW: M(<1 kpc) = 7.14e+07 Msun, Vmax = 31 km/s cored rho0 = 0.0126 Msun/pc^3 r=0.3 kpc: NFW 7.42e+06 vs cored 9.59e+05 Msun r=1.0 kpc: NFW 7.14e+07 vs cored 1.81e+07 Msun 300 pc 内质量比 NFW/cored = 7.7 倍
7.7 倍——这就是尖点与核心之争的量化形态。旋转曲线测的是中心 300 pc 内的 enclosed mass，观测值系统性靠近平核一侧；若 NFW 为真，最亮的矮星系早该被自己的暗晕扯散。顺带注意 Vmax = 31 km/s 这个数：一个 1e10 Msun 的暗晕恰恰就是"太大而不能倒"阵营的成员——它本该点亮一颗 Vmax 20–30 km/s 的 dwarf，而 Fornax、Sculptor 这些最亮的伙伴只有 12 km/s 上下。
三、辩护席上的三个方案
方案一：重子反馈（在标准框架内自救）。 超新星与恒星风把气体吹出矮晕中心，掏空引力势阱，暗物质随之向外重排——尖点被"熨"成核心。这是 ΛCDM 阵营的主辩护：三道黄牌共享同一个嫌犯，即"暗物质-only 模拟漏掉了重子物理"。支持证据：反馈强的星系恰好核大；反对意见：最暗的 dwarf（Vmax < 20 km/s）几乎不形成恒星，反馈掏不动它的势阱，可它们照样是核——辩护在最低质量端失效。
方案二：温暗物质（改粒子）。 keV 级的惰性中微子（第4章 4.3）用自由流抹掉最小的子晕，卫星数量与太大而不能倒同时缓解。代价：Lyman-alpha 森林把质量下限推到 5 keV 以上，缓解幅度被压到"治标"级别——子晕数目降下来了，幸存大子晕的内区剖面照旧是尖的。温暗物质治得了第二张黄牌，治不透第一张。
方案三：自相互作用暗物质（SIDM，改行为）。 让暗物质粒子之间存在截面约 1 cm²/g 的弹性散射。晕中心密度高，粒子频繁互撞，把尖点"热化"成核心；外围密度低，碰撞罕见，大尺度行为与 CDM 无异——精准打击第一张黄牌。定量检验一下这个截面的物理含义：
rho_kg = 0.1 * 2e30 / (3.0857e16) ** 3 # 晕中心 0.1 Msun/pc^3 换 kg/m^3 sig_m = 1e-4 / 1e-3 # 1 cm^2/g = 0.1 m^2/kg lam = 1.0 / (rho_kg * sig_m) # 平均自由程 m print("晕中心密度 = %.2e kg/m^3" % rho_kg) print("sigma/m = 1 cm^2/g 时平均自由程 = %.1e m = %.0f kpc" % (lam, lam / 3.0857e19)) age_s = 1.38e10 * 3.156e7 v_typ = 2e5 # 中心区典型速度 200 km/s print("宇宙年龄内互撞次数 ~ %.0f 次（中心区，v=200 km/s）" % (age_s * v_typ / lam))
输出：
晕中心密度 = 6.81e-21 kg/m^3 sigma/m = 1 cm^2/g 时平均自由程 = 1.5e+21 m = 48 kpc 宇宙年龄内互撞次数 ~ 59 次（中心区，v=200 km/s）
一个粒子在矮晕中心的一生要互撞约 60 次——足以把速度分布"搅拌"成各向同性的等温核，这正是 SIDM 熨平尖点的机理。而在星系团中心，密度低一个量级、速度高五倍，净散射率反而更低，大尺度的轨道结构与子弹星系团证据得以保全（第2章）。SIDM 的尴尬在于需要截面随速度变化：矮星系要求 ~10 cm²/g，星系团压到 <1 cm²/g，粒子物理模型得专门为此搭台阶。
四、悬而未决的判词
三张黄牌的现况可以一句话概括：反馈方案在亮度居中的矮星系上节节胜利，却输在最低质量端；温暗物质与 SIDM 各治一张牌，没有单一方案全票通过；而"太大而不能倒"在计入检测偏差后是否还成立，2010 年代吵了整整十年，至今没有终审。对暗物质悬案而言，这里正是翻案申请（第6章）的翻供入口——替代引力理论的第一句台词，永远从小尺度开始。
本节要点回顾
- 三道黄牌：尖点-核心（300 pc 内质量差 7.7 倍）、卫星缺失（上万子晕对几十个亮卫星）、太大而不能倒（大子晕 Vmax 30–40 对矮星系 12 km/s）。
- 量化对质：旋转曲线测内区质量，观测偏平核；NFW 尖点在最低质量端与运动学冲突。
- 三个辩护：重子反馈（主流，低质量端失效）、温暗物质（治卫星数不治剖面）、SIDM（1 cm²/g 热化中心，需截面随速度演化）。
- 战局：无全票方案，小尺度是悬案最后的开放战线。
案情重构到此结束。暗物质模型赢了大局、输了几场小仗；它到底是谁，还得回到实验台上对质。下一章开庭审讯——三条探测路线与一份翻案申请。