这篇研究论文由Benjamin P. Moster（通讯作者，时任马克斯·普朗克天文研究所及马克斯·普朗克天体物理研究所）、Andrea V. Macciò（马克斯·普朗克天文研究所）、Rachel S. Somerville（空间望远镜科学研究所及约翰斯·霍普金斯大学）、Thorsten Naab（马克斯·普朗克天体物理研究所）以及T. J. Cox（卡内基天文台）共同完成。其研究成果于2011年发表在《英国皇家天文学会月刊》（Monthly Notices of the Royal Astronomical Society）第415卷第4期，具体发表日期为2011年。
该研究属于天体物理学领域，特别是星系形成与演化以及星系相互作用的数值模拟研究。研究的背景在于，尽管宇宙学流体动力学模拟和观测表明旋涡星系通常由五个部分组成：暗物质晕、恒星盘、恒星核球、气体盘和气体晕，但在过去对星系并合（尤其是主并合，即质量相近的星系合并）的数值模拟研究中，热气体晕（hot gaseous halo）的影响通常被忽略。然而，这个热气体晕在星系维里半径内所含的气体质量可达星系总气体质量的80%，是一个巨大的气体储库。通过冷却和吸积过程，热晕可以在并合前后补充冷气体盘。这些冷气体随后可以形成恒星，从而影响并合遗迹（remnant）的形态和动力学。以往的研究往往只关注盘中的初始冷气体，而忽略了持续从热晕冷却补充气体的过程。因此，本研究的主要目标是，首次在盘星系主并合的流体动力学模拟中引入一个旋转的、可冷却的、扩散的热气体晕，并系统地研究其对并合过程中及并合后恒星形成效率、并合遗迹的形态结构及动力学特性的影响。研究旨在验证一个关键假设：在考虑了更真实的、持续供气的热晕后，星系并合是否总是导致恒星形成的显著增强？并合遗迹的性质（如旋转速度、椭率、表面亮度轮廓）会发生何种变化？
该研究的工作流程可分为三个主要阶段：首先，建立和校准包含热晕的孤立星系模型；其次，基于校准后的模型，进行一系列星系主并合模拟；最后，对模拟结果进行分析，重点对比有热晕和无热晕两种情形。
第一阶段：校准热晕模型参数。研究人员使用GADGET-2（一款并行树形SPH代码）进行所有模拟。他们首先构建了一个典型的银河系（Milky-Way）质量（总质量约10^12倍太阳质量）的孤立星系模型，包含暗物质晕、恒星盘、恒星核球、冷气体盘以及新加入的热气体晕。热晕的密度分布采用了观测上常用的β-剖面（β-profile），并假设其处于流体静力学平衡和旋转状态。这里有一个关键的自由参数：热晕的自旋因子α，它定义了热晕的比角动量相对于暗物质晕比角动量的倍数（jhg = α jdm）。为了确定符合实际的α值，研究团队让该孤立星系模型从红移z=1演化到z=0（历时约76亿年），并将模拟结果与两个观测约束进行比较：星系在不同红移下的恒星质量（来自Moster等人2010年的经验关系）和恒星盘尺度长度（来自GEMS巡天数据）。他们测试了α = 1, 2, 4, 8以及无热晕的情形。模拟发现，没有热晕时，孤立星系的恒星形成率（SFR）呈指数衰减，无法匹配观测到的恒星质量增长。当α值太低（如1或2）时，由于热晕角动量过低，冷却气体过于集中于星系中心，导致恒星质量增长过快、盘尺度太小。当α值太高（如8）时，冷却气体分布在过大半径上，导致恒星质量增长不足、盘尺度太大。最终，α=4的模型能够最好地同时匹配恒星质量和盘尺度随时间的观测约束。因此，研究确定α=4为后续所有并合模拟中热晕的基准自旋参数。这部分工作为后续的并合研究提供了物理上可信的初始条件。
第二阶段：执行星系主并合模拟。在确定了热晕模型后，研究团队转向主并合研究。他们采用了Cox等人（2008）使用过的一个典型盘星系模型（称为“G3”模型），其恒星质量约为5×10^10倍太阳质量，并相应地加入了质量为1.1×10^11倍太阳质量、α=4的热气体晕。为了全面探究影响，他们设计了一个多层次的模拟对比方案： 1. 基础对比：对基准G3模型（初始盘气体分数约23%），进行了四组并合模拟，以隔离热晕和星风（stellar winds，一种恒星反馈机制）的影响：(g3g3) 无热晕、无星风；(g3g3w) 无热晕、有星风；(g3g3h) 有热晕、无星风；(g3g3wh) 有热晕、有星风。 2. 气体分数依赖性：对于有热晕和无热晕的情况（均不考虑星风），额外构建了初始盘气体分数分别为0%、40%、80%的星系模型进行并合模拟，以研究初始气体条件的影响。 3. 热晕质量影响：还测试了热晕质量减半（“最大”质量的一半）的情形（g3hx），以考察效应是否随热晕质量线性变化。
所有并合均采用相同的近径向、近共面轨道参数（偏心率0.95，近心距13.6千秒差距，初始分离250千秒差距），其中一个星系的盘与轨道面平行，另一个倾斜30度。每个模拟都运行50亿年，以确保并合遗迹达到动力学的弛豫状态。
第三阶段：数据与结果分析。研究人员对模拟数据进行了详尽分析，主要包括： 1. 恒星形成历史与星暴效率：计算并合过程及孤立演化下的恒星形成率随时间的变化。定义了“星暴效率”（burst efficiency）E = m_burst / m_cold，其中m_burst是因并合额外形成的恒星质量（并合系统新形成恒星质量减去两个孤立演化星系新形成恒星质量之和），m_cold是并合前夕系统中的冷气体总质量。这个效率量化了并合过程将可用冷气体转化为恒星的额外效率。 2. 遗迹的形态学分解：在模拟结束后，使用基于角动量的方法将遗迹中的恒星粒子分解为盘成分和球状成分，计算了核球占总恒星质量的比例（B/T）。 3. 表面密度与三维密度剖面：计算了恒星（分为“老”星和并合期间及之后形成的“新”星）、冷气体、热气体和暗物质的面密度或三维密度分布。 4. 运动学性质分析：通过从1000个随机视角投影遗迹，计算了每个视角下的椭率（ε）以及沿主轴测得的旋转速度与中心速度弥散的比值（V_maj/σ）。通过将结果与各向异性椭圆体模型的理论曲线对比，评估了遗迹的旋转支撑程度和各向异性。
该研究取得了多方面的主要结果。
关于恒星形成，模拟结果清晰地显示，热晕的引入极大地改变了并合前后恒星形成的图景。对于孤立星系，包含热晕的模型能够维持一个近乎恒定的恒星形成率（约每年5个太阳质量），因为热晕的持续冷却和吸积不断补充着冷气体盘。而无热晕的孤立模型，其恒星形成率则迅速呈指数衰减。在并合过程中，热晕的存在显著降低了星暴的“相对”增强幅度。在无热晕的基准模拟（g3g3）中，并合峰值恒星形成率可达孤立星系之和的约30倍，星暴效率E高达0.68。然而，在包含热晕的模拟（g3g3h）中，峰值恒星形成率仅为孤立星系之和的约5倍，星暴效率降至0.51。更引人注目的是，在某些情况下（如g3g3h模拟），并合遗迹在并合结束后的总恒星质量，甚至低于两个原星系在孤立演化下达到的恒星质量之和。这是因为并合事件通过两个机制“损害”了后续的气体吸积和恒星形成能力：第一，并合过程将轨道角动量转移给了热晕气体，使其比角动量增加，从而抬高了气体冷却和沉降到星系中心的离心势垒；第二，并合时两个热晕碰撞产生的激波加热了气体，延长了其冷却时标。这两个效应共同导致并合遗迹从热晕吸积冷气体的效率低于孤立星系，因此并合后的恒星形成率反而可能低于孤立演化的情况。这一发现挑战了半解析星系形成模型中“并合必然增强恒星形成”的普遍假设。
关于遗迹形态，研究发现，包含热晕会使并合后形成的恒星盘质量更大，从而降低了遗迹的核球占比（B/T）。例如，在基准模型中，无热晕遗迹的最终B/T稳定在~0.9，而有热晕遗迹的B/T则从并合后的~0.9逐渐下降到~0.8左右，原因是更多吸积的冷气体形成了一个更显著的恒星盘。不过，即使初始气体分数高达80%，所有主并合遗迹在本质上仍是核球主导的（B/T > 0.8）。在表面亮度分布方面，所有遗迹在大尺度（> 1千秒差距）上都符合经典的R^{1/4}（de Vaucouleurs）轮廓。但在小尺度（< 1千秒差距），有热晕的遗迹显示出更高的表面密度，形成了更密集的恒星核。这不仅是由于新形成的恒星盘集中在中心，还因为新恒星加深了中心势阱，将老恒星也拖曳向中心，从而在中心产生了类似“额外光”（extra light）的陡峭轮廓，与一些观测相符。
关于运动学，各向异性图分析表明，热晕的加入增加了快速旋转体（fast rotators）的比例。无热晕、无星风的并合遗迹旋转缓慢（V_maj/σ最大值约0.25），主要靠速度弥散支撑，各向异性较强。而包含热晕的遗迹显示出更高的旋转速度，V_maj/σ值更大，其分布更接近各向同性旋转体模型。这种更强的旋转支撑主要归功于并合后从热晕吸积形成的大质量、动力学上较冷的恒星盘。
此外，研究还发现星风反馈会延长星暴的持续时间，因为它将部分气体吹离星系中心，这些气体随后再冷却回落。关于初始气体分数的依赖性，在无热晕模型中，星暴效率随初始气体分数轻微下降；但在有热晕模型中，效率随初始气体分数增加，因为高气体分数模型在并合时保留了更多初始的、高度集中的冷气体盘可供消耗。当同时考虑星风和热晕时，极高的初始气体分数（80%）反而导致效率最低，因为星风阻止了并合前形成一个足够质量的恒星棒来有效地将气体驱入中心触发星暴。
该研究的结论具有重要的科学意义。首先，它证明了在星系并合模拟中考虑热气体晕是至关重要的。热晕不仅作为一个巨大的气体储库，持续影响星系并合前、中、后的恒星形成历史，而且其动力学状态（角动量、温度）会因并合事件而发生改变，从而对并合遗迹的长期演化产生深远影响。其次，研究挑战了并合必然导致净恒星质量增加的简单图景，指出在某些条件下（特别是当热晕获得额外角动量并被加热后），并合甚至可能抑制长期的恒星形成，导致遗迹质量低于孤立演化之和。这要求半解析模型等理论框架必须修正其关于并合触发星暴的简化处理。第三，研究表明，热晕的存在有助于形成更多具有显著旋转、中心密集、且拥有R^{1/4}轮廓小尺度核的椭圆星系，这为理解观测上发现的、数量众多的快速旋转椭圆星系（特别是中等质量星系）的起源提供了新的物理机制。以往模拟需要很高的初始盘气体分数才能产生快速旋转遗迹，而现在，通过热晕的持续供气，即使初始盘气体分数适中，也能在并合后形成显著的盘结构，这更符合实际星系的观测属性。
本研究的亮点突出体现在以下几个方面：创新性：这是首次在盘星系主并合的流体动力学模拟中系统性地引入并研究一个旋转、冷却的热气体晕的效应，填补了该领域的研究空白。方法的严谨性：研究没有随意设定热晕参数，而是先通过孤立星系演化匹配观测约束来校准关键参数（自旋因子α），确保了初始条件的物理真实性。系统性对比：通过精心设计的多组模拟（有无热晕、有无星风、不同气体分数、不同热晕质量），清晰而全面地分离并揭示了热晕的独立影响及其与其他物理过程（如反馈）的相互作用。重要发现：提出了“并合可能抑制长期恒星形成”这一反直觉但物理上合理的新见解，并揭示了热晕角动量增加和加热是导致此现象的关键机制，对星系形成理论模型有直接修正意义。广泛的关联性：研究结果不仅涉及恒星形成，还关联到并合遗迹的形态结构、运动学分类、以及椭圆星系的整体形成图景，具有广泛的学术价值。