跳转到内容
← 返回知识库
星系天文学约 20 分钟阅读

暗物质晕与星系旋转曲线

暗物质晕旋转曲线NFW轮廓Vera RubinWIMP轴子直接探测修正引力

概述

"宇宙中大多数质量是不可见的"——这是 20 世纪天文学最重大的发现之一。暗物质不与电磁力耦合（因此不发光、不反射光，对光学、X 射线、射电观测均"透明"），但通过引力效应无处不在地影响着宇宙的结构。今天我们知道，暗物质约占宇宙总能量密度的 27%，远超普通重子物质的 5%（Planck 2018）。星系旋转曲线是识别暗物质最清晰、最直接的观测工具。

旋转曲线：最直接的暗物质证据

开普勒式下降与观测事实

如果星系质量全部集中在可见的恒星和气体中，外区星系的旋转速度 $v(r)$ 应按照类似行星绕太阳运动的规律下降——即开普勒式下降：

v(r)=GM(r)r→r→∞∝r−1/2v(r) = \sqrt{\frac{GM(r)}{r}} \xrightarrow[r \to \infty]{} \propto r^{-1/2}

但自 1960 年代末起，天文学家用 21 厘米氢线（中性氢 HI）追踪旋涡星系外区（远超恒星盘的区域）的旋转速度，发现旋转曲线保持平坦，而非下降：

v(r)≈const∼100–300 km/sv(r) \approx \text{const} \sim 100 \text{--} 300 \text{ km/s}

平坦旋转曲线在距银心 15–80 kpc 的区域均被观测验证（对于银河系，利用 HI 21 厘米线和 Cepheid 变星等示踪体；对于外部星系，利用 HI 和 HαH\alpha 发射线成像）。

Vera Rubin 的历史贡献

天文学家 Vera Rubin（1928–2016）和 Kent Ford 在 1970 年代系统测量了多个旋涡星系（最初以仙女座 M31 为目标，随后扩展至 60 多个星系）的光学旋转曲线，为平坦旋转曲线提供了不可辩驳的统计证据（Rubin et al. 1978, 1980）。尽管早期研究（包括 Zwicky 1930 年代关于星系团质光比的工作，和 Bosma 1978 的 HI 旋转曲线博士论文）已有相关发现，Rubin 的系统性工作是促使天文学界正视"暗物质问题"的关键节点。

平坦旋转曲线要求 M(r)∝rM(r) \propto r——即星系外区仍有质量随半径线性增长，这只能由延伸至远超可见星系盘的暗物质晕来解释。

从观测策略看，旋转曲线之所以特别有力，是因为它把动力学问题压成一维：$v(r)$ 直接约束半径 $r$ 以内的总质量，几乎不依赖恒星演化模型的细节。光学示踪体（HαH\alpha、行星状星云、球状星团）与射电 HI 盘互补：前者分辨率高、适合内盘与隆起；后者可延伸到光学盘数倍之外，正是"平坦外区"最干净的舞台。对单个星系而言，一旦气体盘足够延展且倾斜改正可靠，晕的存在与大致质量尺度便可从曲线形状读出，而不必先假定粒子候选是什么。

暗物质晕的理论描述

NFW 轮廓

数值模拟（Navarro, Frenk & White 1996, 1997）显示，纯暗物质的 Λ\LambdaCDM 宇宙学模拟所形成的暗物质晕，其密度轮廓有近乎普适的形式——NFW 轮廓：

ρNFW(r)=ρs(r/rs)(1+r/rs)2\rho_\text{NFW}(r) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2}

  • 内区（r≪rsr \ll r_s）：ρ∝r−1\rho \propto r^{-1}（尖峰，cusp）
  • 外区（r≫rsr \gg r_s）：ρ∝r−3\rho \propto r^{-3}
  • 尺度半径 rsr_s 和特征密度 ρs\rho_s 通过浓度参数 c=r200/rsc = r_{200}/r_s 相关（其中 r200r_{200} 是平均密度为临界密度 200 倍处的半径）

对于银河系类型的星系，rs∼20r_s \sim 20 kpc，c∼10c \sim 10–$20$，暗物质晕总质量约 1012M⊙10^{12} M_\odot。

浓度 $c$ 随晕质量与形成红移系统变化：更早坍缩的小晕通常更"浓"，大晕则相对松散。观测上常用弱引力透镜、卫星星系轨道与恒星流扰动来交叉约束 rsr_s 与外区斜率；内区则对重子收缩与反馈更敏感，因此同一条 NFW 公式在拟合时往往要允许内坡可变，或改用 Einasto、cored isothermal 等经验形式。重要的是：晕轮廓是统计普适性加上重子改造后的结果，不是单一解析函数的教条。

星系的完整质量分解

一个旋涡星系的旋转曲线是各组分引力的合力，可以分解为：

vobs2(r)=vdisk2(r)+vbulge2(r)+vgas2(r)+vDM2(r)v_\text{obs}^2(r) = v_\text{disk}^2(r) + v_\text{bulge}^2(r) + v_\text{gas}^2(r) + v_\text{DM}^2(r)

通过对旋转曲线做多组分拟合，可以分别约束恒星质光比、气体分布和暗物质晕参数。这是暗物质直接探测（间接：通过引力）最成熟的方法之一。

分解的 degeneracy 也很现实：提高恒星盘的质光比可以"吃掉"一部分内区暗物质贡献，这就是著名的 disk–halo degeneracy。突破口通常来自独立质量探针——气体盘的面密度可由 21 cm 直接换算、恒星动力学的垂直速度弥散可约束盘面密度、强透镜或外盘示踪体可钉住总质量。只有把这些锚点叠在同一条 $v(r)$ 上，晕参数才从"可调函数"变成可 falsify 的物理描述。

暗物质的粒子候选

暗物质的粒子本质是现代物理学最重大的未解之谜。主要候选者：

WIMP（弱相互作用大质量粒子）

质量在 ∼1\sim 1–$1000$ GeV，与标准模型粒子有弱力量级的相互作用。WIMP 在早期宇宙热平衡中产生，冻出后留下与观测一致的丰度（"WIMP 奇迹"，thermal relic mechanism）。

探测途径：

  • 直接探测：WIMP 与探测器中的原子核弹性散射；目前最强约束来自 LZ 实验（LUX-ZEPLIN）：首批结果给出自旋无关 WIMP-核子截面 <9.2×10−48< 9.2 \times 10^{-48} cm²（质量 ∼36\sim 36 GeV 处，Aalbers et al. 2023），2024 年 4.2 吨·年曝光进一步收紧至约 2×10−482 \times 10^{-48} cm²
  • 间接探测：WIMP 湮灭产生 γ\gamma 射线、正电子、中微子，Fermi-LAT、AMS-02 等设备寻找信号
  • 对撞机产生：LHC 寻找超对称粒子（自然的 WIMP 候选）

目前 WIMP 的所有探测尝试均以上限告终，大质量 WIMP（$> 100$ GeV）参数空间已被严重压缩，但低质量 WIMP（$< 10$ GeV）仍有大量未探索空间。

轴子（Axion）

轴子是由 Peccei-Quinn 对称性自发破缺产生的极轻标量粒子（质量 ∼10−6\sim 10^{-6}–10−310^{-3} eV），最初被提出用于解决强 CP 问题，但同时也是优秀的暗物质候选者。轴子冷暗物质通过相干振荡而非热产生。探测依赖轴子-光子耦合（a→γγa \to \gamma\gamma）：ADMX、HAYSTAC 等腔体实验已排除部分参数空间；SKA 等射电望远镜期望探测宇宙磁场中轴子转化产生的射电信号。

原初黑洞（PBH）

宇宙早期密度涨落可在暴胀结束后的辐射主导期坍缩形成原初黑洞。若其质量在 10−1610^{-16}–10−11M⊙10^{-11} M_\odot 范围内，引力微透镜约束允许 PBH 构成大部分暗物质（Carr et al. 2021）。LIGO/Virgo/KAGRA 观测到的部分黑洞并合事件（质量 ∼10\sim 10–100M⊙100 M_\odot）是否包含 PBH 成分，目前有争议。

需要强调：粒子候选章节回答的是"晕由什么粒子组成"，而旋转曲线与 NFW 拟合回答的是"晕的引力势长什么样"。二者相关但不可互相替代——即使实验室长期未直接探测到某种粒子，晕动力学证据依然独立成立；反之，某个粒子模型若无法生成观测所需的晕浓度与子结构统计，也会在天文侧被排除。 轴子与 WIMP 的实验语言不同：前者搜的是极弱的光子耦合与谐振腔调谐，后者搜的是核反冲与稀有事件率。天文侧则共同面对同一张晕质量表——无论粒子多轻或多重，都必须先生成观测所需的浓度、子结构统计与局域密度。实验室空白不能用来否定旋转曲线；天文拟合也不能代替粒子鉴别。

修正引力理论的挑战

面对暗物质长期未被粒子直接探测，部分物理学家转而考虑修改引力定律：

MOND（修正牛顿动力学）：Milgrom（1983）提出，当加速度 a<a0≈1.2×10−10a < a_0 \approx 1.2 \times 10^{-10} m/s²（约对应宇宙学常数量级）时，引力定律偏离牛顿形式。MOND 能成功预言许多旋涡星系旋转曲线的形状（Baryonic Tully-Fisher 关系），但在星系团尺度（子弹星系团）和宇宙学尺度（CMB 各向异性、大尺度结构）均面临严重困难，需要引入额外的"中微子暗物质"方能拟合。

MOND 已发展出相对论性版本（TeVeS、RMOND 等），但均无法同时解释所有观测证据，目前不是主流理论框架。 把 MOND 与晕动力学并置时，关键不在口号而在尺度迁移：星系盘尺度上拟合得好，并不自动保证在碰撞星系团、弱透镜质量图与原初功率谱上同样成立。子弹星系团的透镜质量远超可见重子，是这类修正必须正面回答的反例；若再引入额外碰撞物质来补洞，理论负担就会重新接近引入暗晕本身。

旋转曲线张力：核-尖峰问题

NFW 轮廓预言的中心尖峰（ρ∝r−1\rho \propto r^{-1}）在矮不规则星系的旋转曲线中往往找不到——观测更支持平坦的核（core）轮廓。这一"核-尖峰问题"是 Λ\LambdaCDM 在小尺度上最持久的挑战之一。

超新星反馈可将中心暗物质从尖峰"加热"为核：恒星形成过程中多次超新星爆发驱动气体剧烈流动，通过引力加热使暗物质从中心散开（Pontzen & Governato 2014）。这一机制在高效恒星形成的矮星系（M∗>106M⊙M_* > 10^6 M_\odot）中效果明显，但对超暗弱矮星系（M∗<105M⊙M_* < 10^5 M_\odot）似乎不足以解释观测。

因此核—尖峰之争本质上是晕动力学与重子循环的耦合问题，而不是单纯"改一条密度幂律"。若反馈能反复抬升中心势阱中的气体，暗物质轨道被加热、相空间密度被摊平，cusp 可变 soft core；若恒星形成被再电离紫外背景过早掐断，反馈次数不够，尖峰应更易保留。把不同恒星质量的矮星系排成序列，就是在对这一耦合做差分检验——旋转曲线形状随 M∗M_* 的系统变化，比单个天体的拟合残差更有判决力。

暗物质晕的结构与合并

在 Λ\LambdaCDM 框架下，暗物质晕通过层级并合从小尺度增长到大尺度。每个大晕含有大量子晕（subhalos）——较小的暗物质团块，是矮卫星星系的居所。$N$ 体数值模拟（Via Lactea II 模拟，Diemand et al. 2008；以及 Aquarius 模拟，Springel et al. 2008）显示，银河系质量晕中有约 10410^4–10510^5 个可分辨的子晕，但观测到的明亮卫星矮星系只有约 50–60 个——这是"缺失卫星问题"的核心。

晕的增长时间由特征时间：

tdyn=3π16Gρˉ≈1.4 Gyr(Δc200)−1/2(H(z)H0)−1t_\text{dyn} = \sqrt{\frac{3\pi}{16 G \bar{\rho}}} \approx 1.4 \text{ Gyr} \left(\frac{\Delta_c}{200}\right)^{-1/2} \left(\frac{H(z)}{H_0}\right)^{-1}

决定，其中 Δc\Delta_c 是相对于临界密度的超密度因子。

子晕的命运还受潮汐剥离与动力学摩擦支配：落入主晕深处的子结构会丢失外层质量，留下更致密的核心；有的几乎被撕成恒星流，只在相位空间留下褶皱。因此"模拟里数得出的子晕数目"与"天上看得见的卫星星系数目"之间，隔着选择函数、潮汐生存与能否形成恒星三道滤网。晕不是光滑背景，而是不断并合、剥离、再沉积的活系统——银河系晕中的巨流与壳层结构，正是这份并合史的化石记录。

跨域连接

  • 矮星系：判决核与尖峰要在暗物质主导、重子干扰最小的系统里做，因此矮不规则星系的中性氢旋转曲线是主战场。推论是超新星"加热"这一解释有一条可检验的边界：恒星质量太低时超新星数目不够，尖峰应当保留——这把调参问题变成了阈值预言。
  • 收敛：把标准密度轮廓积起来会发现总质量随半径对数发散，因此"晕的质量"必须先约定一个边界才有意义。推论是不同文献引用的晕质量常常差得不少，而差别往往只在边界约定上，不在测量上。
  • 万有引力与开普勒定律：开普勒式下降是"质量全在里面"的直接后果；观测到旋转速度近乎不变，等价于半径以内的质量仍随半径线性增加。推论只用到牛顿引力与圆轨道：要么存在看不见的质量，要么引力定律在极低加速度下变了。
  • 分治：多体模拟的可行性来自把两两求和换成树或网格近似——远处一团粒子用它的整体力矩代替。推论是"可分辨子晕的数目"是算法量而不是物理量：换分辨率、换判据，数目就变，因此缺失卫星的比较必须在同一分辨率下做。
  • 稀有气体：直接探测把最重的稳定惰性气体当靶，理由是化学的：闭壳结构使氙几乎不参与反应，因此能提纯到极高纯度并长期保持洁净，而本底控制正是这类实验的成败所在。推论是靶的原子量、核自旋与化学惰性决定它对哪类相互作用敏感。

晕的装配史还留下可观测的各向异性：晚近大并合会抬高速度椭球的径向各向异性，并在外晕堆出壳层与巨流；更平静的吸积则倾向更光滑、更近各向同性的外区。因此，把旋转曲线、晕星运动学与子结构普查放在同一势模型里联合拟合，是当前银河系与近邻星系研究的标准做法——目标不只是一个质量数字，而是一份带时间标签的质量装配年表。

为什么暗物质晕如此重要

暗物质晕是星系形成的骨架。没有暗物质晕的额外引力势阱，宇宙中弥散的重子物质（气体）无法克服膨胀宇宙的稀释效应而冷却坍缩成星系——我们所知的恒星、行星乃至生命都不会存在。暗物质虽然不可见，却是宇宙演化最深处的推手。 次晕动力学与卫星星系分布：CDM 模拟预测银河系质量暗物质晕中应有约 10310^3–10410^4 个次晕（subhalos），而迄今已发现的银河系卫星星系约 60 个（主要来自 SDSS 和 DES 巡天）。即使考虑观测限制（低面亮度系统未被发现），模拟预测的次晕数量仍显著多于观测——这是"消失卫星星系问题"（missing satellite problem）。一个可能的解决方案是次晕大多是暗的（无恒星形成），因为紫外背景辐射在宇宙再电离后阻止了小次晕中的气体冷却和恒星形成（Bullock et al. 2000）。Vera Rubin 天文台（LSST，预计 2025 年后大量科学数据）和罗曼空间望远镜（Nancy Grace Roman Space Telescope）将把已发现的矮星系/次晕数量提高到数百个，这是检验 CDM 次晕预测的决定性机会。 晕尺度上的可观测量：除实验室粒子搜寻外，晕本身提供一组与候选粒子弱相关的天文可观测量——外盘旋转速度、卫星轨道衰减、恒星流缺口、弱透镜质量图与逃逸速度曲线。它们共同约束的是势阱形状与质量装配史；粒子模型必须先通过这些晕动力学门槛，再谈截面或耦合。把"晕质量、浓度、子结构丰度"测准，本身就是缩小暗物质理论空间的独立路径。 暗物质在宇宙大尺度结构形成中的核心作用：宇宙微波背景（CMB）涨落在暗物质晕的引力作用下增长，形成了今天观测到的宇宙网状结构（cosmic web）。没有暗物质，重子物质无法在宇宙辐射主导时期（复合前）聚集（辐射压力阻止重子聚集），宇宙今天将没有星系和恒星。从 CMB 到今天的宇宙大尺度结构的演化，在 Λ\LambdaCDM 框架下，用 N 体数值模拟（如 Millennium Simulation，Springel et al. 2005，Nature；IllustrisTNG，Pillepich et al. 2018，Monthly Notices of the Royal Astronomical Society）能够高度准确地重现观测到的星系分布统计——这是 CDM 最强的支持之一。 综述性展望：该领域在未来 10 年将经历观测上的飞跃——Vera Rubin 天文台（LSST）、Euclid 卫星、罗曼空间望远镜（Nancy Grace Roman Space Telescope）和平方公里阵列射电望远镜（SKA）将分别在弱引力透镜、光谱测量、近红外成像和 HI 21cm 线巡天上提供数量级的数据量增加，推动对这一主题的统计理解进入新阶段。 暗物质晕关键物理量速查：

  • NFW 轮廓特征半径 rsr_s：银河系约 20 kpc；M87 约 200 kpc
  • 银河系暗物质晕质量：约 $0.9$–1.5×1012M⊙1.5 \times 10^{12} M_\odot（Eadie & Juric 2019）
  • 银河系太阳附近暗物质局部密度：约 $0.4$ GeV/cm3^3（Bovy & Tremaine 2012）
  • WIMP 直接探测约束（LZ 首批结果，2023）：σSI<9.2×10−48\sigma_{SI} < 9.2 \times 10^{-48} cm2^2 at 36 GeV（2024 年 4.2 吨·年结果约 2×10−482 \times 10^{-48} cm2^2）
  • 轴子暗物质质量范围：约 10−610^{-6}–10−310^{-3} eV（ADMX 等实验探测目标）

参考文献

  • Rubin, V.C. et al. (1980). Rotational properties of 21 Sc galaxies. Astrophysical Journal, 238, 471–487.
  • Navarro, J.F., Frenk, C.S. & White, S.D.M. (1996). The structure of cold dark matter halos. Astrophysical Journal, 462, 563–575.
  • Akerib, D.S. et al. (LZ Collaboration) (2023). First dark matter search results from the LUX-ZEPLIN experiment. Physical Review Letters, 131, 041002.
  • Pontzen, A. & Governato, F. (2014). Cold dark matter heats up. Nature, 506, 171–178.
  • de Blok, W.J.G. (2010). The core-cusp problem. Advances in Astronomy, 2010, 789293.
  • Milgrom, M. (1983). A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. Astrophysical Journal, 270, 365–370.

正文中引用本条 · 10 · 其中 2 条来自其他领域

宇宙物理暗物质暗物质直接探测的"沉默"