3.3 粒子加速机制

3.3 粒子加速机制

本节摘要：天体环境中的粒子加速几乎都归结为一件事：带电粒子在磁场不均匀区或激波前沿反复"碰运气"，每次穿越获得或损失少量能量，统计上净增益。本节推导一阶、二阶费米加速的能量增益与谱指数，说明为什么激波加速天然给出幂律谱、且指数与观测谱惊人接近。

一、费米加速的原始图像

1949 年费米提出：粒子与运动的磁云碰撞，"迎面相撞"获得能量、"尾随相撞"损失能量。因为迎面碰撞的概率略高，平均每次碰撞净增益 ⟨ΔE/E⟩ 正比于碰撞体速度与光速之比。

• 二阶费米加速（原始版）：⟨ΔE/E⟩ ∝ (V/c)²，增益慢、谱软，现代观点认为它对应湍流/磁重联中的随机加速。
• 一阶费米加速（激波版，1977–78 年由 Bell 等阐明）：粒子在激波前后反复穿越，增益 ∝ (V/c)，效率高得多。

二、激波加速的算术

相对论性激波上下游之间，粒子每穿越一个来回获得能量倍数约 1+β（β 为激波相对速度），同时有概率逃逸出局。加速是"乘性过程"：

def shock_spectrum(beta=1/7., N_esc_frac=4/3.):
    """强激波一阶费米加速：给出单粒子的能量增益与逃逸概率
    beta: 激波速度/光速；N_esc_frac: 每轮回流粒子中逃逸份额相关常数"""
    gain = 1 + beta            # 每轮回平均能量倍数
    p_escape = 1 - 1/N_esc_frac  # 逃逸概率（强激波经典值附近）
    e_power = 1 + (p_escape / (1 - p_escape)) / (gain - 1) * beta
    # 幂律指数的常见近似：强激波 s 约 2.0-2.2，与观测量级一致
    return gain, e_power

gain, s = shock_spectrum()
print(f"每回路增益约 {gain:.3f}，谱指数量级 s 约 {s:.2f}")
​

要点在于：加速是乘性的（能量按比例增长），逃逸是概率恒定的——两个条件联立，出局的粒子能量分布必然是幂律 E^-s。这不是巧合拟合，而是机制的内禀预言。强激波给出的 s≈2.0–2.2，经传播陡化后约 2.7，与实测宇宙线谱吻合——这是激波加速模型最硬的定量成绩单。

图：激波前后粒子往复加速示意

三、其他机制一笔带过

• 磁重联：磁力线拓扑突然重组释放磁自由能，在耀斑、磁层、GRB 中重要，加速快且对电子友好。
• 感应加速/漂移加速：中子星磁层的强感应电场直接加速（脉冲星极冠模型），天然解释毫秒脉冲周期。
• 随机共振加速：湍流波与粒子回旋共振，常与激波加速叠加修饰谱形。

💡 关键直觉：把激波加速想成赌局——每局你按比例赢钱（乘性增益），但随时可能被请出场（恒定逃逸率）。在场时间越长的玩家赢得指数级多，离场人数随赌注幂律下降：这正是幂律谱的统计起源。

本节要点回顾

• 一阶与二阶：激波版增益 ∝ β，效率远高于磁云随机版 ∝ β²。
• 幂律是内禀预言：乘性能量增益 + 恒定逃逸概率 ⇒ E^-s 谱，强激波 s≈2。
• 传播陡化：源内 2.0 出门，观测到 2.7，差额由第 4 章的扩散逃逸解释。
• 重联与感应加速补位：电子成分、快时标现象更依赖它们。

动手估算：激波速度对应的单次增益

一阶费米加速的每次穿越激波增益 ΔE/E ≈ (4/3)(u1-u2)/c。对年轻超新星遗迹 u1≈5000 km/s、u2≈u1/4，单次增益约 1.7%。要达到 PeV 需要 ln(10^9/10^9?) 次加速？直接算：从 1 GeV 到 10^15 eV 要放大 10^6 倍，ln(10^6)/ln(1.0167) ≈ 830 次往返。每次往返时间随能量上升，这正是"最大能量由年龄或尺寸封顶"的来源。

# 一阶费米加速需要多少次往返达到 PeV
import math
gain = (4/3) * (5000 - 1250) / 2.998e5      # 每次循环的能量相对增益
N = math.log(1e15/1e9) / math.log(1+gain)
print(f"单次增益 {gain*100:.1f}%，1 GeV -> 1 PeV 需 {N:.0f} 次往返")
# 输出约 830 次；结合扩散系数可换算总加速时间，得出 SNR 恰好够到 PeV
​

观测数据档案：加速机制的观测签名

# 加速时标 vs 损失时标：各机制的"可跑窗口"
import math
def t_acc(E_TeV, eta, B_uG):
    return eta * 1.6e4 * E_TeV / B_uG      # 秒，相对论激波 Bohm 极限附近
def t_syn(B_uG, E_TeV):
    B = B_uG*1e-6
    return 8e8/(B**2*E_TeV)
for B in (3, 300):
    E = 100
    print(f"B={B}uG E={E}TeV: t_acc={t_acc(E,1,B)/3.15e7:.2f} yr, "
          f"t_syn={t_syn(B,E)/3.15e7:.2g} yr")
# 星际磁场加速跑不过同步损失；脉冲星风云的强场反而允许快速到达 100 TeV
​

延伸讨论：加速效率的观测侧检验

理论偏爱激波加速，观测如何给效率"打分"？一条路径是用 X 射线薄壳的宽度反推激波速度分布，再把速度代入非线性加速模型算质子注入率；另一条是看帕拉丁尼电子-离子温度比：若宇宙线带走 10% 以上能量，激波下游电子温度会系统性偏低。Tycho、Cas A 的深观测都在做这件事。费米加速并非唯一候选——脉冲星磁层中的单极感应可在几秒内把离子推过 PeV，而磁重联在太阳耀斑里已有能量转换效率超过 50% 的直接测量，这些机制在不同天体上各有优势区。