Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

技术摘要：重离子碰撞中重子数易受性（Susceptibilities）的贝叶斯提取

问题陈述
守恒电荷（特别是重子数）的涨落是探测量子色动力学（QCD）相图（尤其是假设的临界点 CP 附近）的敏感探针。虽然第一性原理格点 QCD 计算能够精确确定零重子化学势（$\mu_B = 0$）下的重子数易受性（$\chi_B^n$），但其在有限密度下的行为依赖于外推或其他理论方法。实验上，相对论重离子对撞机（RHIC）的束流能量扫描（BES）计划通过测量金-金（Au-Au）碰撞中净质子数的逐事件涨落来探测这一区域。然而，定量解释这些数据受到多种非平衡态和实验效应的影响：测量的质子仅携带总重子数的一小部分；实验接受度在动量空间受限；以及精确的全局重子数守恒非平凡地改变了相对于巨正则系综期望值的涨落。以往的策略通常是将数据与特定的模型计算（例如，强子共振气 HRG、带排除体积的流体力学或输运模型）进行对比，而非直接推断底层的易受性。

方法论
本研究展示了（据作者所知）首次直接从重离子碰撞数据中对重子数易受性（$\chi_B^2, \chi_B^3, \chi_B^4$）进行的贝叶斯提取。该框架构建了一个从理论易受性到测量的质子累积量（cumulants）的前向映射，包含四个截然不同的步骤：

1. 流体力学背景： 使用 MUSIC 代码生成事件平均的流体力学超曲面，采用 3 维初始条件，并在恒定能量密度（$\epsilon_{sw} = 0.26$ GeV/fm$^3$）下进行粒子化。
2. 最大熵冻结（MaxEnt）： 将超曲面上的局部重子易受性转换为重子-反重子联合累积量。该框架将局部易受性输入参数化为相对于理想强子共振气（iHRG）基准的比率：$\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$。MaxEnt 方案将流体力学场及其涨落映射到强子多重数不可约累积量上。
3. 接受度与滤波： 通过独立的二项采样，利用 Cooper-Frye 接受度和重子-质子滤波，将联合累积量映射到被接受的质子和反质子，从而考虑 STAR 探测器的运动学接受度。
4. 精确守恒： 子系综接受法（SAM-3.0）施加精确的全局重子数守恒，将巨正则累积量转换为适用于与实验数据比较的正则累积量。

贝esian 推断是在 $\sqrt{s_{NN}} = 7.7–27$ GeV（具体为 7.7, 9.2, 11.5, 14.6, 17.3, 19.6, 和 27 GeV）的 BES-II 范围内的 0–5% 中心度 Au-Au 碰撞的 STAR 数据上进行的。研究使用了两种互补的提取方式：(i) 质子阶乘累积量比值（$\hat{C}_n/\hat{C}_1$）和 (ii) 净质子累积量比值（$\kappa_n/\kappa_1$）。对易受性比值应用了平坦先验，并使用统计和系统误差构建了高斯似然函数。

主要结果

• 二阶易受性 ($\chi_B^2$)： 分析对归一化后的二阶易受性（$\chi_B^2/\bar{\chi}_B^2$）给出了紧凑的约束，相对不确定度为 5–10%。提取的值在测量的能量范围内表现出非单调的碰撞能量依赖性：在最低能量（$\sqrt{s_{NN}} = 7.7$ GeV）处出现增强（$\sim 1.06–1.07$）；在中间能量（$11.5–19.6$ GeV）处出现极小值（$\sim 0.73–0.82$）；并在较高能量处表现为平坦的抑制（$\sim 0.8–0.9$）。
• 与格点 QCD 的比较： 当映射到化学冻结线时，提取的 $\chi_B^2$ 值与格点 QCD 基于的（特别是 4D-TExS 状态方程）估计值在 $\mu_B \lesssim 300$ MeV 时具有定量一致性。然而，在较大的 $\mu_B$（对应于 $\sqrt{s_{NN}} = 7.7$ GeV）处，提取的值偏离了格点 QCD 带，显示出高达 $\sim 50\%$ 的增强。
• 高阶易受性： 由于从重子到质子的映射效率损失以及来自重子守恒的畸变，三阶和四阶易受性（$\chi_B^3, \chi_B^4$）仅受到较弱的约束。虽然四阶比值在 $\sqrt{s_{NN}} = 19.6$ GeV 处显示出对负值的局部偏好（$\chi_B^4/\chi_B^2 = -33 \pm 14$），但巨大的不确定性使得无法仅凭这些高阶矩对 QCD 相结构得出确定性结论。
• 相关性主导： 一个假设仅存在不可约二重重子相关性的极简方案（设置 $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$）成功再现了在 $\sqrt{s_{NN}} \approx 11$ GeV 附近观察到的质子阶乘累积量比值 $\hat{C}_3/\hat{C}_1$ 的非单调峰值。这表明观察到的结构源于 $\chi_B^2$ 的能量依赖性与精确重子数守恒之间的相互作用，而非需要不可约的三重或四重重子相关性。

意义与主张
本文声称提供了首次基于数据的、贝叶斯确定的沿化学冻结线的重子数易受性，避免了对特定状态方程（EoS）的假设。作者断言，其框架成功分离了二阶易受性，发现其在较低重子密度下与格点 QCD 一致，但在较高密度（$\mu_B \gtrsim 300$ MeV）下，相对于非临界外推显示出潜在的重子数涨落增强。

研究谨慎地总结道，虽然观测到的现象可以在不存在高于二阶的不可约多重重子相关性的情况下得到描述，但这并不排除它们的存在；相反，目前的测量主要对二重重子贡献敏感。作者指出，在低碰撞能量下观察到的增强值得通过未来的固定靶实验进行审查，并且目前的提取是以固定的流体力学和冻结设置为条件的，这代表了全局分析中不确定性的下限估计。该框架被呈现为一个用于约束有限密度 QCD 状态方程的通用工具，并且一旦建模不确定性得到更好控制，可以扩展到包括固定靶能量、体积涨落修正以及其他守恒电荷。