Nonconformal kinetic theory with transverse dynamics and event-by-event simulations using CoMBolt-RTA

技术摘要：基于 CoMBolt-RTA 的非共形动力学与横向动力学及逐事件模拟

问题陈述
高能核碰撞中热 QCD 物质（夸克-胶子等离子体）的产生传统上使用相对论流体力学进行建模，该方法假设系统处于局部平衡态。然而，小系统（pp, p-Pb）中集体效应信号的观测以及关于流体力学适用性的持续争论，使得寻找替代方法变得必要。虽然存在输运模型，但它们通常依赖于唯象截面或微扰 QCD 散射。此外，标准的动力学方法通常假设准粒子是无质量的，这会导致共形状态方程（EOS），从而偏离格点 QCD 的结果，特别是在声速和相变区域的表现。因此，需要一种能够包含有质量准粒子以实现现实 EOS、包含横向空间动力学，并能架起自由流（free-streaming）与理想流体力学之间桥梁的动力学框架，用以研究从流体行为到粒子行为集体的转变过程。

方法论
本文介绍了 CoMBolt-RTA，这是 CoMBolt 事件生成器的一个新版本，旨在求解具有质量准粒子的玻尔兹曼方程，并采用弛豫时间近似（RTA）。

• 理论框架： 该模型采用了 (2+3)D 玻尔兹曼方程，其中系统具有两个空间维度（纵向不变）和三个动量维度。准粒子质量 ($m$) 被视为一个常数，这区别于之前的无质量版本 (CoMBolt-ITA)。碰撞核定义为 $C[f] = -\frac{p_\mu u^\mu}{\tau_{relax}}(f - f_{eq})$。
• 状态方程 (EOS)： 通过调节简并因子 ($g_B$) 和常数准粒子质量 ($m$)，该模型生成的 EOS 可以逼近格点 QCD 结果（具体为 HotQCD）。文中测试了三组参数集（EOS1, EOS2, EOS3），其中 EOS3 在部分子-强子切换温度（$T_{sw} = 155$ MeV）附近与格点数据最为契合。
• 耗散效应： 弛豫时间 $\tau_{relax}$ 通过函数 $\hat{H}(T_{eff})$ 进行参数化，该函数控制特定的剪切粘度 ($\eta/s$) 和体粘度 ($\zeta/s$)。模型允许使用温度依赖的 $\hat{H}$，从而能够重现非零的体粘度，而体粘度在共形（无质量）极限下会消失。
• 数值实现： 玻尔兹曼方程在 (2+3)D 网格上使用特征线法进行求解。动量空间利用 Gauss-Laguerre 求积进行积分，避免了对动量空间的插值需求，提高了稳定性。
• 混合模拟链： 对于逐事件模拟，代码采用混合方法：
  1. 初始态： 使用 TRENTo 生成（针对氧、氦-3 和 氦-4 的特定核构型）。
  2. 部分子演化： 通过 CoMBolt-RTA 进行求解。
  3. 强子后燃器（Afterburner）： 系统在 $T_{sw}$ 时切换至强子气体，并使用 UrQMD 进行演化。

核心贡献

1. 有质量准粒子动力学： 主要贡献是将 CoMBolt 框架扩展到包含有质量准粒子。这使得非共形 EOS 和非零体粘度的产生成为可能，解决了以往无质量模型的局限性。
2. 横向空间动力学： 模型明确包含了横向空间扩展，允许研究空间各向异性如何在非共形介质中演化为动量各向异性。
3. 不透明度与定标： 研究了“不透明度”参数 ($\hat{\gamma}$)，该参数定义为以微观弛豫长度为单位的系统尺寸，用于量化流体类与粒子类机制之间的转变。
4. 轻离子碰撞预测： 模型应用于氧-氧（OO）碰撞，并扩展到预测轻系统（$^3$He–$^3$He, $^4$He–$^4$He）以及 p-Pb 碰撞的可观测物理量。

结果

• 平衡演化： 对高斯初始态的模拟表明，与无质量情况相比，有质量准粒子表现出更丰富的各向异性动力学（涉及纵向、横向和体反向雷诺数）。更接近自由流（大 $\hat{H}$）的系统在较冷区域会产生更多各向异性，而更接近理想流体力学（小 $\hat{H}$）的系统则平衡得更快。
• 与 OO 碰撞数据的比较：
  • 各向异性流与多重数： 通过手动参数调优（具体为 EOS3 配合 $\hat{H}$ param1），模型较好地描述了 ALICE 关于 OO 碰撞中的椭圆流 ($v_2$)、三角形流 ($v_3$) 以及带电多重数的中心度依赖性，尽管在高中心度下表现有所下降。
  • 能谱与 $\langle p_T \rangle$： 与无质量版本相比，引入准粒子质量显著改善了对横向动量 ($p_T$) 能谱和平均横向动量 ($\langle p_T \rangle$) 的描述。新版本减少了 CoMBolt-ITA 中存在的对 $\langle p_T \rangle$ 的系统性高估，使预测更符合 ATLAS 数据。这一改进归功于由有质量 EOS 导出的更真实的声速。
• 轻系统与不透明度：
  • 定标： 随着不透明度的增加，不同碰撞系统（OO, $^3$He, $^4$He）之间的集体响应（动量各向异性相对于初始空间各向异性）表现出近似的定标关系。
  • 流体 vs 粒子机制： 中心 OO 碰撞仍处于流体类机制中。然而，模型预测中心 $^3$He–$^3$He 和 $^4$He–$^4$He 碰撞处于流体类与粒子类动力学的过渡区域。
  • p-Pb 差异： 模型低估了 p-Pb 碰撞中的各向异性流，这表明在建模质子初始态或极小系统中的集体性描述方面可能存在不准确性或缺失物理机制。

意义与主张
论文声称，CoMBolt-RTA 通过引入有质量准粒子，提供了更真实的 QGP 演化描述，这对于匹配格点 QCD 热力学和实验观测量（如 $\langle p_T \rangle$）至关重要。作者断言，虽然中心 OO 碰撞可以被很好地描述为流体，但像 $^4$He–$^4$He 这样更小的系统，由于其较小的尺寸和已知的核结构，为探测向粒子类集体动力学转变提供了独特的机遇。

作者对其结论保持谨慎，指出目前的结果依赖于手动参数调优而非完整的贝叶斯分析。他们明确指出，需要系统的贝叶斯分析才能对粒子类机制的精确转变范围得出稳健结论，并充分约束模型参数。此外，他们也承认目前的恒定质量假设限制了完全重现温度依赖 QCD EOS 的能力，而这正是未来拟实现的特性。