4.3 层次性并合模型
4.3 层次性并合模型
本节摘要：冷暗物质宇宙的结构形成是"自下而上"的：小晕先坍缩，再通过并合逐级组装成大星系。本节讲述层次性并合模型的物理内核、晕质量函数与并合率的定量工具，并讨论"缺卫星问题""过大核球问题"等模型面临的挑战，以及观测对模型的检验。
自下而上的组装逻辑
若暗物质是"冷"的（粒子运动远低于光速），密度扰动在物质主导期之前就可以在亚星系尺度上增长，小结构先形成，大结构由小结构并合而来——这是层次性并合模型最核心的逻辑，与暗物质是"热"（相对论性粒子、会冲淡小尺度扰动）的情形截然相反。后者预言大结构先形成、星系后碎裂，早已被观测排除。
数学上，给定初始涨落谱，任一质量的暗物质晕的形成时间由该尺度涨落超过坍缩阈值的时间决定。小质量晕对应小尺度、高振幅涨落，更早形成；大质量晕对应大尺度、低振幅涨落，更晚形成。于是一个 10¹² 太阳质量的银河系级晕，其质量的大半来自数十上百个 10⁸-10¹⁰ 太阳质量的前身晕，这些前身在红移 3-10 之间先后形成，再不断并合。暗物质的这一组装过程被 N 体模拟完美再现：晕质量函数的形状、并合树的分支结构、晕内部密度轮廓，全部由初始条件唯一决定。
图：层次性并合的时间线
晕质量函数与并合率
层次性并合的定量描述靠两个统计量。晕质量函数 n(M) 给出每单位体积内质量在 M 到 M+dM 之间的暗物质晕数目。Sheth-Tormen 形式给出：低质量端为幂律 dn/dM ∝ M^(-1.9)，高质量端被指数截断，特征质量 M* 随时间增大——因此高红移时大质量晕稀少，今天则普遍。并合率描述每个晕单位时间吸积的质量：质量增长率大致与质量成正比（dM/dt ∝ M），意味着大星系比小星系长得更快，形成所谓的"整体增长"（downsizing 现象：大质量星系的恒星形成反而在更早红移完成）。
并合率随环境显著变化：星系团中心的大星系几乎每天都有卫星落入，而场星系的环境要平静得多。宇宙学模拟给出的并合率解析式是星系形成半解析模型（SAM）的核心输入。结合恒星形成与反馈子网格模型，SAM 可以从初始条件出发预测今天的星系质量函数、形态分布与恒星形成率——这是第 7 章模拟内容的先声。
并合如何塑造星系形态
层次性并合模型对形态分布给出了可检验预言。主并合（质量比小于 4）把气体盘搅碎、驱动星爆，是椭圆星系形成的主要通道；而吸积角动量大的冷气体形成新盘，则是盘星系在两次并合之间"自我修复"的机制。两条通道的竞争结果由红移与质量决定：高红移、高质量环境主并合频率高，故早型星系比例高；低红移、低密度环境盘星系存续比例高。
银河系自身的结构就写在层次性历史的账本上：银晕中的恒星流、球状星团年龄双峰、核球中的盒状轨道，都是多次次并合与一两次主并合留下的考古学证据。模拟显示，银河系在过去的百亿年里经历了约十次质量比大于十比一的次并合，以及数次更小的吸积事件。
模型面临的挑战
层次性并合模型并非无懈可击，三大问题驱动着理论不断修正。缺卫星问题：模拟预言银河系周围应有上千个暗物质子晕，而观测到的卫星星系只有几十个——解药之一是"反馈压制"：小晕难以留住气体形成恒星，只留下暗物质，因此卫星星系比卫星晕少得多，这需要超新星反馈与再电离加热共同作用。过大核球问题：纯 ΛCDM 模拟预言盘星系频繁主并合会制造过大核球，与观测的薄盘星系核球占比不符；回答涉及并合时机的统计与反馈抑制星爆的效率。核心-尖点与多样性问题：小质量晕的中心密度观测普遍低于模拟预言，除了暗物质自相互作用，反馈驱动重子振荡也可能是罪魁祸首。这些张力并不推翻层次性框架，反而把"重子物理"推到前台——星系形成研究从此从"算暗物质"转向"算气体与恒星"。
Python 算例：并合树的统计特征
import numpy as np def merger_rate(M_main, z): """ 简化并合率公式：单位 Gyr 内并入质量大于 M_frac*M_main 的伴星系次数。 量级依据 Fakhouri & Ma 的模拟拟合式（示意性参数）。 """ M12 = M_main / 1e12 return 0.010 * M12**0.13 * (1 + z)**2.2 * 4.0 # 不同红移下的主并合（质量比>0.3）频率 for z in [0, 0.5, 1, 2, 4]: M = 1e12 # 银河系级晕 rate = merger_rate(M, z) * (1 + z) * 0.01 print(f"z={z:4.1f} 主并合特征频率 ≈ {rate*1e2:6.2f} 次/Gyr") # 平均而言：z=0 的银河系级晕平均约 3-5 Gyr 经历一次主并合， # 而次并合（质量比>0.1）频率高一个量级，这与银河系盘存活至今的观测一致
从模型到观测的闭环
层次性并合模型的检验闭环如今已相当完整：CMB 给出初始条件，N 体模拟给出纯暗物质的预言，观测给出星系性质的约束，而连接二者的正是并合树与重子物理。当 JWST 在高红移看到大量"出乎意料地大"的早型星系时，模型被迫调整并合时标与恒星形成效率；当暗能量巡天测出星系团的丰度时，晕质量函数的形状反过来约束 σ₈。模型与观测的每一次"摩擦"，都是理论前进的动力。下一节将把目光从暗物质转向恒星与元素，讨论星系如何把气体炼成星光与金属。