Time-dependent cosmic-ray escape from wind bubbles: hard spectra formation

技术摘要：风泡中随时间变化的宇宙线逃逸：硬谱形成

问题陈述
由年轻大质量恒星团（YMSCs）、沃尔夫-雷特（Wolf–Rayet）星或超爱丁顿吸积致密天体等源头驱动的风驱动泡，是将其加速至 PeV 能级的宇宙线（CRs）极具潜力的候选者（PeVatrons）。虽然在内层风终止冲击波（TS）处通过扩散激波加速（DSA）被广泛认为是主要的加速机制，但粒子从终止冲击波到前向冲击波（进入星际介质的逃逸边界）的随时间变化的传输过程仍缺乏深入研究。以往的研究通常依赖于定常描述或估算传输时间尺度，而非求解完整的随时间变化的传输方程。本研究旨在解决理解风泡的动态膨胀以及终止冲击波与前向冲击波之间随时间变化的距离如何影响逃逸宇宙线能谱这一空白。

方法论
作者将风泡建模为一个球对称系统，其中终止冲击波（TS）和前向冲击波（FS）以随时间变化的速率膨胀（在减速阶段，$R_{ts} \propto t^{2/5}$ 且 $R_{fs} \propto t^{3/5}$）。宇宙线的传输受一维球对称随时间变化的传输方程（福克-普朗克方程）控制，该方程包含了由冲击后风流引起的平流和扩散。

为了求解该方程，作者利用伊藤引理（Itô's lemma）将其转化为一组随机微分方程（SDEs）。这些 SDEs 使用开源框架 CRPropa3.2 的修改版本进行积分。关键方法论特征包括：

• 注入： 粒子在 TS 处以预加速的幂律谱（$dN/dE \propto E^{-2}$）注入，其高能截止能由希拉斯极限（Hillas limit）和加速时间尺度决定。
• 平流： 在下游区域假设径向速度剖面为 $u(r) \propto r^{-2}$，并归一化至冲击后速度。
• 扩散： 研究考察了三种湍流模型：Bohm 型（$\delta=1$）、Kraichnan 型（$\delta=3/2$）和 Kolmogorov 型（$\delta=5/3$）。扩散系数取决于局部磁场强度和湍流功率谱。
• 磁场： 磁场基于冲击压力的比例进行估算。研究测试了下游磁场的各种径向依赖关系（$B \propto r^{-\alpha}$，其中 $\alpha = 0, -1, -2$）以及磁化水平（$\epsilon_B$）。
• 模拟参数： 模拟运行至 5 Myr，覆盖了在超新星主导能量收支之前的风驱动阶段。

主要结果
随时间变化的处理揭示了与标准预期逃逸谱显著不同的现象：

1. 谱硬化： 在特定的物理配置下，逃逸的宇宙线谱可以比注入的 $E^{-2}$ 谱更硬。这种硬化在具有更强能量依赖性的湍流模型（例如 Bohm 扩散）和较弱磁场依赖性的情况下最为显著。
2. 低能抑制： 一个关键发现是逃逸通量中低能粒子的显著抑制。由于前向冲击波（逃逸边界）随时间向远离源的方向移动，具有较小扩散系数的低能粒子无法追上逃逸边界，从而滞留在下游区域。
   • 对于 Bohm 扩散，低能粒子（$E < 1$ TeV）在逃逸通量中几乎完全枯竭。
   • 对于 Kraichnan 扩散，观察到适度的硬化（$\Delta \gamma \approx 0.1-0.2$）和部分低能抑制。
   • 对于 Kolmogorov 扩散，耗尽现象极小，逃逸谱与注入谱非常接近，尽管最大能量仍然较低（$\sim 20$ TeV）。
3. 参数依赖性： 低能耗尽程度和谱硬化程度取决于湍流模型、磁场的径向依赖关系（较弱的依赖关系导致更硬的谱）以及磁化水平。
4. 观测意义： 当传播至地球时（模拟源距 3 kpc），低能抑制现象依然存在，尽管由于长时标下的银河系扩散，硬化现象被部分平滑化。

意义与主张
论文声称，风泡膨胀的随时间变化的特性是塑造观测到的宇宙线谱的关键因素，这独立于能量损失过程。作者强调：

• 硬谱形成： 只要传输是随时间变化的，风泡可以自然地产生硬宇宙线谱（比 $E^{-2}$ 更硬），而无需需要非标准的加速物理。
• 低能截止： 低能粒子在泡内被“捕获”是该随时间变化模型的稳健预测。这可能导致来自泡内的多信使辐射（伽马射线和中微子）出现可观测特征，并影响宇宙线累积的克拉马（grammage）。
• 与超新星遗迹（SNR）的区别： 所得出的硬谱和低能截止与超新星遗迹通常预期的特征不同，这可能有助于解释“膝部”（knee）附近的宇宙线谱或特定恒星团（如 Westerlund 1）的伽马射线观测。

作者总结道，虽然其模型侧重于风驱动阶段，但这些被捕获的低能粒子随后可能被超新星冲击波重新加速，或者在前向冲击波停止时逃逸，从而可能以修正后的谱形状贡献于低能宇宙线群体。未来的工作将结合强子相互作用，以预测随时间变化的伽马射线和中微子发射。