A HINSA view of cosmic-ray ionization in IC 348 and NGC 1333: evidence for a strong low-energy cosmic-ray disparity

技术摘要：通过 HINSA 研究 IC 348 与 NGC 1333 中的宇宙线电离现象

问题陈述
宇宙线电离率（CRIR, $\zeta_{ion}$）是控制分子云化学与动力学演化的基本参数，特别是通过调节磁场与气体之间的耦合。虽然近期的观测揭示了在大质量恒星形成区及原恒星附近存在升高的 CRIR 值，但低能宇宙线（LECRs, $E < 1$ GeV）的物理起源和加速机制仍不确定。此外，在附近的英仙座分子云复合物中存在一项差异：尽管 NGC 1333 的恒星形成活动更高，但其弥散亚 GeV $\gamma$ 射线残余呈现亏损；而其邻近的 IC 348 则表现出正残余以及更硬的 CR 能谱。这些区域内局部 CR 源在多大程度上压倒了银河系背景，以及两者之间差异的本质，仍需进一步澄清。

研究方法
本研究利用 H I 窄自吸收（HINSA）技术，对 IC 348 和 NGC 1333 的 CRIR 进行了新的估算。该方法包括：

1. 观测： 利用来自 FAST 望远镜（针对 NGC 1333）的 H I 数据和 GALFA-H I 巡天数据（针对 IC 348），结合来自 COMPLETE 巡天和 JCMT 的存档 $^{13}$CO ($J=1-0$ 及 $3-2$) 数据，以及来自 Herschel 银带巡天的尘埃连续谱数据。
2. HINSA 分析： 通过对相对于 $^{13}$CO 速度场的 H I 吸收特征进行二阶导数拟合，推导出原子氢柱密度 $N(\text{HINSA})$。在 $N(\text{H}_2) > 3 \times 10^{21}$ cm$^{-2}$ 的区域计算 HINSA 比例（$f_{\text{H-INSA}} = N(\text{HINSA})/N(\text{H}_2)$），以最小化光解作用的影响。
3. 密度约束： 使用受 $^{13}$CO 线强度约束的非局部热平衡（non-LTE）辐射传输模型（RADEX）以及马尔可夫链蒙特卡洛（MCMC）采样，来估算气体体密度 $n(\text{H}_2)$。
4. CRIR 计算： 从 H$_2$ 在尘埃颗粒上的形成与宇宙线破坏（包括电离、解离及 CR 诱导的 UV 过程）之间的稳态平衡中推导出 $\zeta_{ion}$。
5. 能谱建模： 在有限板片衰减框架内对 LECR 能谱进行建模。该模型结合了从 Fermi $\gamma$ 射线观测中推断的高能 CR 能谱约束，通过改变低能能谱斜率 ($\alpha$) 和动量断点 ($p_0$) 来拟合观测到的 $\zeta_{ion}$ 剖面。

主要结果

• CRIR 差异： 研究发现两朵分子云之间存在显著的 CRIR 差异。在 IC 348 中，加权平均 $\zeta_{ion}$ 范围为 $(5.35 \pm 1.49) \times 10^{-17}$ s$^{-1}$ 至 $(2.33 \pm 0.55) \times 10^{-16}$ s$^{-1}$。相比之下，NGC 1333 的数值低了约一个数量级，范围为 $(3.94 \pm 2.46) \times 10^{-18}$ s$^{-1}$ 至 $(1.68 \pm 0.25) \times 10^{-17}$ s$^{-1}$。
• 柱密度趋势： 在两朵云中，$\zeta_{ion}$ 均随 H$_2$ 柱密度的增加而降低，这与缺乏内部加速源的衰减模型一致。
• 能谱建模： 最优拟合的低能 CR 能谱能够重现观测到的剖面。模型表明，约 10 倍的电离率差异是由低能（$\sim$10 MeV）质子通量的相应差异驱动的。IC 348 的能谱比 NGC 1333 显著更硬且更强，而后者在 GeV 能段相对于银河系背景表现出亏损。
• 密度与结构分析： 作者研究了来自气体密度估算和云结构（通过 3D-PDR 模拟）的潜在偏差。虽然 $n(\text{H}_2)$ 的不确定性和分形云结构可能会使 $\zeta_{ion}$ 的绝对值产生偏差（在低 $N(\text{H}_2)$ 时可能低估，在高 $N(\text{H}_2)$ 时可能高估），但分析得出结论，这些因素无法解释两区域之间一个数量级的差异。

意义与主张
本文主张，观测到的 CRIR 差异是由于云本身的物理性质所致，而非演化时间尺度或测量偏差造成的产物。

• 外部加速源： IC 348 中极高的 CRIR 结合其较硬的 $\gamma$ 射线能谱，表明存在位于致密云核外部的局部 LECR 加速源，这可能不同于主导 NGC 1333 的原恒星活动。作者提出，在演化程度更高、嵌入程度较低的 IC 348 星团（包含更多 B 型星）中，恒星风终止冲击波是一个可能的来源。
• 原恒星贡献： 结果表明，虽然原恒星喷流可能会在小尺度上增强 CRIR（如在 NGC 1333 的制图中所见），但它们并不驱动这些环境中的大尺度 CRIR，尤其是考虑到 NGC 1333 拥有更高的恒星形成活动却表现出较低的全局 CRIR。
• 对恒星形成的影响： IC 348 全局较高的 CRIR 意味着更高的电离度，从而导致更强的气体-磁场耦合。作者指出，这可能会增强核心坍缩过程中的磁制动，这可能解释了为何观测到 IC 348 中的原恒星盘平均而言比 NGC 1333 中的盘更小。
• 方法论验证： 本研究验证了 HINSA 技术作为探测低质量恒星形成区 CRIR 的手段，并强调了未来需要通过 3D 建模来精确绝对值，同时确认相对差异的稳健性。