利用联合喷流-风活动星系核反馈解决星系团冷却流问题

《SCIENCE ADVANCES》：Solving the cooling flow problem with combined jet-wind AGN feedback

　　

活动星系核（AGN）反馈被广泛视为解决星系团中长期存在的冷却流问题的最有前景方案，然而先前模型所设定的喷流性质与吸积物理不一致。研究人员利用MACER框架，对一个无并合历史、处于弛豫态且其他性质类似于英仙座星系团的星系团进行了理想流体动力学模拟，该模拟同时纳入了喷流和风，其性质受黑洞吸积的广义相对论磁流体动力学（GRMHD）模拟及观测约束。这种联合反馈成功复现了关键可观测量，包括冷气体质、恒星形成率、热力学径向轮廓以及黑洞增长，而仅喷流或仅风的模型均告失败。其成功源于喷流-风剪切所驱动的湍流，该湍流增强了动能向热能的转化，相对于仅风和仅喷流情形，加热效率分别提升了3倍和6倍。

该研究聚焦于星系团中经典的冷却流问题。大量星系团展现出冷核结构，中心气体具有高密度、低熵的特征，冷却时间尺度约为十亿年。在此条件下，理论预期将产生强烈的冷却流入流，速率可达每秒数百至数千倍太阳质量，同时伴有大量冷气体在星系团核心累积以及中心最亮星系中极高的恒星形成率。然而，多波段观测发现的冷气体质量和恒星形成率远低于理论预测，这一显著差异即被称为"冷却流问题"。为解释这一矛盾，研究者曾提出多种机制，包括恒星与超新星反馈、宇宙线、磁场以及热传导等，其中活动星系核反馈因其与抑制冷却、X射线腔及低熵核心中多相气体的强观测关联而被视为主要解决方案。然而，现有模拟中的活动星系核反馈模型存在重要物理缺陷：分辨率不足导致无法可靠计算邦迪半径处的黑洞吸射率；喷流或风的参数常被当作自由变量处理，与黑洞吸积与喷流形成理论不一致；尽管理论和观测工作已证明喷流和风通常同时存在，但大多数研究未能同时自洽地纳入两者；此外，许多理想化冷却流模拟仅关注解决冷却流问题本身，而未检验其与其他观测约束（如黑洞增长）的一致性。

研究人员基于MACER框架开展了理想化高分辨率三维流体动力学模拟，研究对象为一个无并合历史、处于弛豫态且性质类似于英仙座星系团的星系团。该框架具有多项 distinctive 特征：模拟域的内边界小于黑洞吸积流的外边界，使得结合内边界质量通量与黑洞吸积理论能够获得可靠的黑洞质量吸积率估计，这对研究活动星系核反馈至关重要，因为吸积率决定了包括喷流在内的各活动星系核输出分量的功率；采用了先进的吸积物理，将黑洞吸积按是否低于或高于爱丁顿率的2%分为热模式和冷模式，在每个模式下考虑所有活动星系核输出，即辐射、喷流和风，且这些输出的参数值并非作为自由参数处理，而是直接采用或基于黑洞吸积的理论与观测研究进行外推。

为评估同时纳入风和喷流的重要性，研究人员设计了三组模拟实验：JetWind（喷流+热风）、JetOnly（仅喷流）和WindOnly（仅热风），在相同质量吸积率下保持总外流功率一致。结果显示，在前2亿年，由于初始条件仅包含热气体，几乎没有冷气体或恒星形成，活动星系核活动微弱。约2亿年后，冷却流发展，通过辐射冷却形成冷气体，恒星形成率相应上升，冷气体的落轴触发各模型中强烈的活动星系核活动与反馈。JetWind模型在抑制冷气体和恒星形成方面最为有效，其冷气体质量和恒星形成率与观测一致；WindOnly模型虽也与观测相符，但处于统计样本的极端情况；JetOnly模型则显著高估了冷气体和恒星形成率。在黑洞质量增长方面，WindOnly和JetOnly模型的黑洞质量在15亿年演化后接近或超过10¹⁰ M_☉，远超英仙座星系团观测值范围（3.4×10⁸ M_☉至1.2×10⁹ M_☉），而JetWind模型的黑洞质量始终保持在合理水平，表明其反馈调节更为有效。

在热力学径向轮廓方面，JetOnly模型在约100千秒差距内熵和温度明显低于观测值、密度高于观测值，表明仅靠喷流无法充分加热星系团核心的星系际介质；WindOnly模型在大部分模拟时间内熵显著高于观测值，因温度过高、密度过低，且100千秒差距以外区域几乎不受影响；唯有JetWind模型在整个模拟域内复现了与冷核星系团观测一致的熵、密度和温度径向轮廓。此外，JetWind模型的冷却时间与自由落体时间之比与观测一致，而JetOnly模型低于观测范围、WindOnly模型高于观测范围。在湍流特性方面，JetWind模型在30千秒差距内的速度功率谱与Hα速度结构函数的短波部分斜率和形状一致，平均速度弥散150至200千米每秒与X射线观测定量匹配，而另两个模型未能复现观测结果。

JetWind模型成功解决冷却流问题的关键在于其最高的能量转化效率。尽管JetWind模型的活动星系核功率在三者中最低，但其冷却也最弱。三者的总活动星系核能量与星系际介质冷却能量分别为：JetWind为1.9×10⁶²尔格和1.8×10⁶²尔格，WindOnly为4.2×10⁶²尔格和6.6×10⁶²尔格，JetOnly为1.4×10⁶³尔格和2.9×10⁶³尔格。仅JetWind中时间积分的活动星系核功率超过累积的星系际介质冷却损失。这一结果源于系统的自调节特性：总冷却与加热能量之比越高，冷气体越多，进而增强吸积和活动星系核加热，降低冷却-加热比。在此反馈环路中，将活动星系核功率转化为星系际介质热能的效率是关键决定因素。JetWind的湍流耗散能量为4.6×10⁶¹尔格，转化效率29.6%；WindOnly为3.9×10⁶¹尔格，效率9.3%；JetOnly为6.0×10⁶¹尔格，效率5.2%。计入激波加热后，总转化效率分别为62.6%、17.9%和10.1%。JetWind的高效能源于其最大的湍流振幅，而湍流主要由喷流、风与环境星系际介质之间的剪切通过开尔文-亥姆霍兹不稳定性产生 induction 产生，级联至更小尺度并将动能转化为热能。

研究采用的主要关键技术方法包括：基于MACER3D框架的三维流体动力学模拟，该框架采用球坐标系求解标准欧拉形式流体方程，计算域从内边界100秒差距延伸至外边界2兆秒差距，峰值空间分辨率达约5秒差距；采用低扩散HLLC黎曼求解器和二阶龙格-库塔时间积分器的网格化Athena++代码实现；根据黑洞吸积理论，将吸积分为热模式与冷模式，热模式下同时纳入喷流和风，参数取自三维GRMHD模拟或基于其外推；喷流参数源自快速旋转黑洞的磁化吸积盘（MAD）模式的GRMHD模拟，包括总喷流能量流（0.9?_BHc²）、喷流速度（0.1c，受英仙座星系团射电观测约束）、质量流率（0.35?_BH）和半张角（约7.5°）；热风性质同样基于GRMHD模拟，质量流和速度由吸积流理论确定，角分布限制在10°至50°和130°至170°范围；冷风的参数则基于活动星系核光度与外流性质的观测统计关系确定。

研究结果可归纳为以下几个部分：

模拟设计与模型比较：基于MACER框架对英仙座型星系团进行理想化高分辨率三维流体动力学模拟。JetWind模型同时纳入喷流和风，JetOnly和WindOnly模型分别禁用热风和喷流，保持总反馈功率相当。设置旨在检验多通道反馈配置对抑制星系际介质冷却的重要性。

冷气体质量与恒星形成率演化：JetWind模型产生最低的冷气体质量与恒星形成率，WindOnly次之，JetOnly最高。JetWind和WindOnly与观测上限一致，JetOnly显著高估。

黑洞质量增长：WindOnly和JetOnly模型黑洞质量过度增长至近10¹⁰ M_☉以上，远超观测；JetWind模型始终保持合理水平，体现有效反馈调节。

热力学径向轮廓：仅JetWind模型在整个模拟域复现与冷核星系团观测一致的熵、密度和温度轮廓。JetOnly核心加热不足，WindOnly核心加热过度且影响范围局限。

冷却时间与自由落体时间比：仅JetWind模型维持与观测一致的t_cool/t_ff值，JetOnly低于、WindOnly高于观测范围。

湍流特性：JetWind模型的速度功率谱和速度弥散与Hα与X射线观测匹配，另两个模型不符。

能量转化效率分析：JetWind的湍流耗散与激波加热总效率达62.6%，分别为WindOnly和JetOnly的约3倍和6倍，根源在于喷流-风剪切驱动的强湍流。

模拟与模型的比较：MACER框架在吸积-反馈模式、喷流开角、质量和速度、风的性质、黑洞吸积率计算、外流与星系际介质/星际介质的相互作用方式及分辨率等方面，与现有理想化模拟和宇宙学模拟存在系统性差异，其多通道物理约束的建模方式更为自洽。

讨论部分，研究人员指出该模拟为理想化设置，未考虑并合和宇宙学流入等过程，与真实星系的比较需谨慎。尽管如此，一些研究表明这些过程在抑制冷却流中仅起次要作用。活动星系核反馈在缓解冷却流问题及更广泛调节星系形成与演化中起基础性作用，喷流-风耦合为将活动星系核动能转化为对宿主系统有效反馈提供了高效机制，尤其在大质量晕中具有持续活动星系核活动的场景下，这一过程值得在未来理论和数值研究中明确考虑。研究结论可总结为：研究人员利用MACER框架，首次在理想化高分辨率三维流体动力学模拟中同时自洽地纳入受GRMHD模拟和观测约束的喷流与风，发现联合喷流-风反馈成功解决星系团冷却流问题，复现冷气体质量、恒星形成率、热力学径向轮廓和黑洞增长等关键观测约束，而单一组件模型均告失败；其物理机制在于喷流-风剪切驱动的开尔文-亥姆霍兹不稳定性产生强湍流，显著增强动能向热能的转化，使加热效率相对于仅风和仅喷流情形分别提升3倍和6倍，为理解活动星系核反馈在星系团尺度上的能量沉积机制提供了新的物理图像。该论文发表于《SCIENCE ADVANCES》。