Out-of-equilibrium contributions to charm hadrons in a fluid-dynamic approach

技术摘要：流体动力学方法中魅强子（Charm Hadrons）的非平衡贡献

问题陈述
由于魅夸克（Charm quarks）具有较大的质量，确保了它们主要在初始硬散射中产生，且在火球寿命内无法通过热生产达到热平衡，因此魅夸克是研究夸克-胶子等离子体（QGP）的独特探针。先前的流体动力学框架通过假设其接近局部热平衡来成功模拟魅夸克的动力学过程。然而，这些模型忽略了两个关键阶段的非平衡修正：碰撞初始时刻（自由流阶段之后）以及冻结面（Freeze-out hypersurface）。忽略冻结修正导致魅强子多重数对空间扩散系数（$D_s$）产生了非物理的依赖性，从而阻碍了对魅输运系数的系统性确定。此外，初始分布通常被假设为处于局部热平衡状态，忽略了初始魅生产的非平衡性质。

研究方法
作者扩展了先前开发的用于描述 QGP 中魅夸克的流体动力学框架，通过显式计算初始阶段和冻结表面的非平衡修正。

1. 初始阶段初始化：

   • 该模型在魅生产时间与流体动力学初始化时间（$\tau_0$）之间引入了一个自由流阶段（Free-streaming phase）。
   • 利用特征线法（Method of characteristics），作者解析地求解了粒子在这一阶段沿直线轨迹运动的动力学方程。
   • 这种方法修改了初始魅密度，并在 $\tau_0$ 时刻产生了一个非零的扩散电流（$\nu^\mu_q$），该电流源于二元碰撞密度（$n_{coll}$）的空间依赖性和初始硬散射的动量分布。
   • 采用兰道匹配条件（Landau matching conditions），根据微分产生截面来确定总魅分布函数的积分。

2. 冻结修正（多物种方法）：

   • 作者在魅强子气体（HRGc）中采用了多物种方法，将系统视为由所有魅强子物种与轻组分相互作用组成的流体。
   • 利用矩方法（Method of moments）和纳维-斯托克斯近似（Navier-Stokes approximation），推导了每个强子物种 $i$ 的分布函数非平衡修正（$\delta f_{i,p}$）。
   • 该修正被展开为不可约矩，并截断在秩 $l=1$（扩散电流），假设魅强子的体压和剪切应力贡献相对于轻组分可以忽略不计。
   • 一个关键的解析结果（公式 23）将修正表示为：
     $$\delta f_{i,p} = f^{(0)}_{i,p} \frac{1}{P_i} \bar{\kappa}_i \nu^\mu_q p_{\langle \mu \rangle}$$
     其中 $\bar{\kappa}_i$ 是取决于物种电荷和密度的归一化扩散系数，而 $\nu^\mu_q$ 是总魅扩散电流。该公式允许仅演化总电流而非单个物种的电流。

3. 演化与谱度：

   • 系统使用 1+1D 流体动力学框架（FluiduM）进行演化，该框架具有恒定的 $D_s T$。
   • 通过 Cooper-Frye 程序计算不变产额和横动量（$p_T$）谱，并纳入推导出的非平衡修正。
   • 模型使用 FONLL pQCD 预测和 ALICE 测量的魅截面作为输入，针对实验数据（ALICE 在 $\sqrt{s_{NN}} = 5.02$ TeV 下的 Pb-Pb 碰撞，0–10% 中心度）进行了测试。

核心贡献与结果

• 解决非物理依赖性： 在冻结阶段引入非平衡修正，解决了魅强子多重数对空间扩散系数 $D_s$ 的非物理依赖问题。若没有这些修正，多重数会随 $D_s$ 显著变化；而加入修正后，在测试的所有 $D_s$ 范围内，其积分产额与无扩散情况保持一致，且与初始魅夸克总数相符。
• 初始扩散电流： 自由流初始化产生的扩散电流在 $\tau_0$ 时刻比后期产生的电流大一个数量级。然而，研究发现，与零初始电流的情景相比，这种初始条件对后续魅场演化及最终不变产额的影响微乎其微。
• 模型与数据对比：
  • 使用 FONLL 截面时，模型在不确定度范围内重现了魅介子的积分产额。
  • 在 $\Lambda_c^+$ 重子产额上观察到偏差，该产额被低估。遵循既往文献，作者通过将激发态魅重子的统计权重增加三倍，应用了“魅重子增强”（Charm-baryon enhancement）。经过此调整并使用 ALICE 测量的截面后，模型在 $p_T \approx 3$ GeV/c 以下实现了与实验产额及 $p_T$ 谱的兼容。
• 适用范围限制： 论文指出了基于分布函数正定性的流体动力学方法的根本限制。
  • 由于扩散电流的径向剖面为负，非平衡修正 $\delta f_{i,p}$ 可能变为负值。
  • 当 $|\delta f_{i,p}| \approx f^{(0)}_{i,p}$ 时，总分布函数变为负值，导致物理阐释失效。
  • 这种 $p_T$ 谱的“零交叉”现象随着 $D_s T$ 的增加而在更低的 $p_T$ 值处出现。对于 $D_s T \gtrsim 0.2$，模型在 $p_T \sim 2$ GeV/c 处会出现负谱。
  • 因此，该框架仅在 $D_s T \lesssim 0.1$ 时在整个 $p_T$ 范围内可靠；或者对于较高的 $D_s T$，仅在受限的低 $p_T$ 区域（例如当 $D_s T \approx 0.24$ 时，$p_T \lesssim 2.5$ GeV/c）内有效。

重要性
本文建立了一个一致的魅强子流体动力学框架，考虑了碰撞初始阶段和最终阶段的非平衡效应。通过修正产额对输运系数的非物理依赖，这项工作为未来系统性提取魅扩散系数（例如通过贝叶斯分析）奠定了必要基础。然而，作者谦虚地总结道，当前的线性化方法受限于摄动假设（$|\delta f| \ll f^{(0)}$）在高扩散系数和高横动量下的失效。他们建议，未来的工作必须探索替代方法，例如最大熵方法，以确保在整个运动学范围内分布函数均为正值。