3.3 动力学质量测量

3.3 动力学质量测量

本节摘要：运动学数据的终极用途是称出星系质量。本节从圆轨道近似出发，介绍球对称系统的 Jeans 方程与维里定理，讨论质量-光度比的量级差异，并给出由速度弥散推算椭圆星系质量的 Python 算例，最后对比引力透镜、X 射线气体等独立手段的交叉验证。

三个层次的公式

第一层：圆轨道近似。 对盘星系，假定气体沿圆轨道运动，半径 r 内包含的总质量满足

M(r) ≈ V(r)² · r / G

其中 G = 4.30×10⁻⁶ kpc·(km/s)²/M☉。把平坦旋转曲线 V ≈ 220 km/s 代入 r = 50 kpc，得到 M ≈ 5.6×10¹¹ M☉，而该半径内恒星加气体的质量往往只有百分之十几——差异就是暗物质晕。

第二层：球对称 Jeans 方程。 对椭圆星系、核球和晕这类弥散支撑系统，旋转曲线不再适用。假设系统稳态、球对称，且速度各向同性，Jeans 方程退化为简单的"维里近似"：轨道平均意义下，动能与引力势能满足 2T + W = 0。引入质量剖面模型（如 NFW），观测到的恒星数密度与速度弥散分布就能约束模型参数。对点质量近似，一维速度弥散 σ 与质量的关系近似为 M(r) ≈ r·σ²/G，这是"把弥散当速度"的快速估计。

第三层：完整动力学建模。 更严谨的做法是用分布函数或轨道叠加模型拟合全部运动学数据（包括高阶速度矩、视线速度分布形状），反演出引力势 Φ(r)，再与暗物质模型对照。这类"质量过剩"分析（mass excess）是目前确定暗物质晕形状、倾角与各向异性的标准方法。

质量-光度比的警示

把动力学质量除以恒星光度（或恒星质量）就得到质量-光度比 M/L。对正常恒星系统，M/L 在可见光波段为 2 到 6（单位 M☉/L☉），而星系整体观测值往往达到 20 到 100，甚至更高。这个量级差异是暗物质存在的又一独立标尺。它提醒我们：星系照片上看到的只是冰山一角，光度学的"瘦身"与动力学的"发福"之间的鸿沟，正是暗物质刻画线。

维里定理与星系团质量

尺度再放大到星系团，维里定理成为主角。对自引力束缚、处于统计平衡的系统，势能 W 与动能 K 满足 2K + W = 0。把星系团成员星系看作无碰撞粒子，测得视线速度弥散 σ，就有

M ≈ 3σ²R_h/G

其中 R_h 为半质量半径。对典型富星系团 σ ≈ 800-1000 km/s、R_h ≈ 1 Mpc，得到质量约 10¹⁵ M☉。而成员星系的总恒星质量只有约 10¹³ M☉，X 射线气体约 10¹⁴ M☉，其余全是暗物质。这个"重子与暗物质之比仅约 1:9"的结论，与宇宙整体物质构成（重子约占物质总量的 16%）完全一致，构成暗物质存在的又一独立证据。

维里定理使用时要警惕偏差。正在进行并合或尚在坍缩的子结构（团内亚结构）会系统性低估或高估质量；投影方向不同也会带来因子级的误差。现代工作普遍用模拟校准的缩放关系（如质量-速度弥散关系、质量-X射线温度关系）来代替简单的维里估计，误差可控制在百分之二十以内。

独立手段交叉验证

动力学质量并非孤证，四种独立测量互相印证。引力透镜：前景星系或星系团的引力弯曲背景星系的光，弧的偏折角直接由投影质量决定，与动力学状态无关。X 射线气体：星系团热气体的静水平衡给出质量剖面，典型团质量达 10¹⁴ 到 10¹⁵ M☉，其中暗物质占八成以上。卫星星系运动学：银河系的卫星星系速度弥散约 100 km/s，给出银河系维里质量约 10¹² M☉。宇宙微波背景与弱透镜巡天：从统计意义上标定暗物质晕质量函数。四类方法在几十 kpc 到几 Mpc 尺度上高度一致，说明暗物质确实存在且主导引力。

Python 算例：从速度弥散估计椭圆星系质量

import numpy as np

G = 4.30e-6   # kpc*(km/s)^2/Msun

def virial_mass(sigma_km_s, R_eff_kpc, alpha=4.0):
    """
    用维里定理估算球对称系统质量。
    alpha 为结构因子：点质量 1，奇异性等温球约 2.5-4。
    M ≈ alpha * R_eff * sigma^2 / G
    """
    return alpha * R_eff_kpc * sigma_km_s**2 / G

# 典型近邻椭圆星系 NGC 1399
sigma = 320.0        # 中心速度弥散 km/s
Re = 3.5             # 有效半径 kpc
M = virial_mass(sigma, Re)
print(f"NGC 1399 动力学质量 ≈ {M:.2e} Msun")

# 与光度对比：V 波段绝对星等 Mv = -22.1
Mv = -22.1
Lsun_v = 10.0**(0.4 * (4.83 - Mv))   # 4.83 为太阳 V 绝对星等
print(f"恒星光度 L ≈ {Lsun_v:.2e} Lsun")

# 若恒星质量-光度比 M/L = 4，则恒星质量远小于动力学质量
M_star = 4.0 * Lsun_v
print(f"恒星质量 ≈ {M_star:.2e} Msun，暗物质占比 ≈ {(1 - M_star/M)*100:.1f}%")

# ---- 附加：弱引力透镜质量估计（简化圆对称模型）----
# 弱透镜由背景星系的剪切信号 gamma_t 反推柱质量密度，此处给示意式
# 临界表面密度 Sigma_crit 由几何决定：Sigma_crit = c^2/(4*pi*G) * Ds/(Dl*Dls)
c_km_s = 3.0e5                 # 光速 km/s
Ds, Dl, Dls = 1000.0, 400.0, 600.0   # 角直径距离，单位 Mpc，示意值
Sigma_crit = c_km_s**2 / (4 * np.pi * G * 1e3) * (Ds / (Dl * Dls))
print(f"临界表面密度 Σ_crit ≈ {Sigma_crit:.2e} Msun/kpc^2")

# 若在 r=100 kpc 处测得切向剪切 gamma_t=0.03，则表面质量
r_l = 100.0
gamma_t = 0.03
M_proj = np.pi * r_l**2 * gamma_t * Sigma_crit
print(f"弱透镜投影质量 M(<100kpc) ≈ {M_proj:.2e} Msun")
print("与动力学估计对比：两类独立方法在误差内一致")
​

从观测数据到质量剖面

实际工作中，从原始数据到质量剖面要经过严格的数据链路。以椭圆星系为例：先由深场图像做二维亮度分布拟合，扣除点扩散函数后得到恒星表面亮度分布 I(R)；再由积分场光谱得到每个空间元的速度分布，拟合高斯轮廓提取平均速度与速度弥散 σ(R)。有了 I(R) 与 σ(R)，加上恒星的速度各向异性假设（通常用 β(r) 描述径向与切向弥散之比），代入球对称 Jeans 方程

d(νσ²)/dr + 2βνσ²/r = -ν · GM(r)/r²

其中 ν 是恒星数密度（由 I(R) 反阿贝尔变换得到），即可迭代求解质量剖面 M(r)。这个方程说明一个朴素的物理直觉：恒星数密度越陡、弥散越大，需要的质量就越多。

关键的系统误差来自投影与各向异性简并。我们只能测量视线方向的速度，若恒星轨道偏径向，同样的 σ(R) 会对应更小的质量；若偏切向，则对应更大的质量。打破简并的办法有三：其一，测量速度弥散的四阶矩（kurtosis），径向轨道系统的高阶矩与切向轨道系统不同；其二，同时使用多种示踪体（行星状星云、球状星团、卫星星系），它们从中心到外圈接力覆盖不同半径；其三，用 Schwarzschild 轨道叠加模型穷举轨道库，直接拟合二维运动学场。目前对近邻椭圆星系的建模已能做到把质量剖面约束到百分之十几的精度。

质量剖面与暗物质晕参数

把动力学质量剖面与暗物质模型对照，就得到晕的结构参数。对 NFW 轮廓 ρ(r)=ρ₀/[(r/r_s)(1+r/r_s)²]，外推或内插观测数据可解出尺度半径 r_s 与浓度参数 c = r_vir/r_s。对低红移盘星系，典型浓度 c 约 8 到 15；浓度随质量缓慢下降，随形成红移升高而增大——形成越早、宇宙平均密度越高，晕越致密。矮星系的情形尤其有趣：它们质量小、引力浅，动力学质量几乎完全被暗物质主导，旋转曲线或弥散剖面直接"测量"暗物质晕本身。因此矮星系是检验暗物质模型（包括核心-尖点之争、自相互作用暗物质、暖暗物质）的最干净实验室。

质量测量的前沿问题

动力学质量测量当前面临三类挑战。一是投影退化：我们只看到视线方向的速度分量，轨道各向异性与质量分布存在简并，需要高阶矩或多色示踪体打破。二是非平衡态：并合中的星系、亚结构丰富的星系团并不满足稳态假设，维里定理会系统性偏差，必须用数值模拟校准。三是修改引力理论的竞争：MOND 等修改引力方案能否在不引入暗物质的情况下重现全部观测，仍是活跃争论。未来，SKA 把 21 厘米旋转曲线推广到高红移，JWST 与三十米级望远镜则把恒星速度弥散测量推进到 z≈2 的星系，质量-光度比的"暗物质缺口"将被系统追踪到宇宙年龄只有二十亿年的时候。

需要强调的是，动力学质量测量并不局限于星系本身。银河系内部，恒星围绕中心超大质量黑洞的克卜勒轨道、球状星团的内部动力学、以及奥尔特云的外围示踪体，共同把"称重"这件事从局域一直延伸到宇宙学尺度。每一次把速度换成质量的换算，都是牛顿力学在极端尺度上的一次检验；任何系统性偏离，都可能指向对引力理论本身的修改。因此本节讲述的三种公式——圆轨道近似、Jeans 方程、维里定理——不仅是工具，更是物理学家审视"物质与时空"关系的基本试金石。