6.1 宇宙线起源之争

6.1 宇宙线起源之争

本节摘要：GeV 到膝盖段（约 10^15 eV）的宇宙线公认主要来自银河系内超新星遗迹，证据链相对完整；争论集中在膝以上——最高能段是河外源吗？成分为什么变重？Auger 的成分证据让"质子一统天下"的简单图像让位于混合图像，但源的指认仍未完成。

一、已确立的部分

• 能谱整体幂律 + 膝 + 踝 + GZK 压低，形态稳固（第 2.1 节）。
• 超新星遗迹中观测到 π⁰ 凸包（W44、IC 443），银河系内质子加速确证。
• 能量记账：超新星率与银河系宇宙线亮度相容。

二、争论焦点：膝以上归谁

膝盖能量以上，候选解释分三派：

Auger 的 Xmax 测量（簇射在大气中的最大深度）显示：能量超过踝后簇射变浅、成分趋重——支持混合图像，即最高能段更可能是铁核而非质子。但 Xmax 的解释依赖强相互作用模型（补救措施：用 LHC 前向物理对撞数据校准，如 LHCf 实验），模型不确定性是当前最大系统误差。

三、指认源：为什么这么难

回顾第 4.1 节：磁场把方向搅散。即把偏转角压到几度，还需要知道路径上的磁场分布——而星系际磁场本身就是未知量（第 4.3 节的级联约束只给出量级）。因此超高能宇宙线的"热点"（Auger 与 TA 各自看到微弱的到达方向成团）只能算线索，不是指认。

def deflection_degree(Z, E_EeV, B_nG, D_Mpc):
    """粗糙估算：磁偏转角（度）。Z电荷数、E能量EeV、B星系际磁场nG、D传播距离Mpc"""
    return 16.0 * Z / E_EeV * (B_nG / 1.0)**0.5 * (D_Mpc / 10.0)**0.5 * 2.0

print(f"60 EeV 质子穿 50 Mpc, 2 nG: 偏转约 {deflection_degree(1, 60, 2, 50):.0f} 度")
print(f"60 EeV 铁核同条件: 偏转约 {deflection_degree(26, 60, 2, 50):.0f} 度")
​

这段粗算说明一切：60 EeV 的质子偏转数度尚可容忍，铁核偏转近百度——成分变重直接摧毁方向指认能力。这是 Auger 成分结果让人沮丧的根本原因，也是中微子（无偏转）在指认问题上被寄予厚望的原因。

图：宇宙线起源的证据版图

本节要点回顾

• 低能段已结案：SNR 加速银河系宇宙线，π⁰ 证据与能量记账双确认。
• 膝以上是主战场：三派解释竞争，Auger 的 Xmax 数据偏向混合图像。
• 强子模型不确定性是成分分析的最大系统误差，靠 LHC 前向测量校准。
• 偏转粗算一秒钟：成分×磁场一变重，方向指认即失效——中微子才是终极证人。

动手估算：膝盖以下的能量预算对不上账

银河系宇宙线的能量密度约 1 eV/cm³，盘体积约 πR²H ≈ π(15kpc)²(0.2kpc) ≈ 2.7×10^67 cm³，总能量约 4×10^66 eV ≈ 6×10^55 erg。超新星率约 3/世纪、每次可用 10^50 erg，注入功率 10^41 erg/s，维持能量密度需要驻留约 10^7 年——账是平的，但每克星际物质需要的注能效率约 10% 恰好位于激波加速理论上限附近：起源之争因此不只是"找源"，而是效率之争。

# 能量预算收支表
import math
u = 1.0                          # eV/cm^3
V = math.pi * (15*3.086e21)**2 * (0.2*3.086e21)
E_total = u * V * 1.602e-12      # erg
sn_rate, e_sn = 3.0/3.15e7/100, 1e50
tau = E_total / (sn_rate * e_sn)
print(f"总能量 {E_total:.1e} erg, 需驻留 {tau/3.15e7/1e6:.0f} Myr")
# 输出约 6 Myr：与 Be-10 时钟一致，但注能效率须接近全部理论的极限
​

争论焦点数据卡

# 用各向异性幅度判断"膝盖以上是否仍河内"
# 河内稳态扩散预言 delta ~ (6D/cH)*(E/ZeV)^delta 而河外各向同性
import math
D = 3e28 * (1e16/3e9)**0.35      # PeV 的扩散系数
H = 4*3.086e21
I = 3*D/(2.998e10*H)
print(f"扩散预言 PeV 各向异性 ~ {I:.1e}（实测 ~1e-3 上升）")
# 实测略高:近源瞬时贡献或轻组分剩余，仍是开放争论
​

延伸讨论：百年争论的裁决机制

起源之争从 1912 延续至今而不崩塌，是因为每个能量段都缺"决定性证据"：膝盖以下缺成分的 Z 分辨（需空间量能器如 HERD），膝以上缺伽马/中微子伴随（LHAASO 与 IceCube 各自在攒统计），踝以上缺到达方向的浅层各向异性（Auger 增强在补 Xmax 测量）。可以预见的三种裁决剧本：LHAASO 在 SNR 壳上看到 π0 特征硬谱；IceCube 在银河面看到弥散中微子；或都不发生，把答案推向更奇特的机制（磁星、原初黑洞）。无论哪种，观测窗口都已就位。

争论还牵动一个常被忽略的工程判断：若膝盖以上仍是河内，下一代广延阵列应加厚 μ子层以分辨成分；若踝以上转为河外，则应加密荧光眼以测 Xmax。两条技术路线之争实质是"起源假设的实验投票"。

最后，这场争论的教学价值超越其本身：它示范了一门成熟学科如何在"每个环节都有数据、整体结论仍未定"的状态下有序推进——每一代实验把假设空间削掉一角，而不是等待一次性的决定性判决。