跳转到内容
← 返回知识库
多信使天文学约 21 分钟阅读

宇宙射线

宇宙射线GZK截断极高能宇宙射线Auger费米加速超新星遗迹

概述

宇宙射线（cosmic rays）是来自宇宙空间的高能带电粒子流，包括质子（约 90%90\%）、氦核（约 9%9\%）、重核（约 1%1\%）和少量电子、反质子。其能量范围跨越约 $14$ 个数量级——从 10910^9 eV（GeV）到最高记录的约 3×10203 \times 10^{20} eV（$300$ EeV，"哦，上帝"粒子——Oh-My-God Particle，犹他大学 Fly's Eye 探测器 1991 年记录）。

宇宙射线打在大气层上产生粒子级联（空气簇射），是 20 世纪前半叶高能物理的主要实验场所：正电子、μ\mu 子、π\pi 介子都首先在宇宙射线研究中发现。然而，宇宙射线最高能量端的来源，至今仍是天文学和粒子物理最大的未解谜题之一。

名字里的“射线”是历史误会：它们是带电粒子，会被磁场反复偏折，到达方向几乎不指向源。大气层本身就是量能器，原初粒子在高空打出空气簇射，地面阵列和荧光望远镜量的是次级粒子和氮气荧光，而不是把那个质子抓在手里。正电子、μ 子和 π 介子之所以先在宇宙射线里被发现，是因为簇射把一台天然加速器的产物送到了海平面，实验室加速器当时还够不到同样的质心能量。

能谱：从 GeV 到 ZeV

宇宙射线能谱约遵循幂律分布 ϕ(E)∝E−γ\phi(E) \propto E^{-\gamma}，但在几个关键能量处发生明显变化：

膝（Knee，$\sim 10^{15}$ eV = 1 PeV）

谱指数从 γ≈2.7\gamma \approx 2.7 变为 γ≈3.0\gamma \approx 3.0。膝的位置对应银河系加速机制的上限：

目前主流认为：超新星遗迹（SNR）中的费米 I 型扩散激波加速（DSA）能将银河系宇宙射线加速到约 PeV 量级（"PeVatron"），而膝以上银河系加速机制基本达到上限，需要河外或更高效的加速机制。

LHAASO 在 2021 年首次明确证认了 PeVatron 的存在：在蟹状星云等 12 个超高能 γ\gamma 射线源中，发现能量超过 $100$ TeV 的光子，间接证明了其中有能量超过 PeV 的宇宙射线加速粒子（LHAASO Collaboration 2021，Nature，594，33）。

膝附近的转折不是能谱自己无缘无故折了一下，而是银河系加速器开始交不出更高能的那些轻核。费米 I 型扩散激波加速把超新星遗迹里的质子推到约 PeV，再往上要么换更重的核、要么换更大的加速区。LHAASO 在蟹状星云等源里看到超过 $100$ TeV 的光子，间接证明有粒子已经被推到 PeV 量级，但它证明的是“存在 PeVatron”，不是“膝以上全部仍由遗迹负责”；成分若在膝附近逐渐变重，能谱变陡也可以来自轻核先到达上限，重核稍后才掉下去。

踝（Ankle，$\sim 3 \times 10^{18}$ eV = 3 EeV）

谱指数回硬到约 γ≈2.6\gamma \approx 2.6。踝被广泛认为是河外宇宙射线（来自银河系以外）开始主导的标志。踝以上的粒子回旋半径超过银河系尺度，无法被银河系磁场约束（回旋半径 rL=E/(ZeB)>1 kpcr_L = E/(ZeB) > 1\,\text{kpc}），必须来自河外。

踝附近谱变硬，通常读成河外部件开始压过银河系那一段贡献，因为回旋半径 rL=E/(ZeB)r_L = E/(ZeB) 在文中给出的约 3×10183 \times 10^{18} eV 处已经大于千秒差距量级，银河磁场关不住这些粒子，它们必须来自河外。关不住不等于加速失败：银河系仍可能贡献到踝附近，只是不再能提供各向同性的约束，因此踝是约束失效的标志，也是成分和各向异性分析必须换一套参考系的地方。

GZK截断（$\sim 5 \times 10^{19}$ eV = 50 EeV）

1966 年，格里森（Greisen）和扎采宾-库兹明（Zatsepin-Kuzmin）独立预言：能量超过约 5×10195 \times 10^{19} eV 的质子在穿越宇宙空间时，会与宇宙微波背景（CMB）光子发生光核反应，产生 Δ+\Delta^+ 共振衰变：

p+γCMB→Δ+→p+π0或n+π+p + \gamma_{CMB} \to \Delta^+ \to p + \pi^0 \quad \text{或} \quad n + \pi^+

这一过程的能量损失自由程约 $50$–$100$ Mpc，意味着我们能观测到的能量超过 GZK 截断的宇宙射线，其来源必须在约 $100$ Mpc 以内。

2007–2008 年，HiRes 实验（Abbasi et al. 2008，Physical Review Letters，100，101101）和 Pierre Auger Observatory（Abraham et al. 2008，Physical Review Letters，101，061101）均独立探测到 GZK 截断（或"压制"）的统计证据，确认了极高能宇宙射线必须来自宇宙学邻域。

GZK 过程用的是无处不在的宇宙微波背景光子当靶，阈值由 Δ+\Delta^+ 共振决定，文中把截断写在约 5×10195 \times 10^{19} eV；自由程只有约 $50$–$100$ Mpc，所以截断以上的粒子是本地样品，不是把全宇宙的加速器积分进来。HiRes 与 Auger 在 2007–2008 年各自看到流量压制，确认的是“有压制”，至于压制完全来自质子光介子反应，还是重核光分解也出了一份力，要看最高能端的成分；成分若偏向重核，纯质子 GZK 图像就要改写，但截断本身仍然是一条距离约束。

加速机制：费米加速

主流理论认为宇宙射线通过费米加速（Fermi acceleration）获得能量：

I 型（扩散激波加速，Diffusive Shock Acceleration，DSA）：粒子在激波面两侧多次穿越，每次穿越获得一次能量提升 ΔE/E≈4vs/(3c)\Delta E/E \approx 4v_s/(3c)（vsv_s 是激波速度）。大量穿越统计上产生幂律谱 N(E)∝E−γN(E) \propto E^{-\gamma}，γ=(r+2)/(r−1)\gamma = (r+2)/(r-1)（$r$ 是激波压缩比，强激波 $r = 4$，给出 γ=2\gamma = 2，与观测基本吻合）。

超新星遗迹（SNR）被认为是银河系宇宙射线的主要加速场所，理由：

  1. SNR 激波速度 vs∼104v_s \sim 10^4 km/s，足以加速质子到 PeV
  2. 银河系 SNR 每 $30$–$50$ 年发生一次，能量输出与观测到的宇宙射线能量密度一致（前提是约 10%10\% 的超新星爆发能量转化为宇宙射线）
  3. HESS、MAGIC、VERITAS 等伽马射线望远镜在多个 SNR 中探测到 TeV 伽马射线，间接证明了高能粒子加速

II 型（随机费米加速）：粒子随机地与磁云碰撞，平均而言也能获得能量，但效率低，加速时间长，对极高能端可能有补充作用。

费米 1949 年的原始图像是粒子与随机磁云碰撞，二阶过程，效率按 (v/c)2(v/c)^2 计，太慢，很难单独支撑膝以上的那一段。激波把相对速度变成有方向的一阶过程，每次穿越的分数增益约为 4vs/(3c)4v_s/(3c)，强激波压缩比 $r=4$ 给出 N(E)∝E−2N(E)\propto E^{-2}，这才和观测到的注入谱大致同阶（Fermi 1949，Physical Review，75，1169）。观测到的 γ≈2.7\gamma \approx 2.7 比注入谱更陡，是因为粒子在银河系里扩散逃逸时还要再丢一截能量，膝附近的变陡则叠加了加速上限。

Hillas 后来用“尺度 × 磁场”画了一张图，说明要把粒子推到踝以上，加速区必须足够大、磁场足够强，超新星遗迹很难单独占据那一块。图是必要条件，不是源的身份证：满足 Hillas 条件的候选仍有活动星系核、伽马射线暴和星暴星系风，观测还没有把票投给其中某一个（Hillas 1984，ARA&A，22，425）。要把质子留在激波附近足够久，还需要激波本身放大磁场；没有放大，遗迹的名义尺度和平均磁场往往不够把粒子送到膝。能量预算与谱形还要同时过关：注入谱太硬会过度生产次级伽马与中微子，太软又填不满观测到的银河系宇宙射线能量密度——多信使上限因此反过来筛选加速效率。

极高能宇宙射线与来源

Pierre Auger Observatory（PAO）

皮埃尔·奥热天文台位于阿根廷门多萨省，由 $1600$ 个水切伦科夫探测器阵列（3000 km23000\,\text{km}^2）和 $24$ 台荧光望远镜（分布在 $4$ 个荧光站，每站 $6$ 台）组成，是世界最大的宇宙射线探测器（Auger Collaboration 2008；建造：1999-2008 年）。

主要发现：

  • 能量超过 $55$ EeV 的事件数目（约 $70$ 个事件，截至 2023 年）
  • 各向异性：极高能宇宙射线（$> 8$ EeV）的到达方向有微弱但统计显著的偏向活动星系核方向（可能是活动星系核的角关联，Auger Collaboration 2017，Science，357，1266），超出各向同性分布置信度约 5.2σ5.2\sigma，但具体来源仍有争议
  • 成分：极高能端（$> 50$ EeV）的宇宙射线成分偏向重核（铁、氧），而非纯质子。这与某些"纯质子"GZK 模型不符，使极高能宇宙射线的来源问题更加复杂

Auger 的关键不是水罐数目，而是地面阵列与荧光望远镜组成的混合探测：地面阵列几乎全天候计粒子，荧光望远镜在晴夜看簇射的纵向发展，后者对 Xmax⁡X_{\max} 敏感，这才谈得上最高能端偏重核还是偏质子。重核在大气里更早发展，平均 Xmax⁡X_{\max} 更浅，质子则穿得更深；强子相互作用模型会平移这条尺子，所以成分结论必须连同模型一起读。各向异性在 $>8$ EeV 相对活动星系核方向超出约 5.2σ5.2\sigma，统计上不像各向同性，但“角关联”不等于已经指认了单个源。能量标在 $55$ EeV 以上的事件仍然稀少，上文写到约 $70$ 个（截至 2023 年），任何新的源模型都先得过这套小样本的统计关。

Telescope Array（TA）

望远镜阵列位于美国犹他州，是北半球最大的极高能宇宙射线实验，与 PAO 互相补充（覆盖北半球天空）。TA 的成分分析倾向于质子或轻核主导，与 PAO 的结果存在一定分歧，可能反映了南北半球观测方向上真实的各向异性，也可能是系统误差。

"热点"：TA 在大熊座方向发现约 $20°$ 的过密度（Abbasi et al. 2014，The Astrophysical Journal Letters，790，L21），超过各向同性预期约 3.4σ3.4\sigma，可能对应一个局域加速源（如∼100\sim 100 Mpc 内的活动星系核），但尚无定论。

望远镜阵列覆盖北天，与南半球那座皮埃尔·奥热天文台并不看同一片天空，成分分歧因此有两种完全不同的读法。一种是南北天的源分布或传播路径本来不同，热点所在的大熊座方向过密度约 3.4σ3.4\sigma，可能对应约 $100$ Mpc 内的局域加速器；另一种是两套实验的能量标尺或强子模型系统误差还没对齐。在重叠天区交换分析方法之前，不宜把某一方的成分结论写成极高能宇宙射线的全球性质，那会把仪器差误读成天空差。

宇宙射线与多信使天文学

宇宙射线本身不是"信使"（它们是带电粒子，会被磁场偏转，不能直接指向来源），但通过多信使关联可以约束来源：

中微子：高能宇宙射线与物质（pp 碰撞）或光子（pγ\gamma 碰撞）相互作用产生 π\pi 介子，后者衰变产生高能中微子（和伽马射线）。IceCube 中微子天文台已探测到宇宙学来源的弥散高能中微子（TeV–PeV 量级，IceCube Collaboration 2013，Science，342，1242），并在 NGC 1068（活动星系核）和银河系平面方向发现了中微子超出信号，为宇宙射线加速提供了直接证据。

伽马射线：中性 π0\pi^0 衰变产生的 γ\gamma 射线光子是中性粒子，不受磁场偏转，可以直接追溯来源。Fermi-LAT 和地面大气切伦科夫望远镜（HESS、MAGIC、VERITAS、LHAASO）探测到的 TeV–PeV 伽马射线源是寻找宇宙射线加速场所的关键。

带电宇宙射线不能当指向信使，中微子和伽马射线却可以沿直线回溯到来向：高能质子在源区或传播中打出 π 介子，带电 π 衰变成中微子，中性 π 衰变成伽马射线，三条链的能量账必须能对上。IceCube 的弥散 TeV–PeV 中微子证明河外有强子加速，NGC 1068 方向的超出把活动星系核核区推进候选名单，但单个源的中微子计数仍然少，电磁对应体可以是核的 X 射线，也可以是被挡住的喷流。中微子探测器和地面切伦科夫望远镜因此不是宇宙射线实验的附件，而是给“源在哪”提供直线轨迹的那一半仪器。

2017 年 IceCube-170922A 与耀变体 TXS 0506+056 的伽马射线活跃期在时间上重合，是把单个高能中微子和具体电磁源放到一起的著名例子，但仍是小统计量上的成协，不能单独把极高能宇宙射线的源问题结案（IceCube Collaboration 2018，Science，361，eaat1378）。多信使天文学在这里要求的是同一方向、同一时间窗口里的电磁 SED，而不是事后把最亮的伽马源硬配给每一个中微子。真正有判决力的是：同一候选源是否能同时解释宇宙射线能谱、中微子流量与伽马射线不透明度；缺一条链，就仍只是“方向碰巧亮了一下”。

为什么这很重要

宇宙射线是从宇宙深处传来的高能物质样本，携带着关于极端天体物理过程（超新星、AGN、GRB）的直接信息。它们的能量远超人类建造的任何粒子加速器（LHC 最高能量约 $14$ TeV，而 GZK 粒子能量约 $300$ EeV），是研究超高能粒子物理的唯一窗口。解开极高能宇宙射线来源之谜，将需要宇宙射线、中微子和伽马射线三种信使的联合分析——这正是多信使天文学的精髓。

地面上的大型强子对撞机把质子加速到约 $14$ TeV，文中记录的最高端约 $300$ EeV，中间这十几个数量级目前只能向天空去借。膝、踝与 GZK 截断是借来的三枚路标：膝问银河加速器的上限，踝问约束何时失效，GZK 问本地宇宙的光子靶。Auger 与 TA 把路标变成可以争论的谱和成分，下一步不是再加一个能量数字，而是让中微子和伽马射线把候选源从统计关联变成可重复的多信使事件。粒子加速器在地面给出的是可以调旋钮的对照实验，天空给出的是无法改参数的同一类高能碰撞。成分与各向异性若长期南北不一致，还会把“全球源类”叙事拆成局域天区问题——那本身也是科学进展，而不是实验失败。

跨域连接

  • 超新星遗迹：膝的位置对应银河系加速器的上限：遗迹激波速度足以把质子推到千万亿电子伏量级，再往上银河系机制基本到顶。推论是：膝不是能谱自身的性质，而是加速器的性质；在若干源里探到能量极高的伽马光子，就间接证明了确实存在这一量级的加速。
  • 随机过程：粒子在激波两侧反复穿越，每次得到一个固定比例的能量增益，同时有固定概率逃逸——"按比例增长加随机停止"这一组合必然给出幂律。推论是：谱指数只由激波压缩比决定，强激波给出的指数与观测基本吻合，因此观测到的指数量的是激波强度，不是源的种类。
  • 磁场：带电粒子的回旋半径正比于能量、反比于电荷与磁场强度，踝以上时它已超过千秒差距，银河系磁场再也困不住。推论是：踝的物理含义是约束失效而不是加速失效，这就把踝以上的来源逼到河外；同一条公式也说明为什么这类粒子无法直接指向源。
  • 粒子加速器：这里的能量损失机制就是固定靶实验——极高能质子撞上无处不在的微波背景光子，质心系能量一旦越过共振阈值就开始丢能量。推论是：能量损失自由程只有几十到上百兆秒差距，于是"我们收到了截断以上的粒子"本身就是一条距离约束：源必须在这个半径以内。
  • 可重复性危机：南北两半球最大的两个实验在成分上给出不一致的结果，一边偏重核，一边偏质子或轻核。推论是：出路只有两条——真实的南北天各向异性，或某一台的系统误差；判据不是把统计误差做得更小，而是两台交换分析方法、在重叠天区上直接比对。

参考文献

  • LHAASO Collaboration (2021). Identification of twelve sources as PeVatrons in the Milky Way. Nature, 594, 33–36.
  • Abraham, J. et al. (Pierre Auger Collaboration) (2008). Observation of the suppression of the flux of cosmic rays above 4×10194 \times 10^{19} eV. Physical Review Letters, 101, 061101.
  • Auger Collaboration (2017). Observation of a large-scale anisotropy in the arrival directions of cosmic rays above 8×10188 \times 10^{18} eV. Science, 357, 1266–1270.
  • IceCube Collaboration (2013). Evidence for high-energy extraterrestrial neutrinos at the IceCube detector. Science, 342, 1242856.
  • Greisen, K. (1966). End to the cosmic-ray spectrum? Physical Review Letters, 16, 748–750.
  • Fermi, E. (1949). On the Origin of the Cosmic Radiation. Physical Review, 75, 1169–1174.
  • Hillas, A. M. (1984). The Origin of Ultra-High-Energy Cosmic Rays. Annual Review of Astronomy and Astrophysics, 22, 425–444.
  • Abbasi, R. U. et al. (Telescope Array Collaboration) (2014). Indications of Intermediate-Scale Anisotropy of Cosmic Rays with Energy Greater Than 57 EeV in the Northern Sky Measured with the Surface Detector of the Telescope Array Experiment. The Astrophysical Journal Letters, 790, L21.
  • IceCube Collaboration et al. (2018). Multimessenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A. Science, 361, eaat1378.

延伸阅读

  • Blasi, P. (2013). The origin of galactic cosmic rays. The Astronomy and Astrophysics Review, 21, 70. — 银河系宇宙射线起源的理论综述
  • Alves Batista, R. et al. (2019). Open questions in cosmic-ray research at ultrahigh energies. Frontiers in Astronomy and Space Sciences, 6, 23. — 极高能宇宙射线开放问题综述
  • Kampert, K.H. & Unger, M. (2012). Measurements of the cosmic ray composition with air shower experiments. Astroparticle Physics, 35, 660–678.