技术摘要：利用 MUSES 计算引擎研究重离子碰撞中的 QCD 物质

问题陈述

热且致密的量子色动力学（QCD）物质的状态方程（EoS）对于描述所有碰撞能量下的重离子碰撞至关重要。然而，确定 EoS 面临两个主要挑战：

1. 理论局限性： 由于费米子符号问题，第一性原理格点 QCD（LQCD）计算被限制在消失或极小的实重子化学势（$\mu_B$）情况下。虽然 LQCD 在高温度和低 $\mu_B$ 下提供了可靠的数据，但对于更高密度的情况，则需要唯象模型。确保这些模型以热力学一致的方式符合 LQCD 约束仍然是一个重大障碍。
2. 动力学复杂性： 重离子碰撞本质上是动力学的、非平衡态的系统。将实验数据与理论进行比较需要流体动力学模拟，这些模拟需要演化多种守恒荷（重子数、奇异性和电荷）。这些模拟要求的 EoS 输入不仅要满足热力学一致性，还要具备可逆性（将 $(T, \mu_X)$ 转换为自然流体动力学变量，如熵密度 $s$ 和电荷密度 $n_X$），并能够处理包括潜在临界点（CP）和一阶相变在内的相变过程。

此外，将来自重离子碰撞的约束与中子星联系起来，需要一个能够处理四维相图 $(T, \mu_B, \mu_S, \mu_Q)$ 的统一框架，因为尽管这些系统可能具有相似的重子密度，但它们探索的是该空间的不同区域。

方法论

作者使用了 MUSES（模块化统一状态方程求解器）计算引擎，特别是全新的 Calliope 版本。该框架提供了一个可组合的工作流管理系统，用于编排一系列软件模块以计算 EoS 和观测物理量。

关键模块与工作流

本文重点介绍了为重离子唯象学设计的模块，将之前的 10 个模块扩展到了 17 个：

• 强子相： 基于 Thermal-FIST 的 HRG（强子共振气体）模块，在完整的四维空间 $(T, \mu_B, \mu_S, \mu_Q)$ 中计算热力学性质。它支持理想、排除体积以及量子范德华（QvdW）相互作用。
• 格点 QCD 展开：
  • 4D-Taylor (BQS)： 使用连续极限估计的易受性，将压强按 $\hat{\mu}_X = \mu_X/T$ 的幂级数进行展开。
  • 4D-TExS： 将 $T'$-展开方案推广到四维，通过将笛卡尔化学势映射到球面坐标，以扩展其有效范围。
• 唯象模型：
  • Ising 2D-TExS： 结合二维格点 $T'$-展开与三维 Ising 临界标度，以模拟 $(T, \mu_B)$ 平面上的临界点和一阶线。
  • 全息理论 (NumRelHolo)： 使用爱因斯坦-麦克斯韦-狄拉顿（EMD）规范/引力对应模型来描述强耦合物质，该模型经过 LQCD 在 $\mu_B=0$ 处的校准，并扩展到有限密度。它包含了输运系数。
• 合成模块： 一种基于统计混合和稳定性判据的新程序，允许将两种不同的 EoS（例如全息理论和 HRG）合并为一个单一且热力学一致的 EoS，从而覆盖更广泛的 $(T, \mu_X)$ 范围。该程序可以产生交叉（crossover）转变、一阶转变或临界点。
• EoS 反转器： 一个数值反转模块，用于将 EoS 表从 $(T, \mu_X)$ 坐标转换为自然流体动力学变量 $(s, n_X)$ 或 $(\epsilon, n_X)$，这对于像 CCAKE 这样的拉格朗日流体动力学求解器至关重要。

模拟框架

作者使用 NuclearConfectionary 框架（具体包括用于初始条件的 AMPT、用于流体动力学的 CCAKE 2.0 以及用于强子后燃阶段的 SMASH）进行相对论粘性流体动力学模拟。

• 碰撞： 中心 Au+Au 碰撞，碰撞能量 $\sqrt{s_{NN}} = 7.7, 19.6,$ 和 $39$ GeV。
• 变量： 模拟追踪 SPH（光滑粒子流体动力学）粒子在相图中的路径，从而允许对动力学轨迹、热图和冻结分布进行分析。
• 输运系数： 切变粘度设定为常数 $\eta T/(\epsilon+P) = 0.08$。体粘度包含一个通过重子易受性 $\chi_B^2$ 模拟的临界标度成分。

核心贡献

1. 发布 MUSES Calliope： 本文引入了一套专门为重离子碰撞定制的新模块和工作流，补充了现有的中子星模块。
2. 扩展的 EoS 覆盖范围： 作者展示了构建能够覆盖 RHIC 能量扫描（BES-II）相关相空间的 EoS，包括 LQCD 展开通常失效的高 $\mu_B$ 和低 $T$ 区域。
3. 热力学合成： 合成（Synthesis） 模块的实现允许无缝合并差异巨大的 EoS 模型（如全息理论与 HRG），并具有可调的临界点，确保了跨相边界的热力学一致性。
4. 反转与稳定性： 测试了 EoS 反转器 处理 $(T, \mu_X)$ 与 $(s, n_X)$ 之间非线性映射的能力，包括处理一阶相变及亚稳态/不稳定区域的策略。
5. 动力学轨迹分析： 通过利用拉格朗日流体动力学求解器，作者绘制了流体单元在 QCD 相图中的特定轨迹，识别了不同碰撞能量和快度下探测到的区域。

结果

• EoS 有效性与覆盖范围：
  • 与标准的 4D-Taylor 展开相比，4D-TExS EoS 在 $(s, n_B)$ 平面上提供了显著更好的覆盖范围，特别是在较低温度和较高重子密度的情况下。
  • Holography+QvdW-HRG 合并 EoS 提供了最广泛的覆盖范围，有效地跨越了强子相和去禁闭相。
  • 在 $\sqrt{s_{NN}} = 19.6$ GeV 的模拟中，使用 4D-Taylor EoS 时，在中快度处仅有约 50-70% 的流体单元保持在有效 EoS 表格内，因此需要使用备份的共形 EoS。相比之下，4D-TExS 和 Holography+HRG 模型在 $\sqrt{s_{NN}} \ge 19.6$ GeV 时，在中快度处保持了 >90% 的覆盖率。
• 临界点敏感性：
  • 临界点的位置显著影响哪些碰撞能量对它的存在具有敏感性。
  • 对于位于 $\mu_B^c = 598$ MeV（与近期理论预测一致）的临界点，只有最低能量（$\sqrt{s_{NN}} = 7.7$ GeV）显示出显著比例的流体单元（$\sim$40%）经过临界区域。更高能量主要在较大的前向/后向快度处探测临界区域。
  • 在模拟中未观察到显著的“临界透镜效应”（等熵轨迹的聚集），这表明临界区域可能沿 $\mu_B$ 轴而非 $T$ 轴水平分布，或者非平衡效应抹平了信号。
• 一阶相变：
  • 包含跨越一阶相变的亚稳态和不稳定区域的模拟并未导致 CCAKE 求解器发生数值崩溃，因为负平方声速区域（$c_s^2 < 0$）极其微小且短暂（$\langle \Delta \tau \rangle \approx 1.4 \times 10^{-5}$ fm/c）。
  • 然而，作者指出，目前对一阶线的处理属于初步探索，并非完整的微观动力学相变理论。
• 观测物理量：
  • EoS 的选择（以及临界点位置）对总 $\pi$ 介子多重数（$dN/dy$）的影响极小，这表明归一化参数主导了这些观测物理量。
  • 研究发现，快度分布的宽度（$dN/d\eta$）强烈依赖于 AMPT 中的初始态弥散长度，而非 EoS 的细节。

意义与展望

本文建立了一个在统一、开源框架内系统研究致密 QCD 物质的基础设施。通过提供可合并且热力学一致的模块化 EoS，MUSES Calliope 引擎使得将重离子碰撞数据与理论模型进行更稳健的比较成为可能。

作者强调，其工作是“第一步”且属于“初步探索”，特别是在处理一阶相变和在输运系数中包含临界标度方面。他们明确指出，在动力学模拟中对一阶相变进行系统性处理需要进一步的理论和数值工作。其主要意义在于能够追踪流体在火球演化过程中采样了 QCD 相图的哪些区域，并量化基于展开式的 EoS 的有限有效性如何影响流体动力学预测。该框架旨在促进未来连接重离子碰撞与中子星约束的联合分析。