暗物质(Dark Matter)是宇宙中一种不发光、不吸收也不发射电磁辐射的神秘物质成分。它占据宇宙总质量-能量密度的约 27%,是可见普通物质的 5 倍以上,但其本质至今仍是现代物理学最大的未解之谜之一。暗物质的存在完全通过引力效应推断——它"看不见",但它的引力"触手"塑造了宇宙从大爆炸至今的每一层结构。
暗物质的认识经历了从猜想到实证的百年历程:
时间
人物
贡献
方法
1933 年
弗里茨·兹威基(Fritz Zwicky)
发现后发座星系团星系速度弥散过大,推断存在"暗物质"(dunkle Materie)
维里定理计算质量
1970 年代
维拉·鲁宾(Vera Rubin)
发现旋涡星系(如 Andromeda)旋转曲线平坦,可见质量远不足以维持
光谱多普勒测量
2006 年
子弹星系团观测
碰撞中可见气体与引力质量分离,证明暗物质与普通物质本质不同
X 射线 + 引力透镜
兹威基在后发座星系团中测量了星系的速度弥散为 σ ≈ 770 km/s \sigma \approx 770\ \text{km/s}σ ≈ 770 km/s 。利用维里定理,星系团的平均引力势能 ⟨ U ⟩ \langle U \rangle⟨ U ⟩ 与动能 ⟨ K ⟩ \langle K \rangle⟨ K ⟩ 满足:
2 ⟨ K ⟩ + ⟨ U ⟩ = 0 2 \langle K \rangle + \langle U \rangle = 0
2 ⟨ K ⟩ + ⟨ U ⟩ = 0
其中动能 ⟨ K ⟩ = 1 2 M σ 2 \langle K \rangle = \frac{1}{2} M \sigma^2⟨ K ⟩ = 2 1 M σ 2 ,引力势能 ⟨ U ⟩ = − 3 5 G M 2 R \langle U \rangle = -\frac{3}{5} \frac{GM^2}{R}⟨ U ⟩ = − 5 3 R G M 2 。代入 σ = 770 km/s \sigma = 770\ \text{km/s}σ = 770 km/s 和 R ≈ 1.5 × 10 22 m R \approx 1.5 \times 10^{22}\ \text{m}R ≈ 1.5 × 1 0 22 m (约 1.5 Mpc)进行数值计算:
M = 5 R σ 2 3 G = 5 × ( 1.5 × 10 22 ) × ( 7.7 × 10 5 ) 2 3 × 6.67 × 10 − 11 M = \frac{5 R \sigma^2}{3 G}
= \frac{5 \times (1.5 \times 10^{22}) \times (7.7 \times 10^5)^2}{3 \times 6.67 \times 10^{-11}}
M = 3 G 5 R σ 2 = 3 × 6.67 × 1 0 − 11 5 × ( 1.5 × 1 0 22 ) × ( 7.7 × 1 0 5 ) 2
M ≈ 5 × 1.5 × 10 22 × 5.93 × 10 11 2.00 × 10 − 10 ≈ 2.2 × 10 45 kg ≈ 1.1 × 10 15 M ⊙ M \approx \frac{5 \times 1.5 \times 10^{22} \times 5.93 \times 10^{11}}{2.00 \times 10^{-10}}
\approx 2.2 \times 10^{45}\ \text{kg} \approx 1.1 \times 10^{15}\ M_\odot
M ≈ 2.00 × 1 0 − 10 5 × 1.5 × 1 0 22 × 5.93 × 1 0 11 ≈ 2.2 × 1 0 45 kg ≈ 1.1 × 1 0 15 M ⊙
而从星系团中可见星系的总发光质量仅约 10 13 M ⊙ 10^{13}\ M_\odot1 0 13 M ⊙ 。可见质量与动力学质量差了约 100 倍 ——这就是暗物质的第一个定量证据。
暗物质的证据并非单一来源,而是来自完全独立的天文观测方法,相互印证构成完整证据链。
旋涡星系中,恒星和气体绕星系中心旋转。根据开普勒定律,如果质量集中在中心,轨道速度为 v ( r ) ∝ 1 / r v(r) \propto 1/\sqrt{r}v ( r ) ∝ 1/ r 。然而实际观测发现旋转曲线在远处保持平坦甚至上升。
具体数据:M33 三角座星系的旋转曲线
距中心距离 r rr (kpc)
可见物质预测 v vv (km/s)
实际观测 v vv (km/s)
需要的暗物质质量比
2
80
85
1.1×
5
50
95
3.6×
8
32
100
9.8×
12
21
105
25×
16
16
105
43×
旋涡星系的旋转曲线可以用以下经验公式描述(由 Persic, Salucci & Stel 1996 提出):
v ( r ) = v 0 [ 1 + ( r r opt ) β ] − 1 2 β v(r) = v_0 \left[ 1 + \left( \frac{r}{r_{\text{opt}}} \right)^\beta \right]^{-\frac{1}{2\beta}}
v ( r ) = v 0 [ 1 + ( r opt r ) β ] − 2 β 1
其中 v 0 v_0v 0 是渐近平坦速度,r opt r_{\text{opt}}r opt 是光学半径,β ≈ 0.5 \beta \approx 0.5β ≈ 0.5 。以 M33 为例,v 0 ≈ 105 km/s v_0 \approx 105\ \text{km/s}v 0 ≈ 105 km/s ,r opt ≈ 8 kpc r_{\text{opt}} \approx 8\ \text{kpc}r opt ≈ 8 kpc 。在 r = 16 kpc r = 16\ \text{kpc}r = 16 kpc 处,可见物质的预测速度仅为 16 km/s 16\ \text{km/s}16 km/s ,而实际观测为 105 km/s 105\ \text{km/s}105 km/s ——这意味着 16 kpc 处 98% 以上的质量来自暗物质。
维拉·鲁宾在 1970-1980 年代对超过 60 个旋涡星系进行了系统测量,发表了标志性论文 "Rotation of the Andromeda Nebula from a Spectroscopic Survey of Emission Regions" (1970),首次系统性地展示了旋转曲线的平坦性。
爱因斯坦广义相对论预言大质量天体会弯曲周围时空,使背景星系的光线发生偏折。引力透镜效应既是暗物质的证据,也可用来绘制暗物质分布图。
引力偏折角公式:
α ^ = 4 G M c 2 ξ \hat{\alpha} = \frac{4GM}{c^2 \xi}
α ^ = c 2 ξ 4 GM
其中 ξ \xiξ 是碰撞参数(光线到透镜质心的距离),M MM 是透镜质量。对于一个典型星系团,如 Abell 1689(z = 0.18 z = 0.18z = 0.18 ),强引力透镜的多重像系统表明其总质量为 M ≈ 1 × 10 15 M ⊙ M \approx 1 \times 10^{15}\ M_\odotM ≈ 1 × 1 0 15 M ⊙ ,而其可见的星系和热气体(X 射线)质量仅约 2 × 10 14 M ⊙ 2 \times 10^{14}\ M_\odot2 × 1 0 14 M ⊙ ,即约 80% 的质量为暗物质。
弱透镜剪切(Weak Lensing)的定量测量:
弱透镜无法直接看到多重像,但可以统计测量背景星系形状的微弱相干扭曲——"剪切"(shear)γ \gammaγ 。剪切信号与表面质量密度 Σ \SigmaΣ 的关系通过临界密度 Σ crit \Sigma_{\text{crit}}Σ crit 连接:
Σ crit = c 2 4 π G D s D l D l s \Sigma_{\text{crit}} = \frac{c^2}{4\pi G} \frac{D_s}{D_l D_{ls}}
Σ crit = 4 π G c 2 D l D l s D s
其中 D s D_sD s 、D l D_lD l 、D l s D_{ls}D l s 分别是源-观测者、透镜-观测者和源-透镜的角直径距离。对于红移 z l = 0.3 z_l = 0.3z l = 0.3 的透镜和 z s = 1.0 z_s = 1.0z s = 1.0 的源,Σ crit ≈ 3 × 10 9 M ⊙ / kpc 2 \Sigma_{\text{crit}} \approx 3 \times 10^9\ M_\odot/\text{kpc}^2Σ crit ≈ 3 × 1 0 9 M ⊙ / kpc 2 。若沿视线方向测量的剪切 ∣ γ ∣ ≈ 0.05 |\gamma| \approx 0.05∣ γ ∣ ≈ 0.05 ,对应表面密度 Σ ≈ 1.5 × 10 8 M ⊙ / kpc 2 \Sigma \approx 1.5 \times 10^8\ M_\odot/\text{kpc}^2Σ ≈ 1.5 × 1 0 8 M ⊙ / kpc 2 。
这是暗物质存在的最直观证据 。2006 年,钱德拉 X 射线望远镜和哈勃空间望远镜联合观测了正在碰撞的两个星系团。
观测结果:
碰撞方向 →
可见物质(X射线气体): 引力质量(引力透镜):
┌─────────────────┐ ★═══════★
│ Hottest X-ray │ ║ ║
│ gas (T ≈ │ ║ DM ║
│ 1.5 × 10^8 K) │ ║ ║
└─────────────────┘ ★═══════★
(被气体曳力减速留在中心) (几乎无碰撞,继续运动)
关键参数:
参数
数值
说明
碰撞速度
≈ 4500 km/s \approx 4500\ \text{km/s}≈ 4500 km/s
目前已知最快的星系团碰撞之一
X射线气体温度
T ≈ 1.5 × 10 8 K T \approx 1.5 \times 10^8\ \text{K}T ≈ 1.5 × 1 0 8 K (约 13 keV)
碰撞使气体被加热
可见气体质量分数
总重子质量的 ≈ 90%
星系团中大部分普通物质是热气体
暗物质与气体偏移量
≈ 200 kpc \approx 200\ \text{kpc}≈ 200 kpc
清晰显示了暗物质与普通物质的分离
子弹星系团排除了两种替代解释:
MOND(修正牛顿动力学) :无法解释两个质量峰分离的现象,因为 MOND 的引力修正依赖于总质量分布
纯粹重子暗物质 :如果暗物质是普通不发光的物质(如 MACHO),它在碰撞中也会与气体相互作用减速,不会分离
CMB 的功率谱包含精确的宇宙成分信息。由于重子与光子紧耦合,而暗物质不与光子相互作用,两者在早期宇宙中的演化完全不同。
Planck 2018 测量的宇宙成分参数:
参数
符号
数值
物理含义
暗物质密度
Ω c h 2 \Omega_c h^2Ω c h 2
0.120 ± 0.001 0.120 \pm 0.0010.120 ± 0.001
冷暗物质的能量密度参数
重子密度
Ω b h 2 \Omega_b h^2Ω b h 2
0.0224 ± 0.0001 0.0224 \pm 0.00010.0224 ± 0.0001
普通物质密度
暗能量密度
Ω Λ \Omega_\LambdaΩ Λ
0.6847 ± 0.0073 0.6847 \pm 0.00730.6847 ± 0.0073
暗能量占比
暗物质占比
Ω c / Ω m \Omega_c / \Omega_mΩ c / Ω m
≈ 0.84 \approx 0.84≈ 0.84
暗物质占总物质的 84%
标量谱指数
n s n_sn s
0.9649 ± 0.0042 0.9649 \pm 0.00420.9649 ± 0.0042
原初扰动谱斜率
哈勃常数
H 0 H_0H 0
67.4 ± 0.5 km/s/Mpc 67.4 \pm 0.5\ \text{km/s/Mpc}67.4 ± 0.5 km/s/Mpc
当前宇宙膨胀率
CMB 功率谱中的声学峰位置特别敏感于 Ω m h 2 \Omega_m h^2Ω m h 2 。第一个峰(多极数 ℓ ≈ 200 \ell \approx 200ℓ ≈ 200 )的位置由角直径距离决定,而后续峰的相对高度比是重子与暗物质密度比的函数。Planck 卫星以 0.8% 的精度测定了 Ω c h 2 = 0.120 ± 0.001 \Omega_c h^2 = 0.120 \pm 0.001Ω c h 2 = 0.120 ± 0.001 ——这是暗物质存在的最精确的大尺度证据。
暗物质还影响了宇宙中星系的分布模式。数值模拟(如 Millennium、IllustrisTNG)表明,宇宙的大尺度结构——由星系构成的"宇宙网"——只能在暗物质主导的宇宙中形成:
特征
有暗物质的宇宙
无暗物质的纯重子宇宙
结构形成时间
红移 z ∼ 20 z \sim 20z ∼ 20 开始形成
推迟到 z ∼ 6 z \sim 6z ∼ 6
宇宙网
清晰丝状结构
弥散不均匀
星系数量
约 2 万亿个(哈勃深场估计)
少 10-100 倍
空洞大小
典型 30-100 Mpc
缺乏大尺度空洞
斯隆数字化巡天(SDSS)和 2dF 巡天绘制的星系分布图与 Λ \LambdaΛ CDM 模拟的宇宙网高度一致,提供了暗物质存在的独立验证。
精确的宇宙成分饼图如下:
成分
占比
已知程度
暗能量(Λ \LambdaΛ )
68.3 % 68.3\%68.3%
完全未知
暗物质
26.8 % 26.8\%26.8%
知其引力效应,本质未知
普通物质(重子)
4.9 % 4.9\%4.9%
已知(恒星、气体、尘埃)
中微子
0.1 % − 1 % 0.1\% - 1\%0.1% − 1%
已知部分性质
光子
0.005 % 0.005\%0.005%
已知(CMB)
暗物质与普通物质的总质量比:
Ω c Ω b = 0.268 0.049 ≈ 5.5 \frac{\Omega_c}{\Omega_b} = \frac{0.268}{0.049} \approx 5.5
Ω b Ω c = 0.049 0.268 ≈ 5.5
这意味着宇宙中不可见的质量是可见质量的 5.5 倍 。在一个典型星系中,总质量为 10 12 M ⊙ 10^{12}\ M_\odot1 0 12 M ⊙ 的暗物质晕包含了约 10 11 M ⊙ 10^{11}\ M_\odot1 0 11 M ⊙ 的可见恒星和气体。但这种比例随宇宙时间演化:在早期宇宙(z ≳ 10 z \gtrsim 10z ≳ 10 ),重子物质占比相对较高;随着恒星形成消耗气体,今天的重子占比持续下降。
重子缺失问题: 在大尺度上观测到的重子数量仅为宇宙学预测的约 50%。这些"缺失的重子"可能在星系间的温热气体中(WHIM,温热星系际介质),温度在 10 5 − 10 7 K 10^5-10^7\ \text{K}1 0 5 − 1 0 7 K ,难以观测。
暗物质的本质是粒子物理的核心问题。以下是最主要的候选粒子:
WIMP 是暗物质最主流、最被广泛研究的候选者。其重要原因在于"WIMP 奇迹":如果早期宇宙中存在一种质量在 GeV-TeV 量级、相互作用强度约为弱相互作用的粒子,其热产生的剩余丰度恰好与观测到的暗物质密度相符。
"WIMP 奇迹"的定量推导:
暗物质的剩余丰度由以下公式给出:
Ω c h 2 ≈ 3 × 10 − 27 cm 3 / s ⟨ σ v ⟩ \Omega_c h^2 \approx \frac{3 \times 10^{-27}\ \text{cm}^3/\text{s}}{\langle \sigma v \rangle}
Ω c h 2 ≈ ⟨ σ v ⟩ 3 × 1 0 − 27 cm 3 / s
其中 ⟨ σ v ⟩ \langle \sigma v \rangle⟨ σ v ⟩ 是暗物质湮灭截面与相对速度的热平均。代入 ⟨ σ v ⟩ ≈ 3 × 10 − 26 cm 3 / s \langle \sigma v \rangle \approx 3 \times 10^{-26}\ \text{cm}^3/\text{s}⟨ σ v ⟩ ≈ 3 × 1 0 − 26 cm 3 / s (弱相互作用的典型值),可得 Ω c h 2 ≈ 0.1 \Omega_c h^2 \approx 0.1Ω c h 2 ≈ 0.1 ,与观测值 0.12 0.120.12 量级一致。这不是巧合,而是暗示暗物质可能就是一种热产生的 WIMP。
WIMP 的属性对比:
属性
数值/范围
说明
质量
10 GeV − 10 TeV 10\ \text{GeV} - 10\ \text{TeV}10 GeV − 10 TeV
比质子重 10-10000 倍
相互作用
弱核力 + 引力
不参与电磁力和强相互作用
自旋
通常假设为 1/2(费米子)
部分模型如超对称中性子
速度
非相对论性(冷暗物质)
v ∼ 200 km/s v \sim 200\ \text{km/s}v ∼ 200 km/s
候选理论
超对称(neutralino)
SUSY 最轻超对称粒子
天然性指数
无需精细调节
WIMP 奇迹提供自然解释
轴子最初为解决强 CP 问题而被佩切伊-奎因理论预言,后来被发现也是优质的暗物质候选。
轴子的关键特性:
质量极低:m a ∼ 10 − 6 − 10 − 3 eV m_a \sim 10^{-6} - 10^{-3}\ \text{eV}m a ∼ 1 0 − 6 − 1 0 − 3 eV
通过普里马科夫效应(Primakoff effect)与光子耦合
在强磁场中可能转化为光子
非热产生机制(真空对齐、拓扑缺陷衰变)
轴子探测实验对比:
实验
方法
质量范围
灵敏度
ADMX
谐振腔(强磁场)
1 − 100 μ eV 1 - 100\ \mu\text{eV}1 − 100 μ eV
可达 DFSZ 耦合
CAST
太阳轴子(磁铁指向太阳)
eV \text{eV}eV 级
已设上限
IAXO
下一代太阳轴子望远镜
eV \text{eV}eV 级
CAST 的 1000 倍
MADMAX
电介质盘
40 − 400 μ eV 40 - 400\ \mu\text{eV}40 − 400 μ eV
新方法
ALPS II
光穿墙实验(Laser)
meV \text{meV}meV 级
测试光子-轴子耦合
轴子暗物质的参数空间:
m a ≲ 10 − 2 eV , g a γ γ ≲ 10 − 10 GeV − 1 m_a \lesssim 10^{-2}\ \text{eV},\quad g_{a\gamma\gamma} \lesssim 10^{-10}\ \text{GeV}^{-1}
m a ≲ 1 0 − 2 eV , g aγ γ ≲ 1 0 − 10 GeV − 1
其中 g a γ γ g_{a\gamma\gamma}g aγ γ 是轴子-光子耦合常数。超出此范围会被恒星冷却观测(如 SN 1987A 的持续时间)排除。
惰性中微子(Sterile Neutrino)是一种假设的第四类中微子,不参与任何标准模型相互作用(包括弱相互作用),仅通过引力与其他物质耦合。质量在 keV 量级的惰性中微子是温热暗物质(Warm Dark Matter)候选。
惰性中微子作为暗物质的参数空间:
质量
产生机制
暗物质类型
观测约束
∼ 1 − 10 keV \sim 1 - 10\ \text{keV}∼ 1 − 10 keV
混合振荡(Dodelson-Widrow 机制)
温热暗物质
莱曼-α \alphaα 森林排除 m s < 10 keV m_s < 10\ \text{keV}m s < 10 keV
∼ 7 keV \sim 7\ \text{keV}∼ 7 keV
退耦后共振产生(MSW 效应增强)
温热暗物质
可能存在 3.5 keV X 射线线
∼ 100 keV − MeV \sim 100\ \text{keV} - \text{MeV}∼ 100 keV − MeV
非热产生
冷暗物质
结构形成限制更宽松
2014 年,两个独立团队(Boyarsky 等、Bulbul 等)在 XMM-Newton 对星系团和 M31 的观测中发现了约 3.5 keV 的未识别 X 射线发射线。如果来自惰性中微子衰变,则对应 m s ≈ 7 keV m_s \approx 7\ \text{keV}m s ≈ 7 keV 的暗物质粒子,寿命 τ ≈ 2 × 10 27 s \tau \approx 2 \times 10^{27}\ \text{s}τ ≈ 2 × 1 0 27 s (比宇宙年龄长 10 19 10^{19}1 0 19 倍)。后续 Hitomi 卫星对 Perseus 星系团的精细观测未确认该信号,因此目前争议中。
原初黑洞是宇宙早期由密度涨落直接坍缩形成的黑洞,质量范围远大于恒星质量黑洞。
PBH 质量
形成时间(大爆炸后)
当前可观测性
10 − 5 M ⊙ 10^{-5}\ M_\odot1 0 − 5 M ⊙
∼ 10 − 23 s \sim 10^{-23}\ \text{s}∼ 1 0 − 23 s
蒸发殆尽(霍金辐射)
10 − 1 M ⊙ 10^{-1}\ M_\odot1 0 − 1 M ⊙
∼ 1 s \sim 1\ \text{s}∼ 1 s
引力透镜微透镜效应
10 − 100 M ⊙ 10 - 100\ M_\odot10 − 100 M ⊙
∼ 10 s \sim 10\ \text{s}∼ 10 s
LIGO 引力波
10 3 − 10 5 M ⊙ 10^3 - 10^5\ M_\odot1 0 3 − 1 0 5 M ⊙
∼ 100 s \sim 100\ \text{s}∼ 100 s
星系核动力学加热
观测约束: 致密晕物体(MACHO)调查通过微透镜事件对 PBH 在 10 − 7 − 10 M ⊙ 10^{-7} - 10\ M_\odot1 0 − 7 − 10 M ⊙ 范围设置了严格上限。目前 PBH 只能占暗物质总量的 ≲ 10 % \lesssim 10\%≲ 10% (在大部分质量范围内)。LIGO 探测到的双黑洞并合事件部分可解释为 PBH,但无法解释全部暗物质。
PBH 作为暗物质的质量上限:
f PBH = Ω PBH Ω c ≲ 0.01 − 0.1 ( 大部分质量范围 ) f_{\text{PBH}} = \frac{\Omega_{\text{PBH}}}{\Omega_c} \lesssim 0.01 - 0.1 \quad (\text{大部分质量范围})
f PBH = Ω c Ω PBH ≲ 0.01 − 0.1 ( 大部分质量范围 )
候选者
质量
探测方式
理论优美性
实验支持
WIMP(中性子)
GeV-TeV
核反冲、湮灭、对撞机
★★★★★
尚未发现
轴子
μ eV \mu\text{eV}μ eV - meV \text{meV}meV
光子转换
★★★★☆
ADMX 在深入扫描
惰性中微子
keV
X 射线衰变线
★★★☆☆
3.5 keV 信号未确认
原初黑洞
M ⊙ M_\odotM ⊙ 级
引力透镜、引力波
★★★☆☆
约束严格,占比极小
直接探测的核心原理:暗物质粒子(假设 WIMP)飞过探测器并与原子核发生弹性碰撞,将微小的反冲能量传递给探测器。
WIMP-核散射率公式:
R = ρ 0 m χ 1 m N ∫ E th E max d σ d E R × η ( v min , t ) d E R R = \frac{\rho_0}{m_\chi} \frac{1}{m_N} \int_{E_{\text{th}}}^{E_{\text{max}}} \frac{d\sigma}{dE_R} \times \eta(v_{\text{min}}, t) \ dE_R
R = m χ ρ 0 m N 1 ∫ E th E max d E R d σ × η ( v min , t ) d E R
其中 ρ 0 ≈ 0.3 GeV/cm 3 \rho_0 \approx 0.3\ \text{GeV/cm}^3ρ 0 ≈ 0.3 GeV/cm 3 是局域暗物质密度,m χ m_\chim χ 是 WIMP 质量,m N m_Nm N 是靶核质量,E th E_{\text{th}}E th 是探测阈值。对于 m χ = 100 GeV m_\chi = 100\ \text{GeV}m χ = 100 GeV 的 WIMP 在氙靶(m N ≈ 131 amu m_N \approx 131\ \text{amu}m N ≈ 131 amu )中,预期事件率仅为每年每吨靶材 1-10 次。
主要实验对比:
实验
靶材料
质量
灵敏度下限
地点
LZ(LUX-ZEPLIN)
液氙 7 吨
美国 Homestake
1.4 × 10 − 48 cm 2 1.4 \times 10^{-48}\ \text{cm}^21.4 × 1 0 − 48 cm 2 (σ S I \sigma_{SI}σ S I )
地下 1500 m
XENONnT
液氙 5.9 吨
意大利 Gran Sasso
2.6 × 10 − 48 cm 2 2.6 \times 10^{-48}\ \text{cm}^22.6 × 1 0 − 48 cm 2
地下 1400 m
PandaX-4T
液氙 4 吨
中国锦屏
3.8 × 10 − 48 cm 2 3.8 \times 10^{-48}\ \text{cm}^23.8 × 1 0 − 48 cm 2
地下 2400 m
CDMSlite
锗
美国 Soudan
5 × 10 − 42 cm 2 5 \times 10^{-42}\ \text{cm}^25 × 1 0 − 42 cm 2
低质量 WIMP 探测
DarkSide-50
液氩 50 kg
意大利 Gran Sasso
1 × 10 − 44 cm 2 1 \times 10^{-44}\ \text{cm}^21 × 1 0 − 44 cm 2
脉冲形状甄别
LZ 实验在 2024 年发表了首批结果(220 天运行数据),在 m χ = 40 GeV m_\chi = 40\ \text{GeV}m χ = 40 GeV 处达到了自旋无关散射截面 σ S I < 1.4 × 10 − 48 cm 2 \sigma_{SI} < 1.4 \times 10^{-48}\ \text{cm}^2σ S I < 1.4 × 1 0 − 48 cm 2 的上限。这是目前对 WIMP 最灵敏的直接探测限制。
WIMP-核散射的质量依赖关系:
对于自旋无关(SI)散射,截面与靶核质量数 A AA 有关:
σ S I ∝ μ 2 A 2 \sigma_{SI} \propto \mu^2 A^2
σ S I ∝ μ 2 A 2
其中 μ = m χ m N m χ + m N \mu = \frac{m_\chi m_N}{m_\chi + m_N}μ = m χ + m N m χ m N 是约化质量。当 m χ ≈ m N m_\chi \approx m_Nm χ ≈ m N 时,匹配效率最佳。举例:40 GeV WIMP 与氩-40 核(A = 40 A=40A = 40 )的约化质量 μ = 40 × 37.4 40 + 37.4 ≈ 19.3 GeV \mu = \frac{40 \times 37.4}{40 + 37.4} \approx 19.3\ \text{GeV}μ = 40 + 37.4 40 × 37.4 ≈ 19.3 GeV ,灵敏度最高。这也是为什么多种靶材料(氙、氩、锗、硅)覆盖不同 WIMP 质量范围。
间接探测寻找暗物质粒子在宇宙中对撞湮灭或衰变产生的次级粒子:高能光子(γ \gammaγ )、中微子、正电子、反质子等。
暗物质湮灭的预期通量:
Φ γ ( E ) = 1 4 π ⟨ σ v ⟩ 2 m χ 2 d N γ d E × J ( Ω ) \Phi_\gamma(E) = \frac{1}{4\pi} \frac{\langle \sigma v \rangle}{2 m_\chi^2} \frac{dN_\gamma}{dE} \times J(\Omega)
Φ γ ( E ) = 4 π 1 2 m χ 2 ⟨ σ v ⟩ d E d N γ × J ( Ω )
其中 J ( Ω ) = ∫ los ρ 2 ( r ( s ) ) d s J(\Omega) = \int_{\text{los}} \rho^2(r(s)) \ dsJ ( Ω ) = ∫ los ρ 2 ( r ( s )) d s 是 J JJ 因子,沿视线方向对暗物质密度平方的积分。银河系中心(GC)具有最大的 J JJ 因子,是最有希望的湮灭信号来源。
主要间接探测实验:
实验/仪器
探测粒子
能量范围
关键观测
状态
Fermi-LAT
γ \gammaγ 射线
20 MeV - 300 GeV
银河系中心 γ \gammaγ 超量
争议中
AMS-02
正电子、反质子
GeV-TeV
正电子分数超量
确认超量,来源不确定
H.E.S.S.
TeV γ \gammaγ 射线
100 GeV - 100 TeV
GC 最严格上限
仍在观测
IceCube
高能中微子
TeV - PeV
太阳方向无中微子信号
设严格上限
DAMPE
电子、γ \gammaγ
5 GeV - 10 TeV
1.4 TeV 电子谱拐折
官方确认
银河系中心 GeV 超量: Fermi-LAT 在 2009 年发现银河系中心方向存在一个额外 γ \gammaγ 射线发射,空间分布呈球对称,能量谱在约 2 GeV 处有峰值。这可以被解释为 m χ ≈ 35 − 50 GeV m_\chi \approx 35 - 50\ \text{GeV}m χ ≈ 35 − 50 GeV 的 WIMP 湮灭到 b b ˉ b\bar{b}b b ˉ 末态信号。但后续分析表明,至少部分信号可能来自银河系中心的未解析脉冲星——这是一个仍未解决的争议。
大型强子对撞机(LHC)可能通过高能质子对撞产生暗物质粒子。由于暗物质不参与电磁相互作用,其在探测器中表现为缺失的能量-动量——即"缺失横动量"(MET, Missing Transverse Momentum)。
LHC 对暗物质的主要探测通道:
p p → χ χ ˉ + ISR(初始态辐射) pp \to \chi\bar{\chi} + \text{ISR(初始态辐射)}
pp → χ χ ˉ + ISR (初始态辐射)
过程举例:单喷注 + MET:
p p → χ χ ˉ + jet (monojet) pp \to \chi\bar{\chi} + \text{jet} \quad \text{(monojet)}
pp → χ χ ˉ + jet (monojet)
LHC 运行 2(13 TeV)对暗物质-核子散射截面设置了对撞机上限:
算符类型
质量范围
对撞机截面上限
对比直接探测
标量算符
m χ < 10 GeV m_\chi < 10\ \text{GeV}m χ < 10 GeV
σ S I < 10 − 40 cm 2 \sigma_{SI} < 10^{-40}\ \text{cm}^2σ S I < 1 0 − 40 cm 2
比直接探测灵敏
轴向矢算符
m χ < 100 GeV m_\chi < 100\ \text{GeV}m χ < 100 GeV
σ S D < 10 − 39 cm 2 \sigma_{SD} < 10^{-39}\ \text{cm}^2σ S D < 1 0 − 39 cm 2
补充直接探测盲区
矢量算符
m χ < 10 GeV m_\chi < 10\ \text{GeV}m χ < 10 GeV
σ S I < 10 − 41 cm 2 \sigma_{SI} < 10^{-41}\ \text{cm}^2σ S I < 1 0 − 41 cm 2
低质量区域最灵敏
对撞机对低质量暗物质(m χ ≲ 10 GeV m_\chi \lesssim 10\ \text{GeV}m χ ≲ 10 GeV )的约束优于直接探测实验,因为直接探测在低质量区间受到探测阈值的限制。
冷暗物质模拟预测星系暗物质晕的密度分布符合 Navarro-Frenk-White(NFW)轮廓:
ρ ( r ) = ρ 0 r r s ( 1 + r r s ) 2 \rho(r) = \frac{\rho_0}{\frac{r}{r_s} \left(1 + \frac{r}{r_s}\right)^2}
ρ ( r ) = r s r ( 1 + r s r ) 2 ρ 0
其中 r s r_sr s 是标度半径,ρ 0 \rho_0ρ 0 是特征密度。对于银河系,典型参数为 r s ≈ 20 kpc r_s \approx 20\ \text{kpc}r s ≈ 20 kpc ,ρ 0 ≈ 5 × 10 6 M ⊙ / kpc 3 \rho_0 \approx 5 \times 10^6\ M_\odot/\text{kpc}^3ρ 0 ≈ 5 × 1 0 6 M ⊙ / kpc 3 。
不同距离处的暗物质密度(银河系,NFW 模型):
距离 r rr (kpc)
密度 ρ ( r ) \rho(r)ρ ( r ) (GeV/cm³)
包含累积质量(10 10 M ⊙ 10^{10}\ M_\odot1 0 10 M ⊙ )
1(银心附近)
∼ 500 \sim 500∼ 500
0.1
8.5(太阳位置)
∼ 0.3 − 0.4 \sim 0.3 - 0.4∼ 0.3 − 0.4
5
20(标度半径)
∼ 0.02 \sim 0.02∼ 0.02
15
100
∼ 0.001 \sim 0.001∼ 0.001
50
200(晕边界)
∼ 2 × 10 − 5 \sim 2 \times 10^{-5}∼ 2 × 1 0 − 5
100
尖点-核点问题: 观测表明矮星系中心的暗物质密度比 NFW 预测的"尖点"(cusp,ρ ∝ r − 1 \rho \propto r^{-1}ρ ∝ r − 1 )平坦,呈"核点"(core,ρ ∼ const \rho \sim \text{const}ρ ∼ const )状分布。可能的解释包括:
超新星反馈将重子物质快速喷射出去,从而"冲淡"了暗物质尖点(在 r < 1 kpc r < 1\ \text{kpc}r < 1 kpc 处,暗物质密度降低 2-5 倍)
暗物质自相互作用(SIDM,自相互作用截面 σ / m ∼ 1 cm 2 / g \sigma/m \sim 1\ \text{cm}^2/\text{g}σ / m ∼ 1 cm 2 / g )使尖点软化
暖暗物质具有更大的自由流程(free-streaming length),天然抑制小尺度结构
数值模拟预测银河系暗物质晕内部存在数千个暗物质子晕,质量从 10 5 M ⊙ 10^5\ M_\odot1 0 5 M ⊙ 到 10 10 M ⊙ 10^{10}\ M_\odot1 0 10 M ⊙ 不等。子晕的质量函数为:
d N d M ∝ M − α , α ≈ 1.9 \frac{dN}{dM} \propto M^{-\alpha}, \quad \alpha \approx 1.9
d M d N ∝ M − α , α ≈ 1.9
这意味着质量越小的子晕数量越多。按此关系,银河系暗物质晕应该包含约:
子晕质量范围
预计数量
已观测到的发光卫星星系
> 10 10 M ⊙ > 10^{10}\ M_\odot> 1 0 10 M ⊙
∼ 1 \sim 1∼ 1
大麦哲伦云
10 9 − 10 10 M ⊙ 10^9 - 10^{10}\ M_\odot1 0 9 − 1 0 10 M ⊙
∼ 5 \sim 5∼ 5
小麦哲伦云等
10 8 − 10 9 M ⊙ 10^8 - 10^9\ M_\odot1 0 8 − 1 0 9 M ⊙
∼ 20 \sim 20∼ 20
人马座矮星系等
10 7 − 10 8 M ⊙ 10^7 - 10^8\ M_\odot1 0 7 − 1 0 8 M ⊙
∼ 100 \sim 100∼ 100
10+ 超矮星系
10 6 − 10 7 M ⊙ 10^6 - 10^7\ M_\odot1 0 6 − 1 0 7 M ⊙
∼ 500 \sim 500∼ 500
只有几个
< 10 6 M ⊙ < 10^6\ M_\odot< 1 0 6 M ⊙
∼ 10 4 \sim 10^4∼ 1 0 4
0
缺失卫星星系问题: 模拟预测的暗物质子晕远多于实际观测到的卫星星系(目前本星系群已知约 60 个矮星系)。最可能的解释是大部分低质量子晕从未吸引到足够的气体形成恒星,因此保持"黑暗"——它们确实是暗物质的天然"纯样本"。
太阳附近暗物质的局域密度是一个关键参数,直接决定了直接探测的实验灵敏度。
当前最佳测量值:
ρ 0 = 0.3 ± 0.1 GeV/cm 3 \rho_0 = 0.3 \pm 0.1\ \text{GeV/cm}^3
ρ 0 = 0.3 ± 0.1 GeV/cm 3
这个数值来自 Gaia 卫星对垂直方向恒星运动的精确测量。如果暗物质全部由 100 GeV 的 WIMP 构成,则局域数密度为:
n χ = ρ 0 m χ = 0.3 GeV/cm 3 100 GeV = 3 × 10 − 3 cm − 3 n_\chi = \frac{\rho_0}{m_\chi} = \frac{0.3\ \text{GeV/cm}^3}{100\ \text{GeV}} = 3 \times 10^{-3}\ \text{cm}^{-3}
n χ = m χ ρ 0 = 100 GeV 0.3 GeV/cm 3 = 3 × 1 0 − 3 cm − 3
即每立方厘米约 0.003 个 WIMP。WIMP 在地球附近的通量(相对速度 v ≈ 220 km/s v \approx 220\ \text{km/s}v ≈ 220 km/s ):
Φ = n χ × v × 4 π ≈ 3 × 10 − 3 × 2.2 × 10 7 × 4 π ≈ 8 × 10 5 cm − 2 s − 1 \Phi = n_\chi \times v \times 4\pi \approx 3 \times 10^{-3} \times 2.2 \times 10^7 \times 4\pi \approx 8 \times 10^5\ \text{cm}^{-2}\text{s}^{-1}
Φ = n χ × v × 4 π ≈ 3 × 1 0 − 3 × 2.2 × 1 0 7 × 4 π ≈ 8 × 1 0 5 cm − 2 s − 1
每秒每平方厘米约有 80 万个 WIMP 飞过——这就是为什么我们需要吨级探测器来捕捉那每年几次的碰撞。
尽管暗物质是占绝对主流的标准模型,少数理论试图通过修改引力定律来解释观测现象。
MOND 由 Milgrom 于 1983 年提出,其核心假设为:当加速度 a aa 小于特征值 a 0 a_0a 0 时,牛顿第二定律修改为:
F = m μ ( a a 0 ) a \mathbf{F} = m \ \mu\left(\frac{a}{a_0}\right) \mathbf{a}
F = m μ ( a 0 a ) a
其中 μ ( x ) → x \mu(x) \to xμ ( x ) → x 当 x ≪ 1 x \ll 1x ≪ 1 ,μ ( x ) → 1 \mu(x) \to 1μ ( x ) → 1 当 x ≫ 1 x \gg 1x ≫ 1 。特征加速度 a 0 ≈ 1.2 × 10 − 10 m/s 2 a_0 \approx 1.2 \times 10^{-10}\ \text{m/s}^2a 0 ≈ 1.2 × 1 0 − 10 m/s 2 (约为 10 − 10 m/s 2 10^{-10}\ \text{m/s}^21 0 − 10 m/s 2 )。
MOND 的一个核心预测——重子 Tully-Fisher 关系:
对一个旋转星系,渐近平坦速度 v f v_fv f 与星系的重子质量 M b M_bM b 满足严格关系:
v f 4 = G M b a 0 v_f^4 = G M_b a_0
v f 4 = G M b a 0
代入 a 0 = 1.2 × 10 − 10 a_0 = 1.2 \times 10^{-10}a 0 = 1.2 × 1 0 − 10 和典型星系质量 M b = 10 10 M ⊙ = 2 × 10 40 kg M_b = 10^{10}\ M_\odot = 2 \times 10^{40}\ \text{kg}M b = 1 0 10 M ⊙ = 2 × 1 0 40 kg :
v f = ( 6.67 × 10 − 11 × 2 × 10 40 × 1.2 × 10 − 10 ) 1 / 4 = ( 1.6 × 10 20 ) 1 / 4 ≈ 112 km/s v_f = (6.67 \times 10^{-11} \times 2 \times 10^{40} \times 1.2 \times 10^{-10})^{1/4}
= (1.6 \times 10^{20})^{1/4} \approx 112\ \text{km/s}
v f = ( 6.67 × 1 0 − 11 × 2 × 1 0 40 × 1.2 × 1 0 − 10 ) 1/4 = ( 1.6 × 1 0 20 ) 1/4 ≈ 112 km/s
这恰好与观测到的旋涡星系旋转曲线渐近速度范围(100-300 km/s)吻合!暗物质模型需要额外假设暗物质晕的质量-速度关系才能解释这一相关性,而在 MOND 中它是天然预测。
MOND 的弱点——无法解释的地方:
挑战
暗物质模型
MOND
子弹星系团
自然解释(分离效应)
需要额外假设
CMB 功率谱
Ω c h 2 = 0.12 \Omega_c h^2 = 0.12Ω c h 2 = 0.12 精确拟合
无法产生声学峰
大尺度结构
Λ \LambdaΛ CDM 模拟完美再现
相对论化困难
星系团质量
维里质量匹配
在团尺度需要额外暗物质
引力波速度
与光速一致
张量-标量理论限制
目前绝大多数宇宙学家认为 MOND 或其相对论化版本(如 TeVeS)不能作为暗物质的完整替代,但仍作为探索性理论被研究。
暗物质为何恰好比重子多约 5 倍?这个比值对于星系形成至关重要:
如果暗物质太多(Ω c / Ω b > 10 \Omega_c/\Omega_b > 10Ω c / Ω b > 10 ):星系及结构形成过快,气体耗散冷却效率降低,导致星系形成效率极低
如果暗物质太少(Ω c / Ω b < 1 \Omega_c/\Omega_b < 1Ω c / Ω b < 1 ):宇宙涨落的引力增长过慢,可能无法在合理时间内形成足够的星系
当前比值恰好使得星系形成既充分又有效。
标准 Λ \LambdaΛ CDM 假设暗物质是"冷"且"无碰撞"的(Cold + Collisionless)。但部分观测暗示暗物质可能存在弱自相互作用(SIDM):
自相互作用截面的观测约束:
观测方法
截面上限 σ / m \sigma/mσ / m (cm²/g)
最佳拟合值
矮星系密度核
σ / m ≲ 10 \sigma/m \lesssim 10σ / m ≲ 10
∼ 1 − 5 \sim 1 - 5∼ 1 − 5
星系团椭圆度
σ / m ≲ 0.5 \sigma/m \lesssim 0.5σ / m ≲ 0.5
≲ 0.1 \lesssim 0.1≲ 0.1
子弹星系团
σ / m ≲ 1 \sigma/m \lesssim 1σ / m ≲ 1
≲ 0.7 \lesssim 0.7≲ 0.7
哈勃望远镜强透镜
σ / m ≲ 0.1 \sigma/m \lesssim 0.1σ / m ≲ 0.1
—
不同尺度给出的约束存在张力:簇团尺度要求 σ / m ≲ 0.1 cm 2 / g \sigma/m \lesssim 0.1\ \text{cm}^2/\text{g}σ / m ≲ 0.1 cm 2 / g ,而矮星系尺度可能偏好 σ / m ∼ 1 − 5 cm 2 / g \sigma/m \sim 1-5\ \text{cm}^2/\text{g}σ / m ∼ 1 − 5 cm 2 / g 。如果暗物质自相互作用速度依赖(类似核散射),这种张力可以得到缓解。
Λ \LambdaΛ CDM 在星系以下尺度面临三个主要的"小尺度危机":
尖点-核点问题 (Cusp-Core):如前文所述,观测的矮星系中心密度比 ρ ∝ r − 1 \rho \propto r^{-1}ρ ∝ r − 1 更平坦
太-然-问题 (Too Big To Fail):模拟预测的银河系卫星星系数量远多于观测,且预测的质量比观测到的大 10-100 倍
数量-比对问题 :模拟预测的本星系群卫星星系比观测多约 10 倍
这些可能通过重子物理反馈(超新星、AGN)和暗物质自身性质(暖暗物质、自相互作用)的组合解决。2020 年代以来,更精细的模拟(如 FIRE、IllustrisTNG、EAGLE)已经大幅缩小了理论与观测的差距。
时间
进展
影响
2022
LZ 实验上线并发布首批结果
WIMP 截面上限推向新纪录
2023
PandaX-4T 发布 90 天数据
对 WIMP 的灵敏度进入 10 − 48 cm 2 10^{-48}\ \text{cm}^21 0 − 48 cm 2 量级
2024
LZ 发布 220 天数据
在 40 GeV 处 σ S I < 1.4 × 10 − 48 cm 2 \sigma_{SI} < 1.4 \times 10^{-48}\ \text{cm}^2σ S I < 1.4 × 1 0 − 48 cm 2
2024
ADMX 完成轴子高质量范围扫描
排除 2.7 − 4.9 μ eV 2.7 - 4.9\ \mu\text{eV}2.7 − 4.9 μ eV 的 DFSZ 轴子
2025
Euclid 卫星发布首批暗物质地图
弱透镜精度提高 10 倍
2026
预计的 JAXA XRISM 卫星数据
高分辨率 X 射线光谱研究星系团暗物质
暗物质是连接宇宙学、天体物理和粒子物理的核心问题:
观测证据 来自五个完全独立的方法:旋转曲线、引力透镜、子弹星系团、CMB 功率谱、大尺度结构,每个方法都一致指向暗物质的存在
占宇宙总质量的 26.8% ,是可见物质的 5.5 倍,主导了星系和宇宙结构的形成与演化
WIMP 是最优候选 ,但经过了 40 年的直接探测,仍未找到确凿信号
轴子 和惰性中微子 提供了互补的探索方向,参数空间正在被逐步覆盖
实验的灵敏度 已进入理论预期范围:LZ、XENONnT 和 PandaX-4T 正以每年 1-10 次事件的统计量搜索 WIMP
小尺度危机 和替代理论 (MOND)仍然是 Λ \LambdaΛ CDM 模型需要回应的挑战,但不构成根本否定
暗物质的发现——无论是通过直接探测、对撞机生产还是天文观测——都将开启一个超越标准模型的物理学新时代,并彻底改变我们对宇宙组成的理解。
Bertone, G., Hooper, D., & Silk, J. (2005). Particle dark matter: evidence, candidates and constraints. Physics Reports , 405(5-6), 279-390.
Planck Collaboration (2020). Planck 2018 results. VI. Cosmological parameters. Astronomy & Astrophysics , 641, A6.
Navarro, J. F., Frenk, C. S., & White, S. D. M. (1997). A universal density profile from hierarchical clustering. The Astrophysical Journal , 490(2), 493-508.
Rubin, V. C. & Ford, W. K. (1970). Rotation of the Andromeda Nebula from a spectroscopic survey of emission regions. The Astrophysical Journal , 159, 379.
Clowe, D., et al. (2006). A direct empirical proof of the existence of dark matter. The Astrophysical Journal Letters , 648(2), L109.
Akerib, D. S., et al. (2024). LZ 2024 limits on WIMP-nucleon interactions. Physical Review Letters .
Milgrom, M. (1983). A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. The Astrophysical Journal , 270, 365-370.