编辑推荐：
专业摘要翻译
微类星体（Microquasars）是我们银河系内高能粒子产生的候选场所。由于它们距离遥远且辐射区域紧凑，限制了直接的观测约束，因此数值模拟对于连接喷流动力学与可能的多信使信号至关重要。本研究对来自相对论性磁流体动力学（relativistic
专业摘要翻译
微类星体（Microquasars）是我们银河系内高能粒子产生的候选场所。由于它们距离遥远且辐射区域紧凑，限制了直接的观测约束，因此数值模拟对于连接喷流动力学与可能的多信使信号至关重要。本研究对来自相对论性磁流体动力学（relativistic magneto-hydrodynamic, RMHD）微类星体喷流的中微子产生过程进行了建模，重点关注加速质子与周围光子场相互作用下的质子-光子（pγ）通道。研究人员承认质子-质子（pp）相互作用在稠密环境中是一个可能的过程，但在本计算中未予建模。研究人员应用了一种射线追踪程序来处理流体动力学模拟的输出，从而从一个固定观测者的视角生成合成中微子图像。研究人员展示了合成的谱图和强度图，并将其与当前及未来中微子探测器的代表性灵敏度进行了比较。
论文解读文章
研究背景、现存问题与研究意义
微类星体是由一颗主序星绕一个致密天体运行的恒星双星系统。物质向致密天体的转移驱动了相对论性喷流，这些喷流从射电到甚高能伽马射线的整个电磁波谱都有辐射，并且可能在强子相互作用场景下产生中微子。尽管已有诸多理论探讨，但由于微类星体距离遥远、辐射区域复杂，对其内部的高能粒子产生和中微子发射过程进行直接的观测限制非常困难。特别是，近期超高能伽马射线观测（如LHAASO的结果）加强了将银河系致密天体喷流系统视为极端粒子加速器（PeVatron）的证据，这进一步激发了人们对强子加速及其伴随中微子产生机制的研究兴趣。然而，此前的研究缺乏将喷流的真实动力学演化与具体的中微子辐射过程紧密结合的、基于时间演化的数值模型。为了弥合喷流动力学理论与多信使观测之间的鸿沟，有必要开展更精细的数值模拟研究。
研究方法与技术路线
为开展此项研究，研究人员构建了一套结合多种数值工具的完整技术路线。首先，利用PLUTO (4.4-patch4) 这一激波捕捉型有限体积RMHD代码，对间歇性的相对论性双极喷流进行三维数值模拟，模拟域聚焦于中微子辐射通常发生的内喷流区域。其次，开发并使用名为nemiss的专用代码，该代码读取RMHD模拟输出的每个网格单元数据，基于明确的质子分布和靶光子场假设，求解质子到中微子的级联过程，计算出每个网格内的中微子辐射率。最后，运用名为rlos (1.15.1) 的相对论性视线成像代码，通过追踪射线穿过四维RMHD模拟数据，考虑相对论性集束效应和时间延迟，从静止观测者视角生成合成中微子图像和光谱。此外，还使用了Veusz (4.2.1) 进行数据可视化。
研究结果
2.1 粒子加速区域、机制与扩散模式：
研究人员将喷流中的加速区域分为三类：A区（内部激波/团块前沿）、B区（再准直/相互作用层）和C区（终端/物质加载区域）。假定在激波区域发生一阶费米加速，并根据不同区域湍流水平的不同，采用了依赖于区域的扩散系数参数化方案，即在激波压缩的团块前沿采用近博姆（near-Bohm）标度，而在下游/茧状区域采用较弱的能量依赖性。
2.2 加速区域中的磁场设定：
磁场值直接取自RMHD模拟输出的每个网格单元。在标记为加速区的网格中，若RMHD计算的磁场低于设定的均分比，则采用一个保守的下限值，即热能量密度乘以一个模型参数，以确保加速区域磁场强度的明确性和可重复性。
2.3 非热质子分布：
非热质子被描述为一个幂律能量分布，谱指数为2，质子能量上限截断为106 GeV。该截断是基于Hillas估计，考虑了内喷流辐射区域有限的驻留时间和能量损失。
2.4 靶光子场的起源：
为了消除歧义，研究人员显式地建模了靶光子场，将其分解为三个组分：伴星辐射场（建模为稀释的黑体谱）、吸积盘/冕成分（可选的热或准热成分）、以及散射/风光子场（由恒星/盘光子在风环境中散射形成）。
3.1 喷流内部产生的中微子：
中微子产生源自质子-光子（pγ）相互作用产生的π介子衰变链。中微子辐射率由稳态π介子密度、π介子衰变时标等参数决定，公式中包含了带电π介子衰变链及后续μ子衰变对中微子产额的贡献。
3.2 团块-风碰撞的中微子辐射率：
为了包含团块与周围光子场碰撞导致的中微子辐射增强，研究人员将实验室系下的静态光子分布通过多普勒因子变换到喷流共动系，以此替换了早期方法中使用的同步辐射光子场。
3.3 靶光子分布与模型中使用的多普勒因子情况：
研究人员采用了幂律形式的靶光子分布，并为相互作用的多普勒因子设定了三种方案：案例(a)恒定多普勒因子、案例(b)基于网格单元的局部多普勒因子、案例(c)正面碰撞近似。案例(d)是案例(b)的一个60个网格的子集，仅作为数值基准。
3.5 质子-光子通道：
详细给出了各向同性光子场中质子-光子相互作用的反应率、π介子注入率、π介子数密度以及最终中微子辐射率的计算公式。
6.1 质子-光子谱与可探测性：
研究结果显示，对于所采用的归一化和5千秒差距的距离假设，模拟的质子-光子中微子分量在探测器相关的能量范围内低于冰立方（IceCube）的近似灵敏度阈值。
6.2 通量归一化、仪器灵敏度和距离依赖的可探测性：
研究提供了一个方便的距离-视界度量，用以评估在不同距离上的可探测性。结果表明，在该模型参数下，模拟的质子-光子分量主要功率分布在IceCube主灵敏度窗口以下的能量，且即使考虑相对论性集束效应，也难以弥合与探测灵敏度之间的差距。
6.3 与先前银河系微类星体中微子通量模型的比较：
通过与文献中其他模型（如SS433暗喷流中的pp相互作用模型）的定量比较，研究人员指出本研究的通量较低，主要原因在于本工作孤立地计算了pγ通道，而其他模型通常假设在稠密介质中存在高效的pp相互作用，后者靶物质密度更高，预期中微子输出也更大。
6.4 伴随的γ射线谱：
研究人员估算了伴随的π0衰变伽马射线产量，其能量通量与中微子分量处于同一量级。但强调这是产生层面的估算，并未包含内部γ-γ吸收和电磁级联过程，因此不能代表逃逸出的GeV-TeV伽马射线谱。
总结与结论
本研究提出了一个时间分辨的模型，用于描述来自相对论性微类星体喷流的质子-光子（pγ）中微子辐射，该模型结合了RMHD模拟与关于加速区域、区域依赖的扩散、磁场处理和靶光子场的明确设定。质子-质子（pp）相互作用在稠密环境中是相关的，此处仅为提供背景而提及，并未在本研究中建模。
在更稠密的环境中，喷流中通过质子-质子相互作用产生的中微子辐射是可能的，但其密度可能过高以至于无法允许π介子衰变的伽马射线逃逸。另一方面，在密度较低的区域，质子-光子相互作用可以产生中微子，同时整体不透明度可能足够低，使得可探测的伽马射线能够逃离系统。
在我们的案例中，对于采用的归一化和5千秒差距的距离，计算得到的质子-光子分量低于冰立方（IceCube）的近似灵敏度参考线，尤其是在探测器相关的能量范围内。有利的观测几何和相对论性集束效应可能会增加表观通量，但当前的比较应被视为一个灵敏度受限的估计，而非发现预测。
尽管预计会产生一个中性π介子伽马射线对应体，但其逃逸的GeV-TeV谱无法在不模拟内部γ-γ吸收和电磁级联的情况下推断出来，因为高能光子可能在源内被吸收并发生级联。根据本文所采用的假设，当前结果表明，在微类星体中，从所建模的pγ通道获得可探测的中微子辐射是具有挑战性的，并且很可能不是产生同时中微子-伽马辐射的可行机制。
未来的工作可能会探索更高分辨率的模拟、更精细的辐射场几何结构以及更详细的特定系统参数集。
