木星及更远区域弓激波处的相对论电子加速

《Nature》：Relativistic electron acceleration at the bow shock of Jupiter and beyond

　　

无碰撞激波在天体空间中普遍存在，被广泛认为是宇宙射线加速的主要场所。主流机制——扩散激波加速（Diffusive Shock Acceleration, DSA）要求粒子反复穿越激波前沿，每次穿越获得能量。最大可实现的能量从根本上受希尔拉斯极限（Hillas limit）的约束，该极限将加速器的物理尺度与最大粒子能量关联起来。然而，加速器位置直接观测约束的匮乏限制了研究人员预测大多数天体物理系统中最大粒子能量的能力。在此，研究人员利用美国国家航空航天局（NASA）朱诺号（Juno）航天器的数据，展示了木星弓激波上游相对论电子加速（≥1 MeV）的直接证据，其能量由大尺度前激波瞬变（foreshock transient）提供。研究人员利用这些观测以及互补的太陽系观测，提出了希尔拉斯极限的一个普适标度律，该标度律经验性地将瞬变的可观测尺度与最大粒子能量联系起来。将此标度律应用于从行星弓激波到原恒星喷流和超新星遗迹等多种环境，得到了一个关于最大可获取粒子能量的简单模型，其能量范围从MeV量级到约几十GeV，再到约几十TeV，这为约束天体物理激波处的最大宇宙射线能量提供了一种基于观测的方法。

无碰撞激波（collisionless shock）是宇宙中普遍存在的结构，被广泛认为是粒子加速至相对论能量并贡献宇宙射线（Cosmic Ray, CR）群体的主要场所。扩散激波加速（Diffusive Shock Acceleration, DSA）作为主导机制，描述了粒子如何通过反复穿越激波前沿获得能量。然而，长期存在的"注入问题"（injection problem）指出，DSA仅对能量足够高、能够跑在激波前面的粒子有效，这一过程取决于激波的倾角和强度，而相关机制在所有参数范围内尚未完全明晰。

有前景的解决方案存在于前激波（"precursor"，天体物理学界的术语）的动态环境中，该区域在激波斜交或准平行几何条件下形成于激波上游，即激波法向与环境磁场之间的夹角θ_Bn ? 45°。在此区域内，被称为前激波瞬变（foreshock transients）的大尺度结构能够高效地将低能超热粒子加速至相对论速度。近期在地球弓激波的观测表明，这些瞬变可通过强化激波加速、俯仰角散射（pitch angle scattering）与几何捕获的协同作用，将电子加速至约1 MeV。由此产生的粒子谱已被证实可用E^?1.5幂律很好地描述，这归因于强激波处非相对论粒子经历的DSA，凸显了该机制的有效性，并将准平行激波定位为卓越的粒子加速器。关键的是，这些瞬变（在日球物理学界被称为热流异常（Hot Flow Anomalies, HFAs）、前激波气泡（foreshock bubbles）和自发HFAs（SHFAs）等）是无碰撞激波的基本属性，在不同等离子体环境中形成，同时随宿主系统及其前驱体和前激波的特性而标度化。

这些瞬变过程的普遍性已被整个太阳系的观测所证实，前激波瞬变已在包括水星、金星、火星、地球、木星和土星在内的多颗行星上被识别。研究表明，这些瞬变的物理尺度与主行星弓激波的大小直接相关。应用等价于希尔拉斯极限的原理，该极限将加速器尺度与最大粒子能量关联起来，这种标度化暗示了激波系统全局尺度与其可产生最大粒子能量之间的潜在联系，这一假设进一步得到了动力学模拟的支持，这些模拟在各种参数下持续再现了这些结构。

在此项发表於《Nature》的研究中，研究人员展示了木星弓激波上游相对论电子（>1 MeV）的直接证据，确凿地将它们的加速与前激波瞬变联系起来。研究人员证明了观测到的能量与标度律预测一致，该标度律将系统尺度（S）与前激波瞬变中的加速区域（L）以及随后的能量上限（E_max）联系起来。通过多行星数据验证此标度律，研究人员将分析扩展至原恒星喷流和超新星遗迹等天体物理对象。

研究人员开展此项研究的背景在于，目前存在观测约束不足的问题：尽管希尔拉斯极限理论上有助于预测最大粒子能量，但加速器位置的直接观测约束稀缺，限制了预测能力。此外，DSA的注入问题尚未完全解决，需要探索前激波瞬变作为高效加速器的可能性。研究人员的目的是建立基于观测的框架，将行星尺度的就地观测与天体物理尺度的现象联系起来。

研究人员得出以下主要结论：首先，朱诺号观测证实木星弓激波上游的 foreshock transient 可将电子加速至≥1 MeV，其能谱幂律指数P ≈ ?1.85 ± 0.2，介于非相对论DSA的?1.5极限与相对论极限?2之间；其次，建立了激波驻点距离（S）与加速区域尺度（L）之间的幂律标度关系；第三，将此标度律外推至HH 211原恒星喷流、SN 1987A和SN 1006超新星遗迹，预测其最大粒子能量可达约几十TeV量级，其中SN 1006的预测与观测到的约100 TeV一致。研究的重要意义在于提供了首个基于行星就地观测、可外推至天体物理尺度的宇宙射线能量上限约束框架。

研究人员用到的主要关键技术方法包括：利用NASA朱诺号航天器上JEDI（Jupiter Energetic Particle Detector Instrument）和JADE（Jovian Auroral Distributions Experiment）仪器进行就地等离子体和粒子测量；通过最小方差分析（Minimum Variance Analysis, MVA）确定瞬变结构边界；采用单航天器时序法结合太阳风速度估算瞬变空间尺度；基于非参数Bootstrap分析确定能谱拟合的95%置信区间；构建激波系统尺度（S）与加速区域（L）的经验幂律标度关系；结合希尔拉斯极限与Bohm扩散假设推导最大能量的解析表达式；将行星标度律外推至天体物理系统时使用95%预测区间上限。

**木星弓激波处的相对论电子**

2023年10月1日，朱诺号航天器在 inbound 轨道穿越木星弓激波时，于日木轨道坐标系（Jupiter Solar Orbital, JSO）约[+12.9, +60.5, ?59.3]木星半径处获得了观测数据。在约11:00至13:00 UTC期间，数据记录了一系列局域扰动，包括离子密度和磁场的变化以及高能电子的存在。该活动由两个 distinct 的前激波瞬变事件组成，间隔约1小时：第一个为发育不充分的小尺度瞬变，第二个为更大、信号更清晰的 event。

聚焦于更显著的第二个事件，最强电子增强将瞬变持续时间约束在约12:30–12:50 UTC。单航天器分析得出其特征尺度长度L ≈ 3 × 10⁵ km，相对于行星激波的传播速度v ≈ 320 km·s^?1。该瞬变的特征包括：太阳风的扰动、明显的等离子体加热、局域密度降低以及具有局域极小值（核）和压缩边缘的磁场结构，这些都是前激波瞬变的明确特征。其尺度达几个木星半径，与基于地球和其他行星环境观测的预期一致。值得注意的是，观测到的瞬变期间峰值电子能量比随后的弓激波穿越处高一个数量级，支持了最高效的粒子加速可在上游前激波区域而非仅在主行星弓激波前沿本身发生的观点。

图2展示了JADE（低能）和JEDI（高能）仪器在最强通量期间的联合电子能谱，覆盖从eV到1 MeV的范围。为隔离加速电子群体以与理论直接比较，研究人员减除了环境谱信号。前激波瞬变的能谱在所有能量道次上均显示明显超过3σ背景水平的信号。环境太阳风背景减除后，高能部分电子能谱被P = ?1.5的幂律良好限定，但在更高能量处谱变软趋向?2，最终得到P ≈ ?1.85 ± 0.2的斜率。这是与近期地球前激波瞬变发现一致的关键结果，证明了这种高效加速机制可延伸至木星弓激波环境的MeV能量。相邻弓激波穿越的谱虽显示增至几十keV的强度，但缺乏更高能量的加速，表明高能尾巴仅存在于前激波瞬变观测中。

**从行星尺度到天体物理尺度**

为推广研究发现，研究依托于关键前提：行星激波环境中最高效的粒子加速主要发生在大尺度前激波瞬变内。在此模型中，粒子并非如经典DSA那样严格限制于垂直于激波前沿的加速，而是在这些沿激波前沿传播的结构内被加速。关键的是，加速区域（L）的尺度不是任意的局域参数，而是与激波驻点距离（S）这一全局系统尺度内在关联。这一标度化提供了希尔拉斯极限所需的关键物理尺度（L），从而建立了激波全局尺度与最大粒子能量（E_max）之间的直接经验约束关系。太阳系丰富的就地观测提供了独特机会来确定这一经验关系。

在此框架内，研究人员构建了从行星就地观测得出的经验幂律标度，将激波全局尺度（S）与加速区域（L）联系起来（图3a）。应用希尔拉斯极限并结合此约束尺度和局域激波参数，得到最大粒子能量，其与地球、土星以及本项朱诺号木星新观测结果显示出显著一致性，验证了希尔拉斯极限在太阳系尺度内的适用性。这一在行星尺度的成功验证为外推至天体物理激波奠定了基础，尽管这必然涉及超出直接观测约束的假设。研究人员将此标度关系应用于具有高阿尔芬马赫数（Alfvén Mach number, M_A）和非相对论激波的原恒星喷流HH 211及超新星遗迹SN 1987A和SN 1006。特别地，SN 1006作为已知将粒子加速至几十TeV的γ射线明亮遗迹模型，为模型预测提供了有界检验。

**讨论**

朱诺号观测与早期地球发现共同标志着对行星弓激波处粒子加速理解的转变。研究人员直接捕获了并非仅在激波前沿、而在其扩展激波区域——始于上游前激波——被加速的相对论电子，其中前激波瞬变驱动的电子强度比相邻激波穿越处高一个数量级。瞬变的巨大尺度（跨越数个木星半径）及其至MeV能量的干净幂律谱，证实了大尺度前激波瞬变在强散射条件下提供了高效加速的理想条件。这促使研究人员重新审视经典DSA图像，将其扩展至涵盖沿扩展激波区域的加速。

为超越木星进行推广，研究人员汇编了太阳系内的前激波和激波观测，建立了激波驻点距离S与加速区域尺度L之间的清晰标度。在地球（M_A ~ 4）和木星（M_A ~ 20），前激波扰动的尺度由上游马赫数设定。将此经验标度律外推至具有大数个数量级马赫数和全局系统尺度的原恒星喷流（HH 211）和超新星遗迹（SN 1987A和SN 1006），预测加速区域为10⁸–10¹⁰ km量级。虽然这些尺度无法在天体物理距离直接测量，但支持此推断的物理论述在方法部分有详细说明。此方法最终提出了无碰撞激波全局几何与设定最大粒子能量的瞬变驱动加速之间的直接联系。使用该模型，在乐观但统计上合理的场景下，SN 1987A等系统可达到?10 TeV量级的最大粒子能量，而SN 1006的预测与从?10 TeV γ射线观测推断出的?100 TeV最大能量匹配，支持了所提出的框架。

研究强调，此工作提供了将就地行星观测与大规模天体物理现象联系起来的框架。模型的标度律直接来源于数十年行星数据。朱诺号观测显示跨越数个行星半径的前激波瞬变作为高效粒子加速器，可产生至MeV能量的干净幂律谱。尽管这一观测证据有力，但外推至天体物理激波必然仍是尝试性的。因此，所呈现的框架代表了通过统一的、由观测驱动的方法来约束粒子加速上限的尝试。此外，模型假设Bohm扩散能很好描述环境（即磁涨落强于背景场，δB/B₀ ? 1），且电子有足够时间加速至最大可获取能量。因此，导出的能量极限应被解释为稳健但相对保守的边界。这种有意选择的简洁方法使研究人员能直接从标度律和基本原理进行预测。更复杂的模型可能产生更高或更精确的最大能量，但计算要求高且需要不可获得的天体物理系统上游条件的先验知识。

**研究结论**

研究人员呈现了木星弓激波上游与前激波瞬变相关的相对论电子的直接证据，证明了这些激波产生的瞬变作为高效粒子加速器的作用——这一过程先前在地球确立，现首次在木星展示。观测结果与瞬变所规定空间尺度下强散射环境中激波加速的预测高度一致。利用这些观测与其他行星环境的历史数据，研究人员将标度关系外推至原恒星喷流HH 211和超新星遗迹SN 1987A及SN 1006。应用这些非相对论激波的已知参数空间表明，这些基本过程可将粒子加速至约几十TeV，这一能量范围可能对宇宙射线谱中观测到的电子截断有贡献。最重要的是，SN 1006的预测最大能量与观测一致，为模型提供了直接支持。在此基础上，下一步重要的工作是评估所呈框架在其他天体物理目标上的应用，并确定各种天体物理加速器对宇宙射线通量的相对贡献，这需要跨学科的扩展模拟和观测努力。