3.3 暗物质晕解剖：NFW剖面与质量浓度

3.3 暗物质晕解剖：NFW剖面与质量浓度

本节摘要：暗物质在星系尺度上组织成球状晕，数值模拟给出其普适密度剖面——NFW 形式：ρ ∝ 1/[r(1+r/r_s)²]，内区尖点、外区陡降，由特征半径 r_s 与归一化两个参数确定。习惯上改用总质量 M200 与浓度 c = r200/r_s 描述。本节从定义反解银河系晕的参数，复算内区质量，并说明浓度-质量关系如何成为连接观测与模拟的桥——尖点行为则为第5章小尺度危机埋下引线。

形态学入门：晕是什么形状

星系不是"一团暗物质裹着一个发光盘"这么简单，晕有自己的形状语言。大到银河系这样的旋涡星系，暗物质晕是近球形（略扁或三轴），延伸到两三百 kpc；盘、核球、球状星团团、卫星星系都泡在里面。晕不是一个刚体，而是等温的粒子云：粒子各自在椭圆轨道上翻滚，整体没有清晰边界，天文学只好用"平均密度降到宇宙临界密度某倍数"来人为划界。

划界的行规是 M200 与 r200：平均密度等于 200 倍临界密度的球内质量与半径。200 这个数不神秘——对应坍缩时间的球模型恰好给出几十到几百倍，行内取整而已。在这套语言里，一个晕由两个数完全刻画：质量 M200 与浓度 c ≡ r200/r_s。

NFW 剖面来自纯引力成团的普适性：纳瓦罗、弗伦克与怀特 1996 到 1997 年在当时的模拟中发现，从星系到星系团跨四个量级质量的晕，密度剖面套进同一个双参数公式竟都能对上：

ρ(r) = ρ_s / [x(1+x)²]，x = r/r_s。

内区 x ≪ 1 时 ρ ∝ 1/r（尖点，cusp），外区 x ≫ 1 时 ρ ∝ 1/r³（陡降）。对剖面积分得内含质量 M(<r) = 4πρ_s r_s³ [ln(1+x) − x/(1+x)]。

图 NFW 剖面与观测核心的分歧

动手验算：给银河系晕做解剖

**第一步：从 M200 解出 r200。**定义式 4πr200³/3 × 200ρ_c = M200，直接开立方：

import math
G, KPC, MS = 6.674e-11, 3.0857e19, 1.989e30
H0 = 67.36e3 / (3.0857e22)
rho_c = 3 * H0 ** 2 / (8 * math.pi * G)     # 今日临界密度
M200, c = 1e12, 10                          # 银河系量级质量与浓度
r200 = (3 * M200 * MS / (4 * math.pi * 200 * rho_c)) ** (1 / 3) / KPC
rs = r200 / c
def fNFW(x):
    return math.log(1 + x) - x / (1 + x)
rho_s = M200 / (4 * math.pi * (rs * KPC) ** 3 * fNFW(c))      # Msun/m^3
rho_kpc = rho_s * KPC ** 3
print("M200=1e12 Msun，c=10：r200=%.0f kpc，r_s=%.1f kpc" % (r200, rs))
print("rho_s = %.2e Msun/kpc^3 = %.4f Msun/pc^3（约 %.2f GeV/cm^3）"
      % (rho_kpc, rho_kpc / 1e9, rho_kpc / 1e9 * 38.0))
​

输出：

M200=1e12 Msun，c=10：r200=212 kpc，r_s=21.2 kpc
rho_s = 5.64e+06 Msun/kpc^3 = 0.0056 Msun/pc^3（约 0.21 GeV/cm^3）
​

太阳位置（8 kpc ≈ 0.38 r_s）的局部密度算出来约 0.3 到 0.4 GeV/cm³，与 3.1 节的动力学共识值衔接——模型自洽的小小交叉验证。

**第二步：尖点往里走质量怎么攒。**小 x 时 f(x) ≈ x²/2，于是 M ∝ r²——尖点的意思是越往里质量攒得越狠：

import math
M200, c = 1e12, 10
rs = 21.2                               # 接上块：r200=212 kpc 除以 c=10
def fNFW(x):
    return math.log(1 + x) - x / (1 + x)
for r_in in [0.1, 0.3, 1.0]:
    x = r_in / rs
    print("M(<%.1f kpc) = %.2e Msun" % (r_in, M200 * fNFW(x) / fNFW(c)))
​

输出：

M(<0.1 kpc) = 7.43e+06 Msun
M(<0.3 kpc) = 6.60e+07 Msun
M(<1.0 kpc) = 7.03e+08 Msun
​

银心一 kpc 内按 NFW 该有七亿太阳质量的暗物质。观测（核星团的动力学与微透镜）支持这个量级——但换到矮星系，内区质量被卡到只有 NFW 预言的十分之一以下，尖点问题就在那里爆发（第5章 5.4 节）。

解读：浓度参数是宇宙学的温度计

浓度 c 不是任性参数，它记录晕的形成历史：早形成的晕在宇宙更致密的时代坍缩，攒得更紧，浓度更高。模拟给出干净的规律——质量越大，浓度越低（大晕晚成），红移越高，同质量晕浓度越高（宇宙当年更挤）。银河系量级 c 约 10 到 15，星系团级 c 约 4。观测测得的浓度-质量关系与 ΛCDM 模拟预言吻合到几十个百分点，这是模型在"晕形态学"上拿到的关键分。若未来巡天测出浓度关系系统性偏移，那就是悬案翻供的线索。

变式：把解剖刀用在矮星系与子晕上

同一套公式缩小规模即可解剖矮星系晕（M200 约 10¹⁰ Msun，c 可到 20 以上，内区尖点更显眼），也能算子晕的最大圆速度（5.4 节的"too big to fail"要用）。NFW 之外还有拟合族的竞争：观测者爱用带核心的 Burkert 剖面，模拟者坚持尖点；两者的分歧不是技术细节，而是冷暗物质在小尺度的正面审判——那是第5章的主戏。

两问两答

问：NFW 是从哪来的"普适"？会不会只是拟合函数碰巧通用？
答：它是模拟的涌现结果而非输入假设——无论初始条件、宇宙学参数怎么调，引力成团后的晕都收敛到这个形式，原因是"球壳坍缩"过程的渐近自相似：内区由最早坍缩的物质铺成，外区由晚期吸积主导，两层机制衔接自然给出 1/[x(1+x)²] 的形状。它对流体或自相互作用的暗物质会变（那是第5章的辩点），对"纯引力无碰撞"的假设则相当稳健。

问：浓度 c 测错会怎样？它对探测实验有影响吗？
答：影响是直接的。直接探测的局部密度与速度分布、间接探测的 J 因子、微透镜的光学深度，全都押在晕的结构参数上。银河系 c 的不确定度（8 到 15）给局域密度带来约两成的摆动，实验限值曲线的换算随之浮动——所以晕形态学看似"纯理论"，其实是探测实验的标定环节。

本节要点回顾

• 普适剖面：NFW 形式 ρ ∝ 1/[x(1+x)²]，两参数（M200 与 c）刻画跨四个量级质量的晕。
• 银河系解剖：M200=10¹²、c=10 给出 r200=212 kpc、r_s=21 kpc、ρ_s≈0.006 Msun/pc³。
• 尖点含义：内区 M ∝ r²，一 kpc 内七亿太阳质量；矮星系观测与此冲突。
• 浓度即历史：c 与质量反相关、与红移正相关，观测吻合度是 ΛCDM 的形态学分。

晕的形态清楚了。剩下的问题是：这些粒子从哪儿来、为什么今天恰好这么稀？下一节回到宇宙最初三分钟之后，看"冻结"怎么定下今天的丰度。