2.2 X射线气体：另一种称重方法的证词
2.2 X射线气体：另一种称重方法的证词
本节摘要：星系团的质量大部分不在星系里，而在充满团内空间的千万度热气体中——它发 X 射线。气体靠自身温度提供的压强梯度对抗引力，写下流体静力学平衡方程；测出温度与密度轮廓，就能反推束缚质量。本节用等温 β 模型复算团质量，并清点重子份额：星系团的重子比例与宇宙学比值一致，缺失的份额没有藏进气体，藏进的是不发光的成分。
证人出场：团内介质
肉眼里的星系团是一窝星系；X 射线卫星眼里的星系团是一团泡在千万到上亿度等离子体里的星系。乌呼鲁卫星 1966 年前后的巡天首次发现团内弥漫 X 射线，爱因斯坦天文台随后成象确认：气体质量比团内恒星总质量还高几倍。气体的来历不神秘——成员星系在团内互相穿行、剥离，累积的星际气体落入围内空间，被引力势加热到维里温度。
对查案者来说，这团气体是完美的示踪物兼温度计。它靠轫致辐射发光（自由电子被离子散射减速，辐射软 X 射线），光度密度正比于密度平方，谱线（尤其铁的 6.7 keV 线）给温度与丰度。温度、密度两到手，称重公式就能开动。
勘验记录：流体静力学平衡称重
一团静止（无大规模流动）的气体，每一点上压强梯度与引力平衡：
dP/dr = −ρ_g·GM(r)/r²。
理想气体 P = n kT，代入并把温度轮廓的导数并进来，得到称重公式：
M(<r) = −(k T r / (G μ m_p))·(d ln ρ_g/d ln r + d ln T/d ln r)。
负号的意思：密度向外降、温度向外降（或平），梯度符号自动让质量为正。最常用的工作模型是等温 β 模型，密度 ρ_g = ρ₀[1+(r/r_c)²]^(−3β/2)，温度取常数，质量有闭式解：
M(<r) = 3β·(kT r)/(G μ m_p)·(r/r_c)²/(1+(r/r_c)²)。
import math G, KEV, MPC, MS = 6.674e-11, 1.602e-16, 3.0857e22, 1.989e30 T_keV, beta, rc, r = 8.0, 0.7, 0.25, 1.5 # 温度、beta、核半径、称重半径 kT = T_keV * KEV mu_mp = 0.6 * 1.673e-27 # 平均分子量×质子质量 fac = (r / rc) ** 2 / (1 + (r / rc) ** 2) M_hyd = 3 * beta * kT * r * MPC / (G * mu_mp) * fac / MS print("等温beta模型：T=%.0f keV，beta=%.1f，rc=%.2f Mpc，r=%.1f Mpc" % (T_keV, beta, rc, r)) print(" M(<r) = %.2e Msun" % M_hyd)
输出：
等温beta模型：T=8 keV，beta=0.7，rc=0.25 Mpc，r=1.5 Mpc M(<r) = 9.10e+14 Msun
约 9×10¹⁴ 太阳质量——与上一节透镜秤（10¹⁵ 量级）以及第1章维里秤（1.05×10¹⁵）互相咬合。三杆秤的物理各不相干：透镜用广义相对论，维里用牛顿动力学，流体静力学用气体热压强。三票投出同一个数，这正是"缺失质量是实的"最有力的形态。
数量级复核：重子到底占几成
气体既是证人也是被告——会不会"缺失的质量"其实就是没算够的气体？把账摆开：
ob, om = 0.0493, 0.3153 # Planck 宇宙学参数 print("宇宙重子份额 Omega_b/Omega_m = %.3f" % (ob / om)) print("星系团观测 f_gas 约 0.13-0.16，与宇宙比值同量级")
输出：
宇宙重子份额 Omega_b/Omega_m = 0.156 星系团观测 f_gas 约 0.13-0.16，与宇宙比值同量级
星系团重子占比（气体加恒星约 15%）与全宇宙的重子占比基本一致——团没有特殊癖好，它按宇宙配方攒物质。要让气体补上缺口，得把气体质量虚报六七倍，而气体质量由 X 射线光度与谱线双路测定，虚报空间早已封死。真正的缺失成分必须几乎不发任何波段的光，且主要不与气体纠缠——这份"证词摘要"直接写进第3章的嫌疑画像。
💡 关键直觉：星系团是宇宙物质的"公平抽样"。它引力阱深到几乎不漏重子，所以重子份额约等于宇宙均值；份额之外的八成五，就是暗物质的全宇宙比例在局部的投影。
解读：等温假设有多伤
β 模型的等温假设是最粗的简化。真实团有温度梯度：中心凉、外围热，或反之（冷锋、激波扰动）。忽略 d ln T/d ln r 项会带来 10% 到 30% 的系统偏差，方向不定——所以现代工作用 XMM-牛顿、钱德拉的详细温度轮廓逐环称重，或干脆用透镜做基准交叉标定。有意思的是，若干团显示流体静力学质量比透镜质量低 10% 到 20%，缺口被归因于湍流与非热压强支持——这不是暗物质的危机，而是气体"体重秤"的已知系统误差，测量越细，账越清。
变式：从团到群再到星系
同样的方法下推到小系统：星系群温度约 1 keV，质量 10¹³ 量级；再往下，单星系晕的气体太稀太冷，X 射线称重让位给旋转曲线。随系统变"小"，重子占比一路走低——矮星系里跌到几个百分点。这个趋势在结构形成理论里有名有姓（详见第5章：小晕形成早、反馈吹掉重子更狠），反过来说明重子占比的变化本身不是误差，而是物理。
完整案例：英仙座团的解剖报告
把方法走一遍全过程。英仙座团（Abell 426）是北天最亮的团级 X 射线源，距离约 2.4 亿光年。第一步测温度：钱德拉的谱仪逐环读铁线与连续谱，核心区约 4 keV，外围爬到 7 到 9 keV——中心那口"凉井"正是 AGN 空腔的证据（黑洞喷流把气体凿出两个直径各两万光年的空洞，空腔的 pV 功率还能反过来标定喷流功率，这是一段与本案无关却同样精彩的花絮）。第二步拟合密度轮廓：β 模型给 rc 约 0.25 Mpc、β 约 0.65。第三步代入称重公式逐环积分，r = 1 Mpc 处 M 约 6×10¹⁴ Msun，r = 2 Mpc 处约 1.5×10¹⁵。第四步对账：气体约 8×10¹³、恒星约 1×10¹³，合计重子一成半，其余八成五没有 X 射线、没有星系光、没有 21 厘米信号——流体静力学秤、透镜秤、维里秤三票一致。整个流程没有任何一步需要"暗物质"这个假设参与测量，它只在最后对账时作为差额出场——这就是旁证比嫌疑人自述更有力的原因。
本节要点回顾
- 气体即证人：千万度等离子体发 X 射线，温度与密度轮廓直接反推束缚质量。
- 三秤咬合：流体静力学、维里、透镜三种独立物理给出同量级团质量（约 10¹⁵ Msun）。
- 重子封顶：团内重子约 15%，与宇宙重子份额一致；缺失成分不可能是气体。
- 系统误差有名有姓：等温假设与湍流支持带来 10% 到 30% 偏差，属方法学而非证据危机。
第二位证人的证词与第一位严丝合缝。下一节是一场真实的"对撞实验"——子弹星系团，让质量与气体当场分家。