First detection of ultra-high energy emission from gamma-ray binary LS I +61 303

LHAASO 合作组利用其数据首次探测到伽马射线双星系统 LS I +61 303 的超高能（>100 TeV）辐射，并观测到 3.9σ 置信度的轨道调制特征，为极端粒子加速及轻子与强子过程共同作用的复合模型提供了有力证据。

要点

这是一篇关于天文学重大发现的科普解读。简单来说，中国科学家利用“拉索”（LHAASO）望远镜，第一次在宇宙中一个特殊的“双星系统”里，捕捉到了能量极高、极其罕见的伽马射线。

为了让你轻松理解，我们可以把这篇论文的内容想象成一场**“宇宙超级马拉松”的现场报道**。

1. 主角是谁？（LS I +61°303）

想象宇宙中有一对**“怪咖搭档”**：

• 大哥：一颗巨大的、年轻的恒星（Be 型星），它像是一个正在疯狂喷水的消防栓，不断向外喷射强烈的恒星风和辐射。
• 小弟：一个看不见的“幽灵”，科学家推测它可能是一颗快速旋转的中子星（脉冲星）。它像是一个在消防栓旁边高速奔跑的**“超级短跑运动员”**。

这对搭档手拉手（互相绕转），每 26.5 天转一圈。当“小弟”跑到大哥喷出的“水雾”（恒星风）最浓密的地方时，就会发生剧烈的碰撞，产生高能粒子。

2. 发现了什么？（超高能伽马射线）

过去，我们只能看到这对搭档在“中等能量”下的表现（比如几万亿电子伏特，TeV）。但这次，“拉索”（LHAASO）望远镜——一个建在四川稻城海拔 4400 米高山上的超级“捕网”——捕捉到了**“超高能”（UHE）**的信号。

• 比喻：以前我们只能看到这对搭档扔出的“网球”（普通高能粒子），这次我们竟然看到了他们扔出的**“加特林机枪子弹”**（能量超过 100 万亿电子伏特，即 100 TeV 以上）。
• 意义：这证明了宇宙中存在一种极其强大的“粒子加速器”，能把粒子加速到人类无法想象的极限速度。这就像是在一个小小的太阳系里，发现了一个能制造“宇宙级核弹”能量的工厂。

3. 怎么发现的？（拉索的“三眼”）

“拉索”望远镜由三部分组成，就像有三个不同功能的眼睛：

• WCDA（水切伦科夫探测器阵列）：像一个巨大的**“浅水池”**，专门捕捉中等能量的“网球”。它发现信号非常强（9.2 倍的标准差，相当于在人群中一眼认出了那个穿红衣服的人）。
• KM2A（平方公里阵列）：像一个巨大的**“地面捕网”**，专门捕捉极高能量的“子弹”。它虽然只抓到了 16 个超高能事件（背景噪音只有 5 个），但这足以证明它们确实存在（3.8 倍标准差）。
• 最厉害的一点：他们发现，能量最高的那些“子弹”，并不是随时都有的，而是有节奏地出现。

4. 有趣的节奏（轨道调制）

这对“搭档”在绕圈跑的时候，发射高能粒子的规律很有趣：

• 低能量时：在轨道的某个位置（远日点附近）发射得比较猛。
• 超高能量时：科学家发现，那些能量最高的“子弹”，似乎更喜欢在轨道的另一侧（近日点附近，也就是离大哥最近、碰撞最激烈的时候）出现。

比喻：想象那个“短跑运动员”在绕圈跑。当他跑到大哥喷水的“死角”或者“最浓密的水雾”里时，他不仅跑得更快，还能把水花甩得更高、更远。这种**“能量越高，出现时间越特殊”**的现象，是解开宇宙谜题的关键线索。

5. 他们是怎么做到的？（两种“引擎”）

科学家在讨论：到底是什么机制产生了这么高的能量？

• 纯电子引擎（轻子模型）：就像用电子去撞击光子。但在超高能下，这就像**“用乒乓球去撞保龄球”**，效率会急剧下降（康普顿散射的克莱因 - 仁科效应限制），很难解释为什么会有那么高的能量。
• 质子引擎（强子模型）：就像用**“质子（原子核）”**去撞击。质子更重、更硬，能扛住更高的能量。
• 结论：科学家认为，这很可能是一个**“混合双打”**。
  • 低能量时，主要是电子在起作用（像轻快的乒乓球）。
  • 超高能量时，主要是质子在起作用（像沉重的铁球）。
  • 这对搭档不仅是个“加速器”，更是一个**“银河系内的 PeV 加速器”（PeVatron）**。这意味着它有能力把质子加速到 100 万亿电子伏特（PeV）级别，这通常是只有超新星遗迹或黑洞喷流才有的能力，而这里只是一个双星系统！

6. 总结：为什么这很重要？

这篇论文告诉我们：

1. 宇宙比我们想象的更狂野：在一个只有两颗恒星的小系统里，竟然藏着能制造“宇宙级子弹”的工厂。
2. 新工具的力量：中国的“拉索”望远镜就像一把新钥匙，打开了通往“超高能宇宙”的大门，让我们看到了以前看不到的景象。
3. 物理学的突破：这挑战了我们对粒子加速机制的旧认知，告诉我们宇宙中的粒子加速可能比我们想象的更复杂、更高效。

一句话总结：
中国科学家利用“拉索”望远镜，在宇宙中一对“跳舞”的恒星系统里，第一次抓到了能量极高的“宇宙子弹”，证明那里藏着一个能把粒子加速到极限的超级工厂，而且这个工厂的运作规律非常有趣，可能是电子和质子“混合双打”的结果。

技术摘要

以下是基于论文《First detection of ultra-high energy emission from gamma-ray binary LS I +61◦303》（首次探测到伽马射线双星 LS I +61◦303 的超高能辐射）的详细技术总结：

1. 研究背景与科学问题 (Problem)

• 研究对象：LS I +61◦303 是典型的伽马射线双星系统，由一颗大质量 Be 型恒星和一个致密天体（近期观测表明可能是一颗旋转的中子星/脉冲星）组成。该系统具有约 26.5 天的轨道周期，并在多个波段表现出轨道调制。
• 科学挑战：
  • 尽管该系统已被 MAGIC 和 VERITAS 等望远镜在 TeV 能段（<100 TeV）探测到，但其粒子加速机制的上限尚不清楚。
  • 现有的轻子模型（逆康普顿散射）在超高能（UHE, >100 TeV）区域面临 Klein-Nishina 截面的抑制和同步辐射冷却的限制，难以解释如此高能光子的产生。
  • 需要确认此类双星系统是否具备将质子加速至 PeV（$10^{15}$ eV）能量的能力（即是否为银河系内的“PeVatron"候选体）。
• 核心问题：LS I +61◦303 能否发射能量超过 100 TeV 的超高能伽马射线？其辐射机制是纯轻子、纯强子还是混合模型？

2. 研究方法 (Methodology)

• 观测设备：利用高海拔宇宙线观测站 (LHAASO) 的数据。LHAASO 包含三个主要阵列，本研究主要使用了：
  • WCDA (水切伦科夫探测器阵列)：覆盖 1.4–30.5 TeV 能段。
  • KM2A (平方公里阵列)：覆盖 25 TeV 至数百 TeV 能段，特别用于探测 >100 TeV 的超高能事件。
• 数据分析流程：
  • 数据选取：WCDA 数据（2021 年 3 月 -2024 年 7 月）和 KM2A 数据（2019 年 12 月 -2024 年 7 月）。
  • 背景抑制：使用“直接积分法”估计宇宙线背景，并利用 $\gamma$/强子分离参数（PINCness）进行筛选。
  • 统计方法：采用 3D 似然拟合（3D Likelihood Fitting），同时拟合源的位置、光谱指数和归一化。考虑了银河系弥散辐射（GDE）和邻近源 1LHAASO J0249+6022 的影响。
  • 轨道调制分析：将事件按轨道相位（$\phi$）分箱（每轨道 10 个相位箱），利用似然比检验（Likelihood-ratio test）分析通量随轨道相位的变化。
  • 能谱拟合：测试幂律谱和指数截断幂律谱模型。

3. 主要贡献与关键结果 (Key Contributions & Results)

A. 首次探测到超高能 (UHE) 辐射

• 显著性：
  • WCDA (1.4–30.5 TeV)：探测显著性为 9.2$\sigma$。
  • KM2A (25–267 TeV)：探测显著性为 6.2$\sigma$。
  • >100 TeV 能段：在 KM2A 数据中识别出 16 个 类光子事件，背景估计为 5.1 个，对应 3.8$\sigma$ 的显著性探测。
• 最高能量事件：探测到一个能量估计为 267 ± 61 TeV 的光子，距离源中心仅 0.09°，与探测器的点扩散函数（PSF）一致。
• 能谱特征：
  • 从 1 TeV 到数百 TeV 的能谱可以用单一幂律模型很好地描述：$dN/dE \propto E^{-\Gamma}$，光子指数 $\Gamma = 3.00 \pm 0.05$。
  • 未发现显著的能谱截断（Cutoff），表明粒子加速机制能够持续产生极高能粒子。
  • 积分通量（>1 TeV）约为蟹状星云通量的 3.8%。

B. 轨道调制与能量依赖性

• 低能段 (WCDA, 1.5–30.5 TeV)：辐射覆盖大部分轨道相位，峰值出现在 $\phi \approx 0.6$ 附近（近拱点之后），显著性约为 2.6$\sigma$。
• 中能段 (KM2A, 25–100 TeV)：辐射更集中在 $\phi = 0.3–0.6$ 区间，显著性提升至 3.9$\sigma$。
• 超高能段 (>100 TeV)：由于统计量有限，数据分为两个相位区间。结果显示 >100 TeV 的光子倾向于聚集在 $\phi = 0.6–0.2$（即远拱点附近），这与低能段的分布模式不同，暗示了能量依赖的轨道调制（尽管目前显著性较低，仅为 2.4$\sigma$ 的提示）。

C. 物理机制解释 (Discussion)

• 轻子模型的局限性：纯轻子模型（逆康普顿散射）在 UHE 区域受 Klein-Nishina 效应抑制，且同步辐射冷却限制了电子的最大能量（通常 <100 TeV），难以单独解释 >200 TeV 的光子。
• 强子与混合模型：
  • 强子过程：在轨道相位 $\phi \approx 0.6$（近拱点/盘穿越期），致密天体穿过 Be 星的高密度星周盘，质子 - 质子 ($pp$) 碰撞效率显著提高，产生 $\pi^0 \to \gamma\gamma$。
  • 混合场景：论文提出一个轻子 - 强子混合模型。
    • 低能段（<25 TeV）主要由轻子（电子）通过逆康普顿散射主导。
    • 高能段（>25 TeV，特别是 >100 TeV）主要由强子（质子）主导。
  • PeVatron 候选体：>100 TeV 光子的探测表明 LS I +61◦303 能够将质子加速至 PeV 能量，是银河系内潜在的 PeVatron 候选体。
  • 能量预算：混合模型所需的总注入功率（约 $4.2 \times 10^{33}$erg/s）与脉冲星的自转减慢功率（$\dot{E} \approx 6 \times 10^{33}$ erg/s）相容，解决了纯强子模型面临的能量预算危机。

4. 科学意义 (Significance)

1. 突破能量上限：首次将伽马射线双星的探测能段延伸至 100 TeV 以上，挑战了传统的粒子加速模型，证明了双星系统可以产生极端的超高能辐射。
2. 揭示加速机制：通过观测到能量依赖的轨道调制，为理解双星系统中轻子和强子过程的竞争与协同提供了关键证据，支持了混合辐射机制。
3. PeVatron 候选：确认 LS I +61◦303 具备加速质子至 PeV 能级的能力，为研究银河系宇宙线的起源提供了新的候选源。
4. 多波段协同：结合 LHAASO 的宽能段覆盖与 Fermi-LAT、MAGIC、VERITAS 等历史数据，构建了从 MeV 到 PeV 的完整能谱，为未来多信使天文学研究奠定了基础。

总结

该论文利用 LHAASO 的大视场和高灵敏度，首次确凿地探测到了伽马射线双星 LS I +61◦303 的超高能（>100 TeV）辐射。这一发现不仅扩展了该类天体的能谱范围，还通过轨道调制特征揭示了其内部复杂的粒子加速环境，支持了轻子与强子共同作用的混合模型，并确立了其作为银河系 PeVatron 候选体的地位。